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on the telemetry under S. A. Dabrowski, the analog simu- 
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Appendices B and C were written by H. R 0 Cowan and 
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SUMMARY 


This report presents mathematical models of rate gyros, servo-accelerometers, 
pressure-transducers and an FM/FM telemetry system used with Scout 
vehicles. These models permit prediction of instrument response to transient 
signals, more particularly to signals which are of a duration short compared 
with normal instrument response time and which reach magnitudes large com- 
pared with the nominal range (i. e. , imperfect impulses). The models are 
derived theoretically for all components; and, whenever possible, their relevant 
parameter values are calculated from data obtained by simple measurements. 
The instruments were tested in the laboratory and afterwards analog computer 
simulations were run to find the parameter values giving the best fit to the 
laboratory results. The theoretical values and those obtained on the simulation 
were then compared. Generally, a good agreement was found to exist. 

A method is derived to reconstruct input signals to the transducers, given the 
recorded output signals, and taking into account the effect of noise. A digital 
computer program, based on this method, was prepared and is described here. 
This inversion scheme is very useful in post-flight analyses to obtain a better 
idea of the true shape of the input signals to the transducer, since for fast 
signals the transducers introduce a considerable dynamic lag and distortion. 

The inversion method developed here is optimum in the minimax sense of 
minimizing the maximum instantaneous reconstruction error. 


n 
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1. INTRODUCTION 


In the Scout vehicles, and more generally in aircraft and space vehicles, 
certain events can cause the sensors and transducers in the flight instrument 
package to receive signals which attain unusually large instantaneous magnitudes 
and which vary noticeably in a time interval that is short compared with the 
normal instrument response time. Catastrophic failures in some components, 
such as the combustion chamber, are likely candidates to generate these signals. 

With such signals it becomes important to know the transient response 
characteristics of the transducers and the telemetry system, in order to 
evaluate dynamic effects (such as lag, delay, distortion) introduced by them. 
Therefore, the first part of this report is concerned with the mathematical 
description of the transient response characteristics of three categories of 
instruments (rate gyros, servo accelerometers, and pressure transducers) and 
of an FM/FM telemetry system. In each case, first a theoretical model is con- 
structed, based on whatever is known about the inner structure of the devices. 

The numerical values of the relevant parameters are then calculated, as far as 
possible, from simple measurements of dimensions, weights, and data supplied 
by the manufacturer. Useful information was obtained by opening the instruments 
(after other laboratory tests were finished) and by inspecting parts and running 
tests of transient response to disturbances applied at interior points. 

Then, the instruments and telemetry components were subjected to tests in the 
laboratory. A large variety of input signals were used, depending on the equip- 
ment employed for particular tests. More particularly, the signals mentioned 
before were defined (in the contract) as signals with amplitudes up to ten times 
the nominal full scale range and occurring in time intervals down to two orders 
of magnitude less than the normal instrument response time. Such signals 
approach impulses and with the definition just given the equivalent impulse 
strength (which is the area under the input signal or duration times average 
magnitude) is sufficiently limited to keep the instrument response in its linear 
range. Therefore, linear models could be used for most components. Even for 
a non-linear instrument as the rate gyro, linearized models were generally 
sufficient. In some cases, such as the phase sensitive demodulators and some 
accelerometers and pressure transducers, a linear model is applicable, 
provided some parameter (a damping coefficient) is made dependent upon the type 
of input signal (polarity or speed, respectively). 

Finally, analog computer simulations were made, using the theoretical model 
to find the parameter values giving the best agreement with the laboratory test 
results. These values could then be compared with the ones derived from theory 
earlier. In general the two sets of parameter values were quite compatible. 
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Once models are constructed, it is relatively easy, using analog or digital 
computers, to determine the response of the elements studied, given any 
admissible input signal. However, in many practical applications, such as 
post mortem analyses of in-flight anomalies, it is rather the reverse problem 
that matters: one has recorded output signals coming from transducers and 

telemetry channels and one wants to know what were the true input signals as 
seen by the transducers. For this reason, the second part of this report develops 
the theory of this inversion problem for the three kinds of transducers. A 
digital computer program, based on this theory, prints and plots automatically 
the reconstructed transducer input, given the output signal. This reconstruction 
is, in fact, an optimal one in the minimax sense that the method used produces 
the smallest achievable worst error at any time (this is similar to Chebychev 
approximation, Reference 1)# Since the flight instruments act as signal 
smoothers, the inversion is fundamentally a differentiation problem in the presence 
of noise, and the inversion program is some sort of optimum filter. 

Most of the references listed in the bibliography deal with problems akin to the 
ones treated in the present report. For example, the modeling of gyros can 
be found in many other articles or books than the one quoted in the report 
(Reference 2). On the other hand, the resonance calculation of the Bourdon 
tubes, certain parts of the discussion of the band pass filter response, the 
discussion of interference ripple are, as far as is known, not found elsewhere. 

A particularity of this study is the systematic comparison of values obtained 
from the theory and those obtained through simulation. 

While the use of a minimax criterion in filtering can be considered classical, 
the inversion formulation in this report is new. It can, moreover, be generalized 
to allow inversion of a wide class of linear systems. 

For readers who are not interested in the mathematical details of the modeling, 
Section 7 (Conclusions) gives an expanded verbal description of the most 
important results. 

Note: Time scales on the figures are given in ms/Division. This applies to 

"major" divisions, corresponding with a distance of 1 cm on the oscilloscope 
screen. 


/ 


4 



2. RATE GYROS 


The rate gyros in the Langley Scout package belong to one of two types: 

for pitch and yaw rates: Honeywell GN 90B1 

for roll rates: Honeywell GN 91B1 

Both types are close to the Honeywell GG440 GNAT miniature design. Some of 
their characteristics are specified in the procurement specifications, an extract 
of which is listed below: 


Pitch and Yaw Roll 


Full Scale Rate (deg/ sec) 

40 

100 

Natural Frequency (Hz) 

19-26 

36-42 

Damping Ratio 

0. 49-1. 4 

i—l 

1 

O' 

o 

Scale Factor V rms/ (deg/ sec) 

0. 235 

0. 095 

(with correct load) 

(+ 5%) 

(+ 7%) 


2. 1 Theoretical Model 


2.1.1 Mechanical Part of the Gyros 


The model is based on Newton 1 s angular momentum equation 


= T 
dt 


(2-1) 


where H is angular momentum, T is torque, and which is valid around either 
an inertially fixed point or the center of mass (C. M) of the system under con- 
sideration. For a moving frame F with instantaneous inertial rotation vector 
£3j:, (2-1) is equivalent with: 


(— ) 


o t 


+ u) f x H = T 


(2-2) 


in frame F 


The following development parallels closely the one given by J. M. Slater and 
J. S. Ausman in Reference 2. 


The figure below shows the basic configuration. 
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Frames: (x, y, z) is the gimbal frame with x = input axis, 

y = output axis, z = spin reference axis. 


(x c , y c , z c ) is the case frame 
(x m , Ym» z m) is the mounting frame 
(x 0 , y Q > z 0 ) is the inertial frame 

Rotations: 0 is the rotation angle of the gimbal with respect to the case, 

along the output axis. 

0p is the rotation angle of the gimbal with respect to the case 
along the input axis. 

0 and 0p are small angles and hence are commutative. 

0 C is the rotation of the case with respect to the mounting, 
along the input axis. 

0 Z , 0y. 0x are the Euler angles, in that order, describing 
the rotation of the mounting with respect to the inertial 
frame. (See Figure below). 
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The general transformation matrix from (1 XQ , ,Vp’ lzo) to ^ 1 xm’ Vm> W is 
standard (see Reference 2 again) and is listed below. 


ixm; [ c °s 0 y cos 0 


cos0, sin0 

y z 


-sin0. 


J ■'•y rr ^| = | sin0 x sin0ycos0 z -cos0 x sin0 z sin0 x sin0ysin0 z +cos0 x cos0 2 sin0 x cos0y 
^zmi |cos0 x sin0yCos0 z +sin0 x sin0 z cos 0 x sin0ySin0 z - sin0 x cos0 z cos0 x cos0y 


xo 


yo 


The corresponding transformation matrix for small rotation angles is: 


i Ixm ( 


1 ‘?y 1 

f 

i Ixo ] 

t iy m * 

1 

■K \ 

< 

1 •- 
0 

(. 1 zm . 

) 

*y 'K 1 


l ^o ) 


The other transformation matrices are analogous: 


IK 

_ 

1 

0 

0 

1 

- 1 
9^ 

S 

In' 

I 

9 

- 6 P 



jlxc 


1 

0 

0 

) 

1 1 „ 

> 

0 

1 


< 

i y 
V. i zc 


0 

“ e c 

c 

1 

1 1 
1 l 



The angles 9, 9p and 9 C are struaturally constrained to small values. The 
angles 0 X , 0 y and 0 Z are small in the tests described here. Typically, test 
inputs are pulses of amplitude not exceeding ten times the nominal value, and 
duration at least one hundredth of the basic time constant of the system. 
Therefore, for the roll gyro one can expect an angular rotation of not much 
more than 
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A 0 X = (1, 000 °/ sec) • (0. 75 x 10 -3 sec) = 0. 75° 


which is indeed small enough. 

The different instantaneous rotation vectors are calculated next, up to and 
including second order quantities. 

The rotation vector of the mounting with respect to inertial space is 


lU 


(K - K K) i. 


+ (0y + ^ym 

• ♦ 

+ (0 - 0 0 ) 1~ 

' z x y zm 

The rotation vector of the case with respect to inertial space is 
^ c = + 0 C l xc (which becomes, with the aid of the matrices above) 


= (»C + *X - *> y *>z> ^c 

+ [ 0 y + 0 z (0 x + e c )J l yc 

+ [ 0 2 ^y ($ x + ^ZC 


The rotation vector of the gimbal with respect to inertial space is i/J , (the 
sequence 9p“^ 9 is chosen here, but this choice does not influence the final 
formulas). 


g = [L c + 9p l x + 9 ly + 9 9p 1^ = (use again the transformation matrices) 

= 1 9 p + e c + 0 X -b z (0 y + 9)] T x 

+ 1. 8 + by + K («x + e c + V ] V 

+ L K - K ( 0 p + 6 c + f> x ) + 6 (* e c + 0 p + M T z 
Output axis (or gimbal) equation. 
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For a gimbal with well-oriented principal axes, carrying a rotor with relative 
angular momentum H r , the total angular momentum with respect to the C. M. 
of this subsystem is given by: 


H g = X gx 


oO 


g x 


1 + I 

x gy gy 


i + 

y 


(I gZ c ^gz + 


H r) 1, 


where I is the sum of gimbal z-axis inertia plus wheel -inertia around the 
spin axis; Ig Z and Igy are gimbal -plus -wheel inertias along input and output 
axes. 


All inertias are calculated along axes centered on the center -of -mass of gimbal 
and wheel. If one writes the output axis component of (2-2) and retains only 
those second order quantities which are multiplied by the large rotor momentum 
H r , the following output -equation is obtained. 

■gy •' + B gy ® + K gy 6 = H r * 8 P + K' 

l gy 0 y ^r ^y + + ^gy, residual 

- (eu ew S,(^ A x - ^ \) ( 2 - 3 ) 

The meaning of the symbols is: 


B 

K 

T 


g: 

gy; 

gy : 

gy> 


Mgw‘ 

^ z^x : 

S : 

A X> -^z : 


■p 

reference acceleration, e. g. , g = earth gravity = 32. 17 ft/ sec 
output-damping coefficient 
output spring constant 

residual* ot ^ ier (p a **asitic) torques along output axis, exerted on 
gimbal (such as gravitational and magnetic torques) 
combined mass of gimbal plus wheel 
C. M. offsets away from intersection of xyz-axes. 
maximum C. M. offset ( 1 S y ^ i o z ] 4 cT ) 

components of inertial linear acceleration of the gimbal and wheel 
system, expressed in units of M g" 


Input Axis Equation 

C. M. offsets are neglected here. First, there is the x-axis equation for motion 
between gimbal and case: 

’« * % + V 6 P + K g* e p = - H r <« + V 
- 'ax <«'c + <0 - H r <>z + 8 P> 


g X 
+ T 


gx, residual 


(2-4) 
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where all the notations are analogous to these used in the output-axis equation. 
Secondly, there is the x-axis equation for motion between case and mounting. 


I 0 + B 0 + K 0 = - I 0 

cx c cx c cx c cx x 


+ B gx e p + K gx e p " ( P art of T gx, residual due gimbal-case interaction 

only) 


+ T 


mounting on case. 


x, 


residual 


(2-5) 


with notations analogous to these used above. 


Special cases for Input-Axis Equations: 

A) If the case is rigidly attached to the mounting, K c °° 

0 C 0 and only (2-4) is left with 0 C = 0, while (2-5) is to be omitted. 


B) If the gimbal-case connection is rigid along the input axis, K 

9p 0 and only one equation is left, by addition of (2-4) and^(2-5). 

dcx + Igx) ®c + B cx + K cx e c = - Gc* + : gx) *’* 

- H r (0 + \) - H r 'fi z (f> x + e c ) 

+ T mounting on case, x, residual + <P art of T yx not due to case > 


C) Finally, if input-axis motions are negligible, (2-4) and (2-5) are both 
omitted and 0 C = 9p = 0 is used in (2-3). 

Canonical Form of Equations: 


The output axis equation (2-3) can be rewritten as: 

H 

i 

r 

gy gy 


9 


v a / H \ ’ P 

+ 2 (gy !- t8= (^) <‘* t8 P + e c»- 

gy gy ' w gy 


H 


K 


gy 


0 Z (0 + 9) + T gy. residual J V 

Z y K gy 1 H r J[K gy J \? 


s. 


2L A 

£ 

(2-6) 


Z ) 
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with cO gy = K gy /Ig y 


_S2_ 


gy 2 (I K . 1/2 

2 gy K gy ) 


The input axis equation (2-4) can similarly be rewritten as; 


!z. ^E_ 

“’gx 2 ’^gx 


H 


18 + V- 


© C + fi 


gx 


oJ 


gx 


H 


r 0 (0 +9 + 9 ) + gx, residual 

K z x c p K 

gx gx 


(2-7) 


with 


^ gx = V ' 'gx 


gx 


gx 


2 (I gx K gx) 


1/2 


Stops: There is a mechanical bound on 0: 

[ 9 U 9 


max 


(2-8) 


Note on the Relationship between Input-Axis Stiffness K and Output Axis 


Stiffness K 


gy 


gx 


Kgy results from the torsional stiffness of a bearing member, while K depends 

on the bending stiffness of the same member. The gimbal is pivoted at one end, 

which gives a lever -arm amplification for K (See Figure below). 

gx 


Pivot Gimbal-Axis 



[ 


Elastic Part of Bearing 


I 

I 



n 


The elastic member is, at least in the main part, a circular cylindrical rod. 
Thus, one gets, using classical results of strength of materials (see e. g. , 
Reference 3) 


g* = 2 Ll . 

gy L ‘ 2 


E 

G 



> 21. 3 


since for most metals the ratio E/G of Young* s modulus to torsion modulus equals 
2.6, 8.2, (I^/Ip) = 1/2 (bending area moment of inertia divided by 

polar area moment of inertia for a circular cylinder). 


The true value of the ratio Kg X /Kgy turned out to be much larger than the lower 
bound derived here. 


Some Numerical Results 


From manufacturer's data one can estimate 
M S A 

| E — £> w |^5. 3 x 10” (rad/ sec) per g 

r 

further ® m ax = 3° = 0. 052 

Steady state gain H 

— — = 0. 0647 sec (pitch gyro) 

K 

gy 

pr 

£- = 0. 0248 sec (roll gyro) 

K gy 

Nothing is known about T gy, residual and Tg X , residual. 

One may assume the case-mounting connection to be rigid, i. e. , K c = 00 
and 9 C = 0. 
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2. 1. 2 Electro -Mechanical Output Device 


A Microsyn (variable reluctance) pickoff is used to generate the electrical 
output signal. This device can be viewed approximately as a transformer with 
output -axis -dependent (i. e. , 9 -dependent) transformation ratio nj^. 

Therefore, one obtains the circuit diagram of the figure below. 


Ideal Transformer 



Model for the Microsyn 

This model is now somewhat simplified by shifting the magnetization reactance, 
and the diagram of the figure below is obtained b y adding the load impedance 

Z L = ?L 

1 + R l C l s 



Simplified Model of the Gyro Output Circuit 
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The elements of Z Qut on that figure are given by: 

R o = R 2 + n M R 1 

2 

L o - L 2 + "M L 1 

T 1 _ n 2 
L o “ n M L M 

Strictly speaking, these parameters depend on 0 through the ratio nj^j. If this 
dependence is neglected, or an average value for R Q , L q and L 1 is adopted, 
the electrical output equations are: 


‘/z. 


Uj(t) 

= V s 2 cos 

(w s t + <p g ) 

(2-9) 

V i„ W = C M 8 W 

V = 

u x (t) 

(2-10) 


dt 


(2-11) 

V. = 
m 

V T + R i T + L„ di L 
L ° L ° dt 

(2-12) 

*l = 

+ V L 

+ C T dV L 

(2-13) 

L'o r L 

L dt 


Numerical Results for the Electro-Mechanical Output Device 

Steady state (i. e. , constant input rate) experiments on the pitch and yaw rate 
gyros show that for 

cO = 2 TX 400 = 2, 513 rad/ sec 

v g = 6. 68 V rms 

R l = 8. 06 X 10 8 ohm; G L = = 1.241 x 10" 4 mho 

C L = 0. 28 x 10" 6 F; B L = c0 g C L = 7. 04 x 10" 4 mho 
v-j^ and v are in phase and that the scale factor is: 

— 0 25 ^ ® _ i ao V r m s 

de 

zero phaseshift 


g/sec 




rad/ sec 


(°-Ik 

^ 7 ) 9 


) 


14 



Moreover, for a purely resistive load (C^ = 0), v L leads v in by 50° and the 
scale factor is: 


C 


d V L 

3 b 


) 

Resistive Load 


0 . 18 


V rms 
deg/ sec 


10 . 3 


V rms 
rad/ sec 


Setting 


x„ = X'„ = b t . = c Tj ; g t . = 


s ^ o’ U L 


■'Ll* L 




1 =Y l =G t +j(B L - i— ) 



one has 


V L 



1 

1 + z o y L 


with 

Z o Y L = R o G L + X o - B l) + j tG L X o - R o <57 - B l)] 

The output-impedance is given by 

7 = (H 0 + jX 0 )jX' o = R 0 X ' 0 2 + jX- Q [ r/ +X o (X Q +X' 0 )| 

out 7 7 

R 0 + j< X o + X> 0 ) K q +(X - 0 + X o ) 


According to the manufacturer, = (600 + j 1, 000) ohm. Use of this value 

of Z ou t, together with the fact that there is zero phase shift for and R^ as 
given above leads to 

R q = 2, 260 ohm 

= 72. 8 ohm; L = 0. 029 henry 
o o 

l - 4 

X ' 0 = 1, 413 ohm; -^y- = 7. 08 x 10 mho; L ! o = 0. 562 henry 


As a verification 


4- (Z Q Y l ) = 0°; (Z Q Y l ) 


0.281 


Z . = 616 + j 1, 008 
out J 


For the pureLy resistive load, Vl leads V^ n by 50°, agreeing with the tests. 
The proportionality factor Cjyj is now given by 


C M 



< 


< V- > < "1 > 


where the symbol < • /> is used for steady state amplitude ratios in the zero 
phase shift case. 

Then 


C 


M 


< v T 


< 


V in 


/ #x> 

> V s 


volts rms 
rad/ sec 


H 


(r-) 


gy 


This gives 

(Cm) = 9 rad’* = 0. 80 deg’* 

pitch 
yaw 

(see also experimental results later) 


Assuming that the same microsyn is used for the roll gyro, one has 


(C M ) 

roll 


(%) 


pitch, 


yaw 


As a verification, one obtains the following scale factors for the specified 
load of Rj^ = 10^ ohm; = 0. 3 x 10“^ farad, nominal conditions. 

V rms 

Scale Factor Pitch-Yaw = 0. 244 — - 

deg/ sec 


Scale Factor Roll 


0. 093 


Y rms 
deg/ sec 
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The specifications give 0. 235 and 0. 095 — (+ 5%) 

deg/sec deg/ sec — 

respectively, showing that the obtained values are within the prescribed 
tolerances. 

The phase shift is found to be 5.2°, less than 10° which is acceptable. 

Summarizing, the nominal values for the circuit are: 

V s = 6. 3 V (rms) 

= 2 TT 400 = 2, 513 rad/ sec 
= arbitrary phase 
R 0 = 2, 260 ohm 
L Q = 0. 029 henry 
L l 0 = 0. 562 henry 
= 10^ ohm 

C L = o. 3 X 10' 6 farad 

Cm = 45. 9 (rad) -1 = 0. 80 (deg) -1 



2.1.3 The Gyro Mounting Bracket 


The three rate gyros are assembled on a, mounting block as shown in the figure 
below. The frame (x m , y m , and z m ) of section 2. 1. 1 is attached to this 
mounting block. The block is then firmly fastened to a mounting bracket of 
stainless steel which in its turn is fastened to the rocket frame. 



Mounting block and rocket frame are to be considered as rigid when compared 
with the bracket. Therefore, an elastic analysis of the bracket, loaded by 
mounting block and gyros is in order, to determine the resonances. The rather 
complicated structure of the bracket is replaced by a simple beam model as 
indicated on the figure below. The attachments to the rocket frame may be 
considered to be nearly ideal clamps. The central load of block and gyros is 
symmetrically located; the bracket itself is replaced by an '’equivalent' 1 uniform 
beam. 


Rocket 

Frame 


A 


'A <2> Z ' C ' 


Xb 


4 


Mounting Block & Gyros 

, — IP / - 

j P iTCH f /4W 


Bracket 


o 

aotL 





R ocket F rame 


L 


Beam Model for Calculation of Bracket 
R esonances 
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The methods explained in Reference 4 produce the partial differential 
equation for free bending vibrations: 


E I b 



with 


+ M*L 


<Djrb 

d t 2 




la) 2 v b 
2 b y 


0 


E = Young's modulus 

Ik = area moment of inertia (for bending) 
k = mass per unit length in xj-, direction 
w zb = i nert ial angular velocity along axis 
= bending deflection 


The different modes are obtained by separating the variables x^ and t. The 
time dependent part of y^ becomes T (t) with differential equation: 

T (t) + (w 2 i - w 2 zb (t)) T (t) = 0 

while the numbers w^ are the characteristic values for the space-dependent 
part ^ (xb) with equation: 


EI b 



- w2 i ^ L ^ i = 0 


The boundary conditions for symmetric rro des are: 


<e j (o) = o 

QS 'i (°) = 0 

<£ 'i ( £ ) = o 

(!)+ = 6 

The last two conditions are to be replaced by other ones for antisymmetric 
modes, viz: 

«!(£)+ <£ ’i (£ ) - 2 ~ -S = o 

and a moment - equation 
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M c is the mass of block and gyros, augmented by the mass of the central part 
of the bracket. 


The lowest mode, i = 1, is symmetric and corresponds with linear motions 
along the roll gyro input axis and along the pitch and yaw gyros output axes. 


Integration of the spatial equation and its boundary conditions leads to the 
following relation for the lowest w^: 

V /2 

W1 " — > 


where the dimensionless number ^ is the smallest positive root of 


(sin cos h ol + cos c*. sin h oc ) + 


M r 




oijcos^cos h oc - 1 ) = 0 


The resonant frequency is wj as shown by the time dependent equation for T(t) 
provided 


I W z b (t) I << Wj 

The satisfaction of this condition will be verified later. Generally, w^ is high 
and, therefore, higher modes (which would affect angular rates) may be 
neglected in this study, since they are suppressed by the slower gyro-response. 

As a result, in conjunction with the various gimbal equations (2-6), one must 
use for yaw and pitch gyros 


a x = A -’ a 


A 

x> “z 


for roll gyro 


A x = 


1+2 > — + C—Y 

* br w br ^ wbr 


, A z a 


(2-14) 


with w]-j r = wj . ^br a dimensionless damping coefficient, the Laplace 
operator s is d/dt, and a is the inertial linear acceleration of the supporting 
rocket frame. 
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Numerical Results for the Gyro Mounting Bracket 


For the stainless steel bracket, the following values are used in the equivalent 
beam model. 

g = 0. 282 lb/in^ 

E = 2. 9 x 10 7 lb/in 2 
lb = 0. 020 in 4 
g M c = 1. 70 lb 
gpt-jj= 0. 102 lb/ in 
L = 8. 5 in 

£ = 2. 875 in 

Then 


, 1/2 4 2 

(E I b /p- L ) = 4. 7 x 10* in Z /sec 

smallest root oL ^ = 1.35 

and oO ^ = 10, 400 rad/ sec 

Since one can expect | cO ^(t) | ^ 7. 0 rad/ sec (i. e. , ten times the pitch gyro 

nominal input rate), the condition I cO , \ « oj is easily satisfied and one has 

. , , zb 1 1 

indeed 

^ br = = 10, 400 rad/ sec 

Since oO ^ is already large, higher modes can be neglected in this study. 

The value for gJ agrees with a very rough range prediction that can be 
obtained from the theory of plates. 

The damping ^ b r is unknown, but may be assumed to be smaller than 1 (no 
large structural damping). 
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2. 2 Laboratory Tests and Analog Verifications 


Several types of tests were performed on the gyroscopes. See Appendix B 
for a general description. 

2.2.1 Opened Instruments 

A pitch gyro and a roll gyro were opened, the damping fluid was removed and 
the gimbals were rotated mechanically and then released. The output waveforms 
are shown on Figures 2-1, 2-2, clearly indicating the very slightly damped 
motion of the gimbal. 


The experimental results are: 

Roll gyro ^ gy = ^80 rad/sec, 
Pitch gyro oO gy = 158 rad/ sec, 


f gy = 44. 5 Hz 
f gy = 25.2 Hz 


2.2.2 Transient Tests on Special Yoke, and Analog Simulation 

Figure 2-3 shows one test on the roll gyro. Knowing the total angle of rotation 
(1°) and the shape of the (upper) accelerometer trace, one can deduce the input 
waveform to the gyro. Figure 2-4 shows this input and the analog computer 
simulation response, as well as intermediate signals. The agreement between 
test (Figure 2-3) and simulation (Figure 2-4) is good (see e. g. , the heavy marks 
on the simulation showing experimental values). 


The parameter values used were, for the best fit, 

tO = 264 rad/ sec, f„ = 42. 0 Hz (roll gyro) 

gy &y 

\ gy =0. 84 

40 gX = 40 ^ gy 

Microsyn parameter values were listed above. 


Figure 2-5 and Figure 2-6 show, respectively, the test results and the simulation 
for another, more complex, input consisting of four succeeding pulses of 
alternating sign. The same parameter values were used and the agreement is 
acceptable between test and simulation. 
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FIGURE 2-1 OPENED ROLL GIRO 


TIME SCALE 10 ms/DIVISION 



FIGURE 2-2 OPENED PITCH GYRO 


TIME SCALE 10 ms/DIVISION 





FIGURE 2-3 TRIANGULAR INPUT TO ROLL GYRO 
TIME SCALE 5 ms/DIVISION 
UPPER TRACE: ACCELEROMETER 

LOWER TRACE: ROLL GYRO OUTPUT, 

7. 5 VOLTS/ VERTICAL DIVISION 
MAXIMUM ANGULAR ROTATION IS 1° 



FIGURE 2-5 MULTI-PULSE INPUT TO ROLL GYRO 

MAXIMUM ANGULAR ROTATION IS 2. 5° 

UPPER TRACE: ACCELEROMETER 10 ms/ 

DIVISION 

LOWER TRACE: ROLL GYRO OUTPUT 50 ms/ 

DIVISION AND 15 VOLTS/ VERTICAL DIVISION 
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This value for cO is consistent with the one found for undamped motion 
(section 2. 2. 1), within the accuracy obtainable which is of the order of 5 percent 
(see also section 7, Conclusions). 


2.2.3 Gyro Tests on Flight Table and Analog Simulation 


The inner gimbal of a flight table was used to generate particular inputs to the 
gyros. These inputs consisted of pulses alone or superimposed on sinusoidal 
inputs. Some pictures of the tests are shown in Figures 2-7 and 2-8 for the 
roll gyros and 2-9 and 2-10 for the pitch gyros. Upper traces are gyro responses 
to the inputs shown in the lower traces. The corresponding analog computer 
simulations are shown in Figures 2-11 and 2-12 (roll gyros) and 2-13 and 2-14 
(pitch gyros), respectively. Note that the output modulation was omitted, for 
clarity; therefore, the simulations give the envelopes of the gyro output waveforms 
(or in other words, the demodulated output of the gyros). The parameter values 
for the roll gyro are the same as the ones used in the previous simulations 
(section 2.2.2); those for the pitch gyro are 


od gy = 164 rad/ sec (fgy = 26 Hz) 

*gy = 0 ' 75 


The agreement between tests and simulation is well within the tolerances and 
uncertainties. 


2. 2. 4 Summary of Numerical Results for Gyros 

Pitch Gyros (40°/ sec) 


Output axis 


= 164 rad/ sec (26 Hz) 
& y 


Mechanical gain 
Stop 

Full Scale Output 


? gy * 75 


h t 


K 


= 0. 0647 sec 


gy 


0 max =3° 
9. 75 V rms 


0. 052 rad 


Roll Gyros (100°/ sec) 
t-0 = 264 rad / sec (42 Hz) 

l = 0. 84 

gy 

H r 

— : = 0. 0248 sec 

K gy 

9. 3 V rms 


Input axis 


Output Circuit 
Parasitic Effects 

Bracket 


w CT = 40 w 

gx gy 


= 0. 7 

is sufficiently large 


i gx 

(note: w 

See end of Section 2. 1. 2 
| g M gw c $ /H r | i 5. 3 x 10“ 4 


K c = oo 


or Q c = 0 


W br = 10, 400 rad/ sec (f br = 1, 


to be negligible) 
(rad/ sec) per g 
660 Hz) 
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FIGURE 2-7 ROLL GYRO CARCO TABLE TEST: 
PULSE INPUT 

TIME SCALE: 20 ms/DIVISION 

LOWER TRACE: INPUT (RATE) 

UPPER TRACE: GYRO RESPONSE 

10 VOLT /DIVISION 



FIGURE 2-8 ROLL GYRO CARCO TABLE TEST: 

PULSE PLUS 10 Hz SINUSOIDAL INPUT 
TIME SCALE: 20 ms/DIVISION 

LOWER TRACE: INPUT (RATE) 

UPPER TRACE: GYRO RESPONSE 

10 VOLT/DIVISION 
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Figure 2-9 


PITCH GYRO CAR CO TABLE TEST: 
PULSE INPUT 

TIME SCALE: 20 ms/DIVISION 

LOWER TRACE: INPUT (RATE) 
UPPER TRACE: GYRO RESPONSE 
10 VOLT /DIVISION 



FIGURE 2-10 PITCH GYRO CAR CO TABLE TEST: 

PULSE PLUS 10 Hz SINUSOIDAL INPUT 
TIME SCALE: 50 ms/DIVISION 

LOWER TRACE: INPUT (RATE) 

UPPER TRACE: GYRO RESPONSE 
10 VOLT /DIVISION 
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Figure 2-11, Pulse-Response of Roll Rate Gyro 
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3. ACCELEROMETERS 


Two types of single axis accelerometers used on the Scout vehicles are to be 
modeled here. 


Gulton Industries - LA 530250 (+3g) 

Gulton Industries - LA 460250 (-lg to +15g) 

Both types are servo-accelerometers. 

The procurement specifications require: 

an output voltage from 0 to -f-5 VDC into 500 K,Q. 
a natural frequency of 135 Hz (nominal) 
a nominal damping coefficient of 0. 6 to 0. 8 at 75°F 
crosstalk: 0. 010 g/g 

overload stops: 125% of full range 

The manufacturer expects a natural frequency of 90 Hz (-3g, +3g) and 110 Hz 
(-lg, +15g), with a damping around = 0. 7. 

3. 1 Theoretical Model 

3. 1. 1 Electro-Mechanical Part of Accelerometers 


These servo-accelerometers consist of a pendulously supported proofmass, 
whose motion is detected by an eddy-current pickoff. The pickoff signal, after 
demodulation, amplification and compensation, actuates a torque generator. 
Generally, the purpose of the feedback compensation is to obtain a well behaved 
second order response. However, experiments conducted during a Scout vehicle 
failure investigation indicate that above the break-frequency of the basic second 
order system, the system frequency response deviates markedly from ideal 
order behavior: the attenuation is higher (indicating an apparent increase in 

damping) while the phase shift corresponds to an apparent decrease in damping. 
While detailed information on the internal structure (especially damping) and 
the compensation circuits could not be obtained, an incomplete block diagram, 
shown below, can still be constructed. 


L s* 




Detector 

_ _ K P _ 

__ Angle 


Amplifier 


K 


A 




Torquer 


Compensation and Damping Omitted 
Block Diagram for Servo-Accelerometers 


A *~'% 


v.. 
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The notations are: 


Input acceleration (g's), (inertial acceleration of case along 
sensitive axis minus gravitational acceleration on proofmass along 
same axis) 

A QU t : Output reading converted in g's 
V Qut : Direct output reading in volts 

2 

I: Moment of inertia of proofmass with respect to suspension (g cm ) 

L: Radius of gyration of proofmass referred to support (cm) 

K : Gain of detector or signal generator (V/rad) 

P 

Amplifier gain (A/V) 

Krj,: Torquer-constant (dyne cm/A) 

S^: Scaling factor to g’s ( crn /a) 

sec^ 

Sy: Scaling factor to volts ( V/A) 

s: Laplace transform operator 

The linearized transfer function is then 



A in L K T 1+ s 2 I 

K p K A K T 

L K t 

indicating that S^ = _ 

and that the natural (undamped) frequency is given by 

w2 n = K p K A K T /;[ t 3 - 1 ) 
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For a more compiete model, damping and cr os s- coupling must be considered. 
Due to misalignments and structural constraints, cross -coupling can be 
expected, i. e. , case accelerations A crogs perpendicular to the sensitive axis 
will influence the readings. Lack of knowledge of instrument structure prevents 
theoretical modeling of the cross-coupling effect. Therefore, the same dynamic 
behavior is adopted for cross -acceleration, with a scale factor 4 for which 
an upper bound can be found from manufacturer's data. These considerations 
lead to the following model: 


out 




1 + 2 > i — 

1 w n 

where is bounded 


1 + 2 f -J— + 


v n 


w n 


vp I ^ ^ 


max 


(3-2) 


(3-3) 


The "angle 11 ^ max * s P art ^y determined by the maximum angular rotation of 
the proofmass. 


Some Numerical Results for the Servo-Accelerometer s (see also Appendix B) 


A (-lg, + 15g) instrument was opened and open loop information was obtained. 
Measurements of the proofmass and its connecting system led to a calculated 
value of the inertia 

I = 0. 213 g cm 2 (+ 10%) 

and a weight of 0 . 9 gram. 

Static experiments on the torquer (using balancing weights added to the proofmass) 
gave 


Kp = 0. 139 x 10^ dyne cm/ A (+ 5%) 

The signal generator (Kp) and amplifier (K^) were combined in the open loop 
tests and yielded 

K p = 1.22 A/rad (+ 10%) 

Hence, one finds for the (-lg, +15g) accelerometer 

w n = 893 rad/ sec (+ 15%) 

Damping could not be determined after the instrument was opened. 

Manufacturer's data show vi/ = 0. 002 

1 max 
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3.1.2 Accelerometer Bracket 


The accelerometers are mounted on stainless steel brackets which are then 
fastened to the rocket frame, The configuration is shown in Figure 3-1. The 
bracket is very stiff along the sensitive axis and along the perpendicular cross - 
axis, but less stiff along the lateral cross-axis. The purpose of this section 
is to calculate the resonant frequency for transmission along the lateral cross- 
axis. 

First one calculates the stiffness k c of the corners, with definition: 

k _ moment exerted on a corner 

c total additional bending angle of corner (radians) 

Using the theory of bending of curved bars (Reference 4), one finds: 


t/ r c 

1 + (t/2 r c ) 
e 1 - (t/2 r c ) 

where E is Young* s modulus and the other symbols are illustrated in the 
Figure 3-2 (w = width, t = thickness, r c = rounding radius). 
ln 0 is the natural logarithm. Next, one calculates the stiffness kL for the 
linear displacement deformation of Figure 3-2. 

kj^ = force along lateral cross -axis through instrument 

linear displacement (of instrument) 


k ‘ = 5r 


E w t r 


1 - 


kL depends on k c and the properties of the bending members. Application of 
beam theory gives after some manipulations 




+ 


SS-5-) 

Ew t 3 / 


The approximate resonant frequency of the bracket w^ is then given by 
w br = (2 k L /M eff ) 1/2 

where M g ££ = mass of accelerometer plus mass of top plate of bracket. 
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Deflection 
Driving Force 


Elastic Support (kc) 
Bending Member 

L Accelerometer Bracket 
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Then, for accelerations along the lateral cross -axis. 


^■lateral, 


cross 


a lateral, cross 


1+2 ^br 


1 

s 

w br 


+ 


s 2 

w ^br 


(3-4) 


where -A^ atera j cross a parasitic input to the accelerometer, and 
a lateral, cross the inertial acceleration experienced by the support. 

For the other axes, one has: 


sensitive axis A 

. A - a 

perpendicular cross-axis 

The coefficient ^ is the dimensionless damping coefficient for the bracket 
vibrations. 


Numerical Results for the Accelerometer Bracket 


The same bracket is used for both types of accelerometers. 

^ is unknown, but may be expected to lie between 0 and a maximum slightly 
larger than 1. 

For the stainless steel bracket with weight density = (g f ) = 0. 282 lb/ in 

E = 2. 9 x 10 7 lb/in 2 

w =3 in. 

t = 0. 047 in. 

a = 0. 625 in. 

length of top plate = 1. 375 in. 

r c = 0. 1 in. 

one gets 

k_ = 0. 49 x 10 4 in -_ lb - 
c rad 

k T = 1. 5 x 10 4 Ah. 

^ in. 

The weight of the top plates is 0. 055 lb. Then (g M e f£) = 0. 274 lb. 
finally giving w^ r = 6, 500 rad/ sec (f^ r = 1, 040 Hz) which falls within the 
frequency range of interest. 
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3. 2 Laboratory Tests and Analog Computer Simulations 


As for the gyros, several types of tests were performed on the accelerometers. 
See general description in Appendix B. 


3. 2. 1 Transient Tests on Special Yoke, and Analog Simulations 

Both categories of accelerometers (-3g, + 3g) and (-lg, + 15g) were subjected 
to triangular input pulses, shown in Figures 3-2, 3-3, 3-4 and 3-5. The 
corresponding analog computer simulations are shown for comparison in 
Figures 3-6, 3-7, 3-8 and 3-9. The agreement is within acceptable tolerances. 


The best fit was found for the following parameter values: 


(-3g, +3g)* 


(-lg, + 1 5g) 


Resonance w n (rad/ sec) 

Damping 1 n 

Gain 


800 (fast input) to 850 (slow input) 

0. 35 (fast input) to 0. 55 (slow input) 
0. 85 V per g 


750 
0. 35 

0. 34 V per g 


Nominal output range 0 to 5. IV 


*Note : For another ( - 3g, +3g) instrument, w n = 700 rad/sec, Y n = 3 (fast) 

to 0. 45 (slow) was found to give a best fit. This remark gives an idea of the 
possible spread among instruments, of the order of 10% to 2 0%. Such a typical 
spread immediately yields a bound on the accuracy for which one should strive 
in obtaining parameter values for a series of instruments in the present study. 

Of course, for one particular instrument one may try to achieve a higher accuracy 
in identifying parameter values. One also observes a slight non-linear effect 
in the decrease of the damping factor 'f n , when the input gets faster, at least 
for the (-3g, + 3g) type. 


3. 2. 2 Carco Table Tests for Accelerometers and Analog Verification 

Pulses, alone or super -imposed on 10 Hz sine values, were used as inputs to 
the accelerometers. Some experimental results and analog simulations are shown 
in Figures 3-10, 3-11 3rl2,3-13, respectively. The parameter values used are 

the same as before (in section 3.2. 1). There is some discrepancy, the negative 
value of the first undershoot in the simulations, but this is probably due to 
uncertainty caused by ringing in the very first part of the input waveform. 
Elsewhere, the agreement is acceptable. 
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FIGURE 3-2 SLOW TRIANGULAR INPUT TO 
(-3g, +3g) ACCELEROMETER 
TIME SCALE: 2 ms/DIVlSION 

UPPER TRACE: INPUT (PIEZO-ELECTRIC 

ACCELEROMETER ) 

LOWER TRACE: INSTRUMENT RESPONSE 

0. 5 VOLT / DIVISION 



FIGURE 3-3 FAST TRIANGULAR INPUT TO 
(-3g, +3g) ACCELEROMETER 
TIME SCALE: 2 ms/DIVISION 

UPPER TRACE: INPUT (PIEZO-ELECTRIC 

ACCELEROMETER) 

LOWER TRACE: INSTRUMENT RESPONSE 

0. 5 VOLT /DIVISION 
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FIGURE 3-4 SLOW TRIANGULAR INPUT TO 
(-lg, +15g) ACCELEROMETER 
TIME SCALE: 2 ms/DIVISION 

UPPER TRACE: INPUT 

LOWER TRACE: INSTR UMENT RESPONSE 

1 VOLT /DIVISION 



FIGURE 3-5 FAST TRIANGULAR INPUT TO 
(-lg, +15g) ACCELEROMETER 
TIME SCALE: 2 ms/DIVISION 

UPPER TRACE: INPUT 

LOWER TRACE: INSTRUMENT RESPONSE 

0. 5 VOLT /DIVISION 
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Figure 3-7, (-3g, + 3g) 
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Figure 3-8, (-lg, +15g) Accelerometer No. 119 
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FIGURE 3-10 CARCO TABLE TEST FOR 

(-3g, +3g) ACCELEROMETER 

TIME SCALE: 20 ms/DIVISION 

LOWER TRACE: INTEGRAL OF PULSE INPUT 

(A VELOCITY) 

UPPER TRACE: INSTRUMENT RESPONSE 

1 VOLT/DIVISION 



FIGURE 3-11 CARCO TABLE TEST FOR 

{ -3g, +3g) ACCELEROMETER 
TIME SCALE: 20 ms/DIVISION 

LOWER TRACE: INTEGRAL OF INPUT 

(PULSE PLUS 10 Hz SINE) 
UPPER TRACE: INSTRUMENT RESPONSE 

2 VOLT /DIVISION 
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4. PRESSURE TRANSDUCERS 


This section deals with two types of pressure transducers used aboard Scout 
vehicles. 

Giannini - 461319BV-2 (0 to 500 psia) 

Bourns - 2007253703 (0 to 800 psia) 

Both types use a Bourdon - tube as pressure sensing device. The Giannini 
transducer is connected directly to the pressure source. The Bourns transducer 
is connected to the pressure source by means of an asymmetric T-tube, whose 
third arm leads to a switch chamber. 

An extract from the procurement specification follows: 

- Bourdon tube material Ni-Span-C 

- Acceleration sensitivity: + 0. 05%/g along lateral and 

transverse axes, _+ 0. 004%/g along the longitudinal axis. 

- Response time: time for the output to reach 63% of its final 

value after a step input in pressure shall not exceed 

50 milliseconds. 

4. 1 Theoretical Model 

4. 1 . 1 Bourdon - Tube Analysis 

Static elastic analysis of Bourdon tubes is quite complicated, and dynamic 
analysis is even more involved. Therefore, only approximate results can be 
obtained. Both transducers use flattened sections, shown in Figure 4-1, (see 
also Reference 5). Giannini employs a helical tube, while Bourns has a spiral 
tube. Therefore, the treatment of both cases is somewhat different and will 
be given separately. Basically, the resonant frequencies are obtained here by 
an application of Rayleigh* s principle (see Reference 4). One equates approximate 
expressions of maximum potential and kinetic energies for the lowest mode of 
vibration. However, since the trial functions used do not necessarily satisfy 
all boundary conditions, the resonances found may be smaller as well as larger 
than the true values. 

Giannini (or Helical) Tube 


The Giannini tube is a regular helix wound on a circular -cylindrical surface. 

The radius of curvature R of the axis of the tube is a constant; therefore, the 
theory of Reference 5 can be applied directly if the helix is replaced by equivalent 
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Cross-Section of Bourdon Tubes 
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Figure 4-2, Tube Nomenclature 
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circular rings. Application of pressure results in an increase in radius of 
curvature cFR. Reference 5 gives the following relationship between torque 
on tube, T, and the angular deflection of the tube, per unit of the angle 9, 

(see Figure 4-2) cT R/R (assumed constant over the entire axis of the tube, an 
assumption equivalent with the hypothesis of constant bending moment). 


T = k • 

R 



E = Young* s modulus 

v = Poisson*s ratio 

a, b, h, R: See Figure 4-1, geometric quantities describing the tube 

cross-section. 

Fj_, F^: dimensionless coefficients, tabulated in Reference 5. 

The corresponding maximum potential energy of deformation is then given by: 



In a helix with constant length of the tube axis, the displacement u (see Figure 
4-2) is given by 

u^ = ( cTR) 2 \_2 (1 - cos 9) + 9^ - 2 9 sin 9 J 

In harmonic motion with angular frequency wg, the corresponding maximum 
kinetic energy is (for m£ = tip mass) 

^kin = \ W B 2 kin 

2 3 

Z kin = f A r 3 C'^7(-| £+ 2 S f - 4 sin 6- f + 2 9 f cos 9 f ) + m f u f 2 

mass density of tube (metal plus pressurized fluid) 
cross-sectional area given by 

A = 4h a [l +(~- -l) ±] 


where 


f = 

A. = 
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Therefore, the approximate resonant frequency of the undamped tube is 


^ Upot/ Z kin 


The mass of the pressurized fluid should be added as a correction to 
Generally, this is a very small effect. 


When the tube is immersed in a damping fluid, some of this damping fluid will 
oscillate with the metal tube. This phenomenon can be interpreted as an 
additional apparent mass which has to be added to ^ as the correction: 

(* eff f + ^ apparent 


This value of ^ e ff should be inserted in the formulas above for greater accuracy, 
yielding a value e ff. Again, this is normally a small correction. 


Then, the resonant frequency with damping fluid, wgj), is given by: 


W BD 


2 


2 Upot 
Z kin, eff 


The transfer -function for the Bourdon tube with damping fluid now takes the 
form, if the damping is included, and with proper scaling (i. e, , zero steady 
state error). 

Pout = 1 

£ +(JS 

BD W BD W BD 


1 + 2 


(4-1) 


However, the tube is also sensitive to accelerations. The lowest mode for 
accelerations will be close to the lowest mode for pres sure -input. Therefore, 
the following transfer -characteristic can be written down: 


^out 


1+2 f 


BD 


BD 


+(— — ) 

> w ' 


BD 


i=l 


*i PFS 


1+2 


A. 

i 


BD 


+/S 


-) 

BD 


(4-2) 


54 



where 


A. = acceleration (expressed in units of "g") of the transducer 
along one of three perpendicular axes (i = 1, 2, 3) 

= full scale* or nominal, pressure 

^ = sensitivity coefficients (dimensionless numbers) indicating 

fractions of full scale per g along the three perpendicular 
axes (i = l f 2, 3). Because of the directional properties of 
the Bourdon tube, the ^ -coefficients are different for each 
direction; they also depend on the properties of any balancing 
masses that may be used. Their bounds can be deduced from 
manufacturer^ data and from the specifications. 

Numerical Results for the Giannini Transducer 


The tube material is always Ni-Span-C 

7 2 

with E = 2.6 x 10 lb/ in 

(gf) = 0.293 lb/in 3 
2 

g = 386 in/ sec 
v = 0. 3 

From manufacturer's data, confirmed by opening of an instrument, one finds: 

a = 0. 080 in. 

b = 0. 021 in. 

h = 0. 008 in. 

R = 0. 225 in. 

0f = 1 1 TT (i. e. f 5 1/2 turns) 

P-pS = ^00 lb/in^ (absolute) 

Reference 5 gives Fj = 1. 058, F 2 = 1. 451 

One calculates immediately 

k = 40. 3 in. lb. 

A = 2. 94 x 10 -3 in 2 
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At full pressure and a temperature of 535°B (or 75°F), using nitrogen as the 
medium, one has 

mass tube = 125 

mass compressed gas 

and hence the mass correction due to the compressed gas is entirely negligible. 
The damping fluid correction is also negligible. 

Finally, one has 

/ x 1,990 

w bd ( rad / sec ) - ~r TT ,/i - 

m^£ being the tip mass, expressed in grams. 

Opening an instrument showed that the total tip mass is not more than 1. 4 gram. 
The effective tip mass, taking into account mass distribution and support is 
smaller. Hence, one has: 1,260 < Wg£> <. 1,990 (rad/sec). 

A probable value for effective tip mass is about 1 gram, yielding 

W BD = 380 rad/sec. 

The manufacturer's estimate was wg£) = 1,250 rad/sec. Manufacturer -supplied 
data indicate a bound for the acceleration sensitivity coefficients. 

^ 2 x 10 ^ for all axes in fractions of full scale per g. 

These values for I ^ i | agree with (i. e. , are smaller than) a theoretical 
upper limit found by assuming no balancing is available. 

Opening of Giannini Transducer (See ALso Appendix B) 

One transducer was opened and the damping fluid removed. An impulsive 
mechanical displacement was then applied to the Bourdon tube and the response 
is exhibited in Figure 4-3. One deduces from it that wgp) = 1, 250 rad/sec, 
which is compatible with the results above. Shock inputs to the damped tube 
also roughly indicated a damping factor ^ of the order of 0. 2, i. e. , a rather 

lightly damped system. 
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FIGURE 4-3 RESPONSE OF UNDAMPED GIANNINI BOURDON TUBE 
(OPENED INSTRUMENT) 

TIME SCALE: 5 ms/DIVISION 



FIGURE 4-4 RESPONSE OF UNDAMPED BOURNS BOURDON TUBE 
(OPENED INSTRUMENT) 

TIME SCALE: 2 ms/DIVISION 
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If 





Bourns (or Spiral) Tube 


For the spiral tube R varies with 9 and the theory of Reference 5 is not directly 
applicable. The method used here is the well-known quasi- steady approximation, 
i. e. , R is taken constant over short lengths of the tubing axis. The manu- 
facturer's drawings show that the Bourdon tube is made up of consecutive half 
circles with increasing radius. This real shape is replaced by an equivalent 
spiral of Archimedes with polar equation: 

r = q s 6 


taken between two limits: 

0 = 9i n (fixed end, where fluid enters) 
0 = 0£ (free, moving end of tube) 

The coefficient q g is the appropriate average value: q s 

giving 6 in = r in /q s . 


The lower limit: 
and 

The upper limit: 


r f - r in 
0 f - 6 in 


For 9 not too small (i. e. , 0i n >> 1) polar radius r and radius of curvature R 

are nearly equal for the spiral of Archimedes as they must be for the successive 
half circles. 


R =, r 


q s 0 


The bending theory of curved beams (Reference 4) shows that for a constant bending 
moment, one has <5* R/R. (the angular deflection of the gage per unit of its 
length) nearly constant. Therefore, one can write 

T = (kR) av 


where 


(kR )_ = 4 E h b 2 a 

' 'ave 


+ (i - *) ; 



av 
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The maximum potential energy of deformation now is given by 


V = - (kR) 
pot 2 av 


£r‘ 


(R in /R f ) 


2 q s 


Again for © not too small, the displacement u is given approximately by 

l u l 5 ( 0 - «in) 

since, approximately, 

+ 1 U I = R S 0 = - (G - e in ) cT R 

Therefore, the maximum kinetic energy for oscillations with angular frequency 
wg becomes 

V kin * { w2 B Z kin 

with 

z k in = A f / u 2 r d e 

tube + tip mass 



The expression for the tube area A is given above in the discussion for the 
helical tube. The remainder of the discussion for the spiral tube is exactly 
the same as for the helical tube and leads to identical forms for the final 
transfer functions. 
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Numerical Results for the Bourns Transducer 

The tube is made of Ni-Span-C (see Giannini data). According to data from the 
manufacturer, confirmed by opening an instrument (except for tip mass) 

a = 0. 148 in. 
b = 0. 046 in. 

R in =0. 28 in. 

Rf =0.61 in. 

h = 0. 01 4 in. 

0 = 0f - 9 in = 6. 14 IT (or 3 complete turns plus an additional 25°) 
PpS = 800 lb/in^ (absolute) 

Total tip mass between 2 and 4 grams (this is doubtful 
as indicated below) 

Hence, one finds 

q s = 0.0172 in/rad 
0 in = 16. 4 = 5. 2 TC 
0f = 35. 8 = 11.4 rc 

From Reference 5, one obtains Fj = 1. 053, F 2 = 1. 506 
Also: A = 0. 976 x 10" 2 in 2 

The mass correction due to the presence of pressurized gas is again negligible, 
since for nitrogen at 800 lb/ in^ and 535°R 

mass tube _ ^ 

mass compressed gas 

The damping fluid correction is_ also negligible. One finds (k R )^ n = 109 in^ lb 
and (k R)f = 166 in^ lb; an average value (k R ) av = 138 in^ lb is taken. 


1 , 170 

(1" + 0.358 m* f ) iyi 

with m*f expressed in grams 


Finally one finds 

W BD ( ra d/sec) ~ 
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For m*£ = 0. 6, an approximate value found by opening an instrument, one 
has wgQ = 1, 060 rad/sec. 

Also, according to the manufacturer, one finds bounds for the acceleration 

sensitivities [ i \ ( p | 4 5 x 10"^ (most sensitive axes) 

1 j c 

I 4 4x10"^ (least sensitive axis, perpendicular to 
the plane of the spiral) 

expressed as fractions of full scale per g. 

Opening of Bourns Transducer (See Also Appendix B) 

The procedure here is exactly the same as was already explained for the 
Giannini transducer. Figure 4-4 shows the results, from which one deduces 
Wg£j = 1, 050 rad/ sec, in very good agreement with the numbers arrived at above. 
A shock test of the damped (unopened) tube roughly pointed to a f of the 

order of 0. 2. 


4. 1 . 2 Electrical Output Devices 

The displacement of the free end of the Bourdon tube is converted into a (direct) 
output voltage across a resistor by means of a moving wiper. Because of lack of 
detailed information on the output device, its dynamics are neglected. However, 
this does not affect the instrument model much, since the output device responds 
much faster than the mechanical part. 

4.1.3 Connecting Pipes (Bourns -Transducer ) 

The Bourns transducer is connected with the pressure source by means of an 
asymmetric T, shown in Figure 4-5. 

The model for the pipes is the transmission line with longitudinal losses, shown 
in Figure 4-6. The transfer relations for it are derived in Reference 6, and 
the end results are listed below. 


cos h ( r s) 


Z Q sin h ( r s) 


-L. sin h ( r s) cos h (Ts) 
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Switch Chamber Junc+Lo^ 



Length of Tube 1 = A j, Volume = V^ 

Length of Tube 2 = A 2 > Volume = V£ 

Length of Tube 3 = A 3 , Volume = V 3 


Figure 4-5, Connecting T-Pipe 


"Generator -End ' 1 

Pq ° 

*• 

O 

Length = A 


"Load" -End 



Pressure 
Volumetric Flow 


<?oo -» — 'irrr — 'vwv 


i d X 

O 


r c c/x 




-r c 0 c/a 


4 

p (a i-d xj 


X <- 


d x 


-V X ‘f- 0/ x 


Figure 4-6, Lossy Transmission Line Model 



2 Q 

with p = pressure (say lb/in ), Q = volumetric flow (say in /sec), Z Q is the 
characteristic impedance 


Z 

o 




1/2 

) 


and 7T is the "time constant" 

r 1/2 

^ = A ( c o< 1 o + -^ » 

. 4 

& is the length of the pipe, c e is the ’’capacitance" per unit length (say in _ , 

lb 

1 is the "inertance" per unit length (say lb sec^/in^), r Q is the "resistance” 
per unit length (say lb sec/in^), and s is the operator d__ . For a compressible, 

dt 

ideal gas one has the pressure-density relation. 

P * f R T 

The flow is to be taken as isothermal in this case: 


T = constant = T Q 
Then one finds 

- 1 o 7 

c o 1 o = ( R T 0 ) ( sa V- sec / in ) 

c o r c = 2 f (D \ r R t”)' 1 (say, sec/in 2 ) 

where f is the usual dimensionless friction coefficient (defined as ratio of 
pressure drop per quarter inner diameter length to dynamic pressure) and D 
is the inner diameter of the pipe. 


The transducer can be considered as a constant volume device, i. e. , as a 
capacitance C^ n with 


C in 

C 2 


m 

Co^2 
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The switch has two positions, corresponding with two chamber volumes. At 
the pressure of 50 lb/in 2 the chamber goes from a smaller volume (0. 07 in 3 ) 
to a larger volume (0. 14 in 3 ) and stays there. Since 50 lb/ in^ is rather small 
compared with the nominal pressure (800 lb/in 2 ), one can take V sw = constant 
(and equal to the larger volume) for the present approximate analysis. 

Then the switch acts as a capacitance C sw with 


C o A 3 


(Capacitance being proportional to volumes) 


One finds the impedance ratios: 


v r c - Z o C sw ( 
^o u sw s “ 

*3 


Zo C in s 


Repeated application of the basic transmission line formula and of the 
relations 


C sw s 


leads to the transfer-relation: 


^ = (cos h(-r 2 s) + Z Q C in s sin hfT^s)) jcos hfcis) 

p g 

'sin h(r* 3 s) + Z Q cos hfir^s) sin h(r" 2 s ) + ^o^in s cOS 

+ sis Mt^s) h(T 3 s) + Z Q C sv< ?sin h(-c 3 s) cos h(r 2 s) + Z Q Ci n s sin h(x- 2 s 
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This complicated transcendental transfer function is now replaced by the 
low order rational approximation. 


P i» 


1 + 2 h 


(4-4) 


1 + 2 


b 


w 


Best values of w a and w-^ are found by a product-development of the transcendental 
numerator and denominator of pj n /p s * 

The zeros of Pi n /p s (or the poles of P s /pi n ) are given by the equation 


cos h (X 3 S) + Z Q ^sw s sin h (r^s) 


»/z 


If one sets T^s = C-l) v, v is real for the roots of the preceding equation. In 
fact, the zeros v^ are given by 


(tan v n ) 


= V3/ v 


sw 


zeros 


The smallest positive root ( v l) zero yields 


( V 1 ) 


w a 


A 3 ( A oC 0 ) 


1/2 


The damping is given by: 


2 i ^ w a 


The poles of P^ n /p s (zeros of P s /pi n )» including w^, are given by the (real) 
roots of the equation in v: 


/ ^ 1 \f / ^ 2 \ V in . ( ^ 2 \ 

cos ( v )L cos (tt v )- vsm (xr 7 


sw 


v sm 


in vj 


. [A 1 


“l , r .^ 2/3 \ t , v ,w v m Z\ , , 4 2+ a 3 l v int v sw 1 

it v iL s “v-i— V K ‘ — v ) + c °<~ a — y > V J 
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One has for the pipes 

3 

the volumes v\ = 0. 17 in 

= 0. 26 in^ 

V 3 - 1 o 1 3 in^ 

the lengths t\ ^ = 5 . 38 in 

A 2 = 80 27 in 

^ ^ = 36. 0 in 

and inner diameter D = 0. 20 in» 

3 

The switch has a volume v sw = 0. 14 in J 

3 

The Bourdon tube has a volume v. = 0. 17 in 

in 

One takes (for the laboratory experiments) 

T q = 535°R (75°F) 
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Obviously, for real Scout-data, the appropriate (possibly) different temperature 
must be taken. 


The gas constants are 

for nitrogen Rp^ = 2. 56 x 10 


for helium -^He = x 10 


in fc 


sec 

. 2 
in 


sec 


R 


R 


Damping estimates are found to be unreliable; they give only an order of 
magnitude estimate. It is best to rely on experimental evidence to find f a, £ b. 
Since helium was used in the laboratory experiments, the constants were cal- 
culated for this gas. 

One finds successively 

c Q l Q = 1, 04 x 10"^ sec^/ in^ , — - - /2 = 3,1 x 10^ in/ sec 

0o c °v 


^ v l^zero 

= 1. 45 (corresponding with 82. 9°) 

( v l>pole 

= 1.25 (corresponding with 71.6°) 

w 

= 1, 250 rad/sec 

a 


w b 

= 1 , 080 rad/ sec 


One sees that w-^ is only a little smaller than w a . Theoretically, as discussed 
before, one should then take more zeros and poles in the tube-model. However, 
the need for additional terms in the tube -model is obviated by the effect of 
orifices, to be discussed in the description of the experiments later on (Section 
4.2.3). 
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4. 2 Laboratory Tests and Analog Computer Simulation (See also Appendix B 
for Description) 

4. 2. 1 Giannini Transducers 

Figures 4-7 A&B, 4-8 A&B, 4-9 A&B, show inputs and outputs for tests per- 
formed on one Giannini instrument. There is a connecting orifice between the 
transducer and the pressure source, and this connecting orifice accounts for 
a second quadratic term besides the basic Bourdon transfer -function derived 
in Section 4. 1. 1. The total transfer -function is now to be taken as 


out 

P in 


(4-6) 


(1+2 f 


w 

or or 


) (1+2 \ 


or 


BD 


/ BD 


BD 


Full scale gives 4. 94v out. 

Figures 4-10, 4-11 and 4-12 show, respectively, the best fits obtained on the 
analog computer for the preceding inputs. These fits are quite good for the 
range of inputs considered. 


Summarizing, one finds for the 


Bourdon tube 


W BD = 300 rad/sec (fg D = 207 Hz) 

(bd'°- 2 


agreeing well with the theoretical values derived in Section 4. 1. 1. 


For the connecting orifice, one gets v 

w = 1, 300 rad/sec to 1, 600 rad/sec 
or 

= 0. 3 (slow input) to 0. 7 (fast input) 

again showing some non-linear effect in the gas flow, with damping increasing 
as the input gets faster. The spread among instruments of the same category 
reaches at least 10 percent. 
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FIGURE 4-7 A INPUT TO GIANNINI, 195 PSI FINAL 
TIME SCALE: 1 ms /DIVISION 

(POSITIVE LEFTWARD) 



FIGURE 4-7 B GIANNINI RESPONSE 

TIME SCALE: 2 ms/DIVISION 

SMALL FIXTURE WITH 0. 150" ORIFICE 
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FIGURE 4-8 A INPUT TO GIANNINI, 150 PSI FINAL 
TIME SCALE: 1 ms/DIVISION 

(POSITIVE LEFTWARD) 



FIGURE 4-8 B GIANNINI RESPONSE 

TIME SCALE: 2 ms/DIVISION 

SMALL FIXTURE WITHOUT ORIFICE 
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FIGURE 4-9 A INPUT TO GIANNINI, 200 PSI MAXIMUM 
TIME SCALE: 5 ms /DIVISION 

(POSITIVE LEFTWARD) 



FIGURE 4-9 B GIANNINI RESPONSE, 235 PSI PEAK 
TIME SCALE: 5 ms /DIVISION 

SHOCK TUBE WITH FAST RISE 


A 

■ 
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Figure 4-10, Giannini 



lb/ in^ (No. 465-1) 
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Figure 4-12, Giannini 500 lb/in^ (No. 465-1) 







4.2.2 Bourns Transducers without T -Tubing 


In order to test the value of the theoretical calculations of Section 4. 1. 1, the 
Bourns transducer was subjected to tests without T-tubing (in order to avoid 
masking the Bourdon tube response by other components). 

One such test is shown in Figure 4-13 A&B, and the corresponding best analog 
computer fit is on Figure 4-14. As in the case of the Giannini tests of the 
preceding Section 4.2. 1, there is a connecting orifice. The best fit parameter 
values are 


for the Bourdon tube: W BD ~ ** 000 rad/sec (160 Hz) 

\ BD = °* 1 2 

and for the orifice effect: w^„ = 1. 000 rad/ sec 

or 

Full scale corresponds with 4. 94 v ou t. 

The Bourdon resonance value from this test agrees very well with the theoretical 
prediction of Section 4. 1. 1. 


4. 2. 3 Bourns Transducers with T -Tubing 


Figure 4-15 shows the response of the T-tubing and Bourns transducer for what 
was essentially a step input in pressure (using helium as a pressure medium) 
and Figure 4-16 exhibits the best fit obtained on the analog computer. The 
response at the T-tube end is also shown. The entire numerical transfer function 
(best fit) can be written as 


1 + 2 x 0. 8 


1 , 000 


(i f 

K 1 , 000 ' 


1 + 2 X 0. 8 - — + (2 — ) 

200 ' 2nn I 


200 • 


(4-7) 


1 + 2 x 0. 3 - — + 7 - — ) 
400 ^ 400 ' 
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FIGURE 4-13 A INPUT TO BOURNS TRANSDUCER 
TIME SCALE: 5 ms /DIVISION 

(POSITIVE LEFTWARDS) 

700 PSI PEAK 



FIGURE 4-13 B BOURNS RESPONSE (NO T-TUBING) 
TIME SCALE: 5 ms /DIVISION 

SHOCK TUBE WITH FAST RISE 
520 PSI PEAK 
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-15 STEPRESPONSE OF T-TUBE AND BOURKS-TRAMSDUCER 

TIME SCARE: 5 ms/DIVISION (POSITIVE LEFTWARDS) 




While this expression has the theoretical form of Sections 4. 1. 1 and 4. 1.3, 
Equations 4-1 and 4-4, some remarks are in order. The numerical numerator 
above agrees with the theoretical w a calculated in Section 4. 1. 3 (1, 250 rad/sec) 
within the expected accuracy. The denominator factors, however, are 
numerically widely different from w^ and wgj) calculated before. This discrepancy 
is probably due to relatively slow orifice effects between pressure source and 
T -tube on the one hand and between T -tube and Bourdon tube on the other hand, 
whose effect is to mask the faster response connected with w^ and wgp. The 
Bourdon resonance wgp was indeed experimentally verified as valid in the results 
of Section 4. 2. 2. The comparison of the results of 4. 2. 2 and 4 0 2. 3, moreover, 
points out that the numerical transfer functions obtained must be considered as 
global (overall) descriptions of individual systems and that it is dangerous to 
try and split these expressions in an effort to identify individual components. Such 
a conclusion is expected because of the interaction of successive components in 
a fluid flow system. A good example of this interaction is afforded by the 
discussion of Section 4. 1. 3, which clearly shows that Pi n /p s for the T-tube 
depends on the type of termination (here defined by the ''admittance" Ci n s). 
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5. TELEMETRY STUDY 


5. 1 General Description 

The airborne, NASA -supplied equipment which is discussed in this report 
consists of: 

a) An S-band Telemetry Transmitter CTM-UHF-305 (Conic Corp. , 

San Diego, California) 

b) A Lightweight Telemetry Package TDD 1799A (Tele-Dynamics, 
Philadelphia, Pennsylvania) 

The principal Avco ground -equipment used in conjunction with the NASA material 
is an EMR Model 229 Tunable Discriminator made by Electro-Mechanical 
Research, Inc., Sarasota, Florida. Some other equipment was also employed 
and its name is given where applicable. 

For the purposes of this study, the following block diagram of the T/M chain is 
applicable: 



Transducer s 


Airborne Part 



Ground Equipment 
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The various components will be discussed in the logical order they occupy in 
the block diagram. 

The various instruments transmit data over a number of standard channels of 
the IRIG Telemetry Standards for FM/FM (see Reference 7). The relevant 
channels are Listed in the table 5-1 below. 

TABLE 5-1 

IRIG CHANNELS 


Maximum deviation is ■+ 7. 5% of center frequency of each channel. 


Band 

Number 

Center 

Frequency 

(kHz) 

Fc 

Lower 

Limit 

(kHz) 

Upper 

Limit 

(kHz) 

Frequency 
R esponse 
(Hz) 

Instrument and 
Appendages with Lowest 
Proper Frequency in Hz 

5* 

1. 3 

1. 202 

1.399 

20 


6 

1. 7 

1. 572 

1. 828 

25 

Yaw-rate gyro (26) 

7 

2. 3 

2. 127 

2. 473 

35 

Roll-rate gyro (42) 

8 

3. 0 

2. 775 

3. 225 

45 

Pitch rate gyro (26) 

10 

5. 4 

4. 995 

5. 805 

81 

Bourns pressure 

11 

7. 35 

6. 799 

7. 901 

110 

transducer (32) (2nd stage 
headcap pressure) 

Bourns pressure 

12 

10. 5 

9. 712 

t— * 
i — ■ 

ro 

00 

00 

160 

transducer (32) (1st stage 
headcap pressure) 

Also: Giannini pressure 

transducer (3rd stage 
headcap pressure) (207) 

Transverse accelerometer 

13 

14. 5 

13. 412 

15. 588 

220 

(-3g, + 3g) (110 to 135) 
Normal accelerometer 

14 

22. 0 

20. 350 

23. 650 

j 

330 

(-3g, +3g) (110 to 135) 
Longitudinal accelerometer 

15 

30. 0 

i 

27. 750 

32. 250 

, 450 

(-lg. +15g) (119) 


*NOTE: Channels 5 and 15 are not used by the transducers in this study, 

but their characteristics are given for the study of interference 
or of filter selection. 
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Channel 15 is not used by any of the basic instruments studied in earlier reports; 
however, the study of its behavior for different combinations of filters was 
desired by NASA. 

The frequency response in the table is based upon maximum deviation and a 
deviation ratio of 5; it is not necessarily compatible with the actual bandwidth 
of the corresponding transducers. In fact, for the yaw rate gyro (channel 6 ) 
roll rate gyro (channel 7) and for the Giannini pressure transducer (channel 11) 
the transducer bandwidth is certainly larger than the tabulated frequency respons 

Some nomenclature is recalled here: an ideal FM subcarrier oscillator has as 

input signal a variable voltage vi n (t) = S - -*- w^ n starting at time t = 0 , S“* being 
a scale factor (volts X sec/rad) and as output signal a voltage (or current) 

r t 

v out (t) = A cos < w o t + | w in ( tl ) dt ) 

J o 

where A = constant amplitude 

w Q /(2 TC) = f Q = center frequency of the S. C. O. 

2 JT f(t) = w (t) = w Q + wi n (t) - instantaneous pulsation of 

the output signal 

(w-w 0 )/(2 7T) = f - f Q = instantaneous deviation of output frequency 

Now Vj n (t) can be decomposed in its sine-cosine components by Fourier-series 
or Fourier -integral. 

One such component may be: 

v in, I;. ( t ) = - cos -C- t; 0 ^ t i 4 1 

with c £ w = a fixed amplitude 

_C1 : pulsation within the transmission band of the transducer 

Combination of the two preceding equations yields: 

v out, a <*) = A cos < w o 1 + - 77 -^ sin *) 
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Hence, one has the maximum instantaneous frequency - deviation 


The deviation ratio is defined as 


m = 


u w | /^ max 


(and the modulation index is sometimes defined as ] £ A w| / SX ) and m determines 
the side bands in the classical Bessel-function expansion. For the table above, 
m = 5 is used conventionally; however, the actual transducers can in some cases 
possess a max larger than that implied by m = 5 and ]A w| imposed by 
IRIG - Standards. 


The components to be discussed are, in order 

- phase- sensitive demodulators 

- subcarrier oscillators 

- mixer 

- transmitter -receiver 

- tape recorders 

- bandpass filters 

- discriminators 

- lowpass fiLters 

followed by a discussion of the entire chain viewed as a single block. The 
experimental test setup is described, with diagrams, in Appendix C. 

5. 2 Phase-Sensitive Demodulators (P . S, D, ) 

The output of the rate gyro output circuit is an AM-signal. A phase-sensitive 
demodulator is required to recover the envelope with correct algebraic sign. 
Its output is constrained to vary between 0 and 5 V in order to be compatible 
with the S. C„ O. that follows. 

The type used here is a ring demodulator as shown in the figure below. 


84 



3 

Instrument °f 


Output 



Phase -Sensitive Demodulator 


The static input-output relation for P. S 0 D 0 was found to be very linear within 
the normal range. Beyond the normal range there is first saturation and 
ultimately the output dips. Experimental data are shown in Figure 5-1. These 
characteristics provide a bound on over -ranging in channels 6, 7 and 8 and 
thereby help reduce channel interference. 

According to specifications, the P. S. D. sensitivity is bounded by 

maximum sensitivity: 5V out (D. C. ) for 20m V rms input 
minimum sensitivity: 5V out (D. C. ) for 200m V rms input 

Experimental data on sensitivities are shown in the table below. 


TABLE 

Steady State Gains for P. S. D. 


Number 

E in (mV) 

E . 
out 

(V. D. C. ) 

E in ( m V) r - m * s * 

A Volts out = K 

of P. S. D. 

r. m. s. 


for 5V output 

Volts in (r. m. s. ) 

1 

40 

4. 

377 

53. 3 

47. 0 

2 

38 

4. 

IS) 

00 

o 

53. 4 

46. 8 

3 

38 

4. 

350 

51. 4 

48. 6 

4 

41. 5 

3. 

173 

153. 0 

16. 3 

5 

38 

4. 

351 

51. 3 

48. 7 

6 

40 

3. 

107 

165. 0 

15. 1 
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The P. S. D. #4, 6 have low sensitivity, the others have high sensitivity (high 
gain). 

The transient response, here particularly the step response, turned out to be 
sensitive to direction of voltage swing at the output. 

While all these responses can be described with sufficient accuracy by a second- 
order transfer -function (including the ripple -filter) 

A V out = _^ f 5 ' 1 ) 

Yin (algebraic envelope) l+2\ s + 


(remembering the bias of 2. 5 volts at output) 

responses where the output voltage increases show some overshoot (with 
individual variations) while the responses where the output voltage decreases do 
not show any measurable overshoot (and this is so whether the input is in-phase 
or out-of-phase). 

Some experimental results are shown in Figures 5-2 and 5-3. Analog simulations 
are given in Figure 5-4 for the best fit (for which the agreement is quite good). 

The results of best fitting of experimental data are given in the table below. 

TABLE 

Dyna mic Parameters of P. S. D. 
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FIGURE 5-2 STEPRESPONSES OF PHASE - SENSITIVE 
DEMODULATOR #1 (UPPER TRACE) 

TIME SCALE 10 ms/DIVISION 
INPUT ON LOWER TRACE 



FIGURE 5-3 STEPRESPONSES OF PHASE - SENSITIVE 
DEMODULATOR #6 (UPPER TRACE) 

TIME SCALE 10 ms/DIVISION 
INPUT ON LOWER TRACE 
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Therefore, the lag associated with the P. S. D. is comparable to that of the rate 
gyros proper. One also observes some spread between individual P. S. D. 's 
(compare ^ for #1 and 6, output up). 

The parameters for underdamped response (output up) were obtained from 
values of overshoot and time of overshoot. Other points of the response were 
then compared with the theoretical values and found to agree within the accuracy 
of experimental data. 

E. g. , time of first full-scale passage experiment 15 msec 

P. S. D. #1 ( theoretical 14.7 msec 

Similar agreement was obtained for the (over) critically damped responses 
(output down). The fact that ^ depends on direction of output swing while w n 
varies very little can be explained by the observation that a resistance shift 
(in points A, B of the circuit diagram above) is used for the P. S. D. to obtain 
the bias of 2. 5V out. Therefore, the resistance in the up-and-down current 
paths is different, leading to different damping coefficients. 

NOTE 1: No detailed circuit-schematics were available for the P. S. D. , and 

therefore the present second order model is mainly based on 
experimental evidence. 

NOTE 2: Full scale output voltage from the gyro output circuit is approximately 

9. 5V r, m. s. On the other hand, maximum normal input to the P. S. D. 1 s 
varies between 0. 0513 V rms and 0. 165 V rms as shown in the table 
for steady state gains. Therefore, some signal reduction between 
gyros and P. S. D, 1 s is very likely. The amount of this reduction, 
however, is unknown since no complete schematics for this part were 
available. 


5. 3 Sub-Carrier Oscillators 

The S. C. O. 's are voltage controlled oscillators whose output has an instantaneous 
frequency deviation proportional to the input voltage, except for transient 
effects, (see Section 5. 1 for the ideal case). 

Steady State Characteristics ( n Static Gains' 1 ) 

The experimental data on static input-output relationships are shown in Figures 
5-5 through 5-10 as output-frequency (Hz) versus input volts D„ C. It is clear 
that these are very linear relationships, which moreover extend considerably 
beyond the nominal input range (0 to 5 V. D. C. ) without any sign of saturation. The 
edges (0V and 5V) also correspond accurately with the specified IRIG -values. 
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Figure 5-8, Static Gain of S. C. O. 10 (5. 4 kHz) 
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The very wide range of linearity is most clearly exhibited on Figure 5-7, the 
S. C. O, of channel #8, especially on the side of the positive input voltages. 

This figure also gives an idea of the variations in output frequency with voltage 
in. This means that adjacent channels could possibly interfere with each other, 
unless a limitation occurs earlier. For the gyros, the P. S. D. exhibits 
saturation as was explained in the preceding Section 5. 2. The pressure trans- 
ducers have a maximum supply voltage of 5 volts, imposing a bound. The 
accelerometers are also limited around their normal range. Therefore, no 
interference is expected in normal operation with standard filters (see also 
further discussion of interference in Section 5. 11). 

The differential static gains are listed in the table below. 

TABLE S. C. O. Static gains K w 4 (rad/ sec) out Kf A (Hz) out 

(A Volts)in A (Volts) in 


Channel Number 

K 

w 

K f 

6 

322. 

51.2 

7 

435. 

69. 2 

8 

565. 

90. 0 

10 

1020. 

162. 0 

11 

1380. 

220. 4 

12 

1980. 

315. 2 

13 

3740. 

435. 2 

14 

4150. 

660. 0 

15 

5650. 

900. 0 


One has obviously = 2/1 K£ ; comparing with the expression 

in Section 5. 1 one has also 

< K w)ideal = S 

Generally, K f = 0. 03 (F c ) (Hz) 

where F c is the center -frequency of the pertinent 
channel. 
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Transient Response (Step Response) 


A more realistic description than that of Section 5. 1 for the S. C. O. response 
to an input signal V- n ( t ), starting at t = 0, is given by 


V out (t) = A cos (w Q t + K v 


V in (t") g (t'-t") dt" dt') 


This relation assumes a linear, time -invariant dynamic response for the output 
frequency characterized by the normalized input response g(t) (with f g(t) dt = 1) 

g(t) exp (-st) dt, with G(o) = 1 

o 

(the Laplace transform of g (t )) . The idealized model assumed that G(s) = 1 for 
all s„ The present model assumes no changes in the output voltage envelope 
A; unless A is made to depend upon Vi n (t). However, the only dynamics of 
interest here are described by G(s), as a transfer function. 


A w out (s) 

V in (*) 


= K w G(s) 


( 5 - 2 ) 


Experimentally, V^ n (t) can be realized easily as a step (or as a rectangular 
wave with sufficiently low repetition rate to allow settling of the output). The 
output-frequency w out (t), however, is not directly measurable. A complicated 
way of obtaining w 0 ut (t), could be based on a correlation study of V out (t). On the 
other hand, the delay or lag associated with the S. C. O. is experimentally found 
to be quite small compared to other components in the T/M chain, especially 
the low-pass filter. Therefore, a simple approximation to w Qut (t) is sufficient, 
and this is obtained by measuring time-intervals between successive zero-crossings 
of V ou t(t) after the step-change in V^ n (t). The results appear as shown in the 
figure below (exaggerated for clarity). 



S. C. O. Step Response 


98 



The meaning of T ave is also explained in this figure. Typical experimental 
results are shown in Figures 5-11 a, b, c for channel #8. An average value 
at T ave has been used for rise *_nd fall in V^ n (t). Numerical results are shown 
in the table below, together with the experimental uncertainties. 


TABLE 

Settling Times for S. C. O. [s 


Channel Number 
and F c (kHz) 

T ave (msec) 
Experimental 

Uncertainty 
l^ T avel (msec) 

[l.75/ Fc (k H z)j<m s e c) 

Average Fit 

6 

1. 7 

1. 3 

0. 15 

1. 03 

7 

2. 3 

0. 65 

0. 11 

0. 76 

8 

3 0 0 

0. 66 

0. 08 

0. 58 

10 

5. 4 

0. 32 

0. 046 

0. 32 

11 

7. 35 

0. 25 

0. 034 

0. 24 

12 

10. 5 

0.166 

0. 024 

0. 166 

13 

14. 5 

0.112 

0. 026 

0. 120 

14 

22. 0 

0.091 

0. 011 

0. 080 

15 

30. 0 

0.067 

0. 008 

0. 058 


In fact ] A T ave J = %£! 

* c 

These results are also shown on Figure 5-12, indicating for the present 
equipment, channel #7 has unusually fast S. C. O. response (but well within the 
uncertainties and permitted tolerances). Figure 5-12 also indicates a good 
approximate average fit. 

T = 1. 75 

ave Fc 

also shown in the table above 


Inspection of the transient behavior shows that one can take for G(s) a second 
order expression 


G(s) = 


1 + 2 £ 


n 


w 2 n 


(5-3) 
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FIGURE 5-11 A 

STEPRESPONSE OF S. QO. 
#8 (3 KHz), TIME SCALE 
2 ms/DIVISION , INPUT 
IS 0 - 5 VDC SQUARE WAVE 


FIGURE 5-1 IB 

DETAIL OF FIGURE S-liA 
(RISE OF INPUT) 

TIME SCALE 0.2 ms/DIVISION 




FIGURE 5-1 1 C 

DETAIL OF FIGURE 5-llA 
(FALL OF INPUT) 

TIME SCALE 0.2 ms/DIVISION 
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with little or no discernible overshoot (with an accuracy not better than 5% of 
full scale). Since overshoot for ^ = 0. 7 is 4. 5%, a sufficiently accurate 

model is 


^ = 0. 7 for all channels 
and then w„ = hJ. 

II ryi 

1 ave 

(this is the point where the output first reaches full scale) 
So, if the average fit is used, one gets 

w n (rad/ sec) = 1.9 = 1, 900 F c ^ Hz ^ 

f r/ HZ ^ = °* 30 f c <HZ) = 300 F c ^ kHZ ^ 

These results are listed in the table below 


TABLE 


Frequencies of Transfer Functions for 


Channel No. 

S. C. O. 

w n (rad/ sec) 
Experimental 

w n (rad/ sec) 
Average Fit 

f n (Hz) 
A verag 

6 

2, 540 

3, 230 

510 

7 

5, 100 

4, 360 

690 

8 

5, 000 

5, 700 

900 

10 

10, 300 

10, 250 

1, 620 

11 

13, 200 

14, 000 

2, 205 

12 

19, 900 

20, 000 

3, 150 

13 

29, 500 

27, 600 

4, 350 

14 

36, 300 

41, 800 

6, 600 

15 

49, 200 

57, 000 

9, 000 


It is obvious from the last two tables that the dynamic lag due to the S. C. O. is 
indeed quite a small fraction of the lag due to the basic transducers. 
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5. 4 The Mixer 

The mixer takes the different signals from the S. C. O. 1 s of each channel and 
adds them to form a single composite output signal. 

^(output mixer) ” ^ ^(output signals of S. C. O. ) 

all channels 

Figures 5-13 and 5-14 show that the mixer has a nearly perfectly flat gain over 
the frequency band of interest and for normal input voltage amplitudes (steady 
state results). As for the dynamic behavior, excitation with a triangular wave 
of variable repetition frequency showed no lag, delay or distortion at the output 
within the frequency-band of interest (certainly between 0. 5 kHz and 100 kHz). 
Only at very low rates of 100 Hz (quite below the channel 6 frequency of 1, 700 Hz) 
did some distortion appear. 

Therefore, the mixer can be represented by a constant gain K (of 1. 0), same 
for all channels 6 through 15, without any dynamical distortion effects. 


5. 5 Transmitter and Receiver 

The transmitter proper and the receiver proper are very wideband and have a 
very linear output-input characteristic; they are considered to have perfect, 
instantaneous response (compared with the delays and lags of the other components) 
and, for this study, their transfer -function is taken to be unity. The radio link 
in vacuum and atmosphere depends, of course, on a series of widely variable 
and often uncontrollable factors, such as distance and altitude, atmospheric 
conditions, and its study was not a part of the present effort. 


5. 6 Tape Recorders 

For the channels studied in this program with center frequencies between 1, 700 
Hz and 30, 000 Hz, experiments showed that recording on tape and reading off 
tape introduced no discernible distortion, differential lag or delay in the signals 
from the various channels. The setup used to study this problem is shown in 
Figure 5-15. Conventional equipment was used. The counter was employed to 
determine accurately whether any differential change in the ZH(see Figure 5-15) 
could be observed, when the tape-machines were used, compared with the A t 
without tape-machines. Since no such effect was observed, the tape recorder can 
be viewed as an ideal component with transmission equal to unity. 


4 

■ 
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Figure 5-13, Mixer Gain Versus Frequency 
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5.7 


Bandpass Filters (B. P, ) 


For the EMR Model 229 Tunable discriminator, the B. P. filters are active 
two-stage filters. Each stage consists of a feedback amplifier with band -reject 
bridged-T feedback network. 

The conventional transfer -function of a bandpass -filter is defined by 


f bp( s ) 


E out 


E. 

m 


(s) 

W 


(5-4) 


with E Qut (s) and E^ n (s) being the Laplace -transforms of BP output signal 
e out (t) and BP input signal e^ n (t). However, in telemetry e^ n (t) is essentially 
a sum of output signals from S. C. O n f s (see Section 5. 3). 


e in< t ) = A k cos < w o, k* + K W) k / dt 'f 

^o o 

r t 

= \ A k cos ( w o,k t+ J A w in>k (t')dt') 


v in, k< t ")gk( t '- t ") dt ") 


with index k indicating the various IE IG- channels. 


Theoretically, each BP filter lets pass through only the correct k-component 
for the particular channel, and e Qut (t) has the form, for channel k. 


e (t) 
out * ' 


A , cos 
out, k 


< w o, k t + 


/’ 


out, k 


(f) dt') 


The desired transfer -function in this study is not Fgp(s) but the frequency- 
transfer function defined as 


F BP, k< s) 


^ W out,k < s ) 
A w in, k < s ) 


(5-5) 


for each channel k. 

A 

Fgp k( s ) is hard to measure directly (and would require extensive data processing, 
like correlation studies, etc. ), but fortunately two properties permit its indirect 
identification. 
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First Property: Bandpass to Low Pass Transformation 


If a BP filter with (conventional) TF-function Fgp(s) can be regarded as a 
displaced low pass (DLP) filter with TF-function Fggpfs), such that, for 

s = iw (i = V r -T) 

F bp (i w) = e' 1 ^ F DL p(i (w-w Q )) + e 1 ^ F DLp (i (w+w 0 )) 
then the impulse-responses are connected by 

f BP (t) = 1 *DLP W cos o 1 -< E) 

Therefore, the envelope of the response of BP to an amplitude -modulated signal 
(with centered car rier -frequency '^o : m(t) cos (w Q t -y- ) is the response of the 

2*r 

DLP to the envelope m(t), both being given by 



(t - V) m (t 1 ) dt r 


provided the BP filter is relatively narrow (i. e. , its pass band is relatively 
small compared to its center frequency w o or, equivalently, the time-constant 

2 rc 

of fDLP( t ) l ar g e compared to 2 7T ). This property is discussed in Chapter 7 
of Reference 8. w Q 

Second Property: Relationship Between Conventional and Frequency Transfer 

Functions of a BP-Filter 


Generalizing a result found by H. Salinger and reported in Appendix 13 of 
Reference 9, one can show that for a BP filter having the first property, the 
following simple result is valid. 


F BP ^ 


F PLP( S > 

F DLp( o) 


provided the maximum average frequency deviation in the input is very small 
relative to the bandwidth of the bandpass filter. Appendix A of this report 
derives this result and the conditions of its applicability. 
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Conclusion 


To determine Fgp(s) it is sufficient to measure (experimentally) the envelope 
of the transient response of the BP filter to an amplitude -modulated signal with 
centered carrier (say, with a step-change in amplitude). 

Each stage of the BP has the configuration shown below. The transfer function 
for such a stage is easily found to be 


E 0 (s) (I^Cj^s (1 + ^ 3^2 s ) (1 + R 4 C 3 s ) 

E^s) ” " (l + (R 1 C 1 )s) (1 + R 3 C 2 s + R 3 C 2 R 4 C 3 ) 


R 2 C 1 s 

1 "I - R ^ C ^ s 


( 1+ 


R 4 C 3 


-<$r~= 


1 + R 0 C-j s + R^C^R^C^s 


‘ 3 2 



Voltage amplifier, gain = 1 , 
High Input Impedance 
Low Input Impedance 


One stage of the active BP filter. 

In the vicinity of the passband center, defined by 
w 2 o (R 3 C 2 ) (R 4 C 3 ) = 1 

one has the usual RLC -tuned behavior, (normalized): 


1 


1 + e o (— + 

V w Q 



w h ere 

e o = ( r 4 c 3 /R 3 c 2 ) 


The resistances R^, R 4 are used for tuning w Q , while C 2 > C 3 are employed for 
bandwidth choice,, The normalized low pass equivalent is of the form 

1 

1 + T s 


with T being close to 2 0 o /w o . 

Since the BP uses two such stages and since the experimental evidence shows 
no overshoot in the envelope step response, one may assume for practical purposes 
nearly zero or zero staggering (see again Chapter 7 of Reference 8 ), Therefore, 
one has for the displaced low pass filter 

F DLP = 1 (5-6) 

F DLp( o) (1 + T s ) 2 

with step response 1 = exp (-t ) (1 + t ) 

T T 


corresponding with a critically damped second order system 


(i. e. , ^ = 1 . 0 in 


1 

1 + 2 \ _i_ 



) 


(5-7) 


According to the manufacturer, the BP bandwidth (i. e. , the edges) are defined 
as points approximately 2. 5 db (or a factor of 1/1. 334 = 0. 75) down from center 
frequency. For the DLP above this corresponds with 

Wg T = 0. 578 (subscript E denotes edge) 

Lr 

Since We t -p = W v - W 0 , this gives for a + 7. 5% BP filter 

0. 578 
0. 075 w Q 




(+ 7. 5%) 


-1— i , UJT 

0 . 578 


Z'Tt 0. 075 x F,. 
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and for a 4* 15% BP filter 


■(+ 15 %) 


0. 578 

2 TV 0. 15 F 


0. 578 
0. 15 v/ 


where F c = 


2tt 


is, as before, the center -frequency of the pertinent channel. 


Particularly, for channel #6, F c = 1, 700 Hz. 
w Q = 10, 700 rad/ sec 


and 


(+ 7. 5%, ch. #6) 
'(+ 15%, ch. #6) 


= 0. 72 msec. 


0. 36 msec. 


_ O 1 

As a verification, (1 + Ts) is down by a factor 


3. 376 


or 10. 57 db (more than 


10 db imposed by the specifications) for 
w 20 


W E, LP 


7. 5 


corresponding with a 20% shift in carrier -frequency for the BP 

filter. The DLP characteristic (1 + Ts)"* or rather the RLC-BP characteristic 


(1 4“ 0 (_? + Wq )) ^ from which it was derived, also shows a satisfactory 

° — s 


numerical agreement with the experimental frequency response of the BP 
exhibited in Figure 5-16. 

In the table below, for channel #6, we now present a comparison between 
theoretical predictions for a step response using the T-values derived above 
and the experimental evidence, which is also shown in Figure 5-17. 


TABLE 

Channel #6 Step-Response 


Lab. Experiments 

Theoretical Value 

+ 7. 5% filter: (average of 4 measurements) 
Time 1/2 Full Scale is reached 1 . 4 Jh 0.3 

1.7 T = 1.23 msec. 

Time response is fully settled 3. 7 + 0. 3 
(take 97. 5%, midway error) 

5. 5 T = 3. 98 msec. 

+ 15% filter: (average of measurements) 
Time 0. 5 full scale 0. 8 _+ 0. 3 

1.7 T = 0. 615 msec. 

Time fully settled 1. 8 +_ 0. 3 

5. 5 T = 1. 99 msec. 
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Figure 5-16, Frequency Response Bandpass 
Filter 
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FIGURE 5-17A STEPRESPONSE OF 7 . 5 $ BP FILTER OF CHANNEL #6 
TIME SCALE 1 ms/DIVISION 

LOWER TRACE: INPUT TO BP (CENTERED CARRIER 

WITH STEP IN ENVELOPE) 

UPPER TRACE: BP RESPONSE 


: 1 

1 

r r ^ 1 

1 |~ j 

' 1 1 


"1 



J . i 1 4 

4 

* f • • 

^ "V y ^ - 


* f *** oX »V. — s. U.S. «y. 


FIGURE 5 -JJB STEPRESPONSE OF 1 5# BP FILTER OF CHANNEL #6 

SAME TIME SCALE AND COMMENTS AS FIGURE 5-17A 


113 




The agreement is quite satisfactory (within the accuracy of the measurements). 
The smallness of T shows, moreover, that the bandpass filter introduces only 
a small amount of lag and distortion when compared to the low pass filter (if 
the standard filters are used, of course). 

It is interesting to display the non-linear dependence of the BP output on the 
magnitude of the frequency swing at its input (in the signal coming from the 
S. C. O. ) a feature discussed in Appendix A. This is done in Figure 5-18 . 

Channel #6 was isolated there to avoid interference from its neighbors. The 
extra-wide 125 Hz (constant delay) LP filter was used, instead of the standard 
25 Hz LP in order to minimize the distortion in the LP. Trace A of Figure 5-18 
corresponds with a small swing of 2V to 3V (centered around the middle of the 
range at 2. 5V) at the input of the S. C. O. , i. e. , a small swing A w^ n at the 
input to the BP. In fact, with the notations of this section and of Appendix A 

' AW ^ = 1 - 0.2 

W E, LP 5 

Then the small signal approximation of Appendix A is valid, i. e„ , the "second 
property" of this section applies, and the response of Trace A agrees with the 
normalized step-response derived above as | L w* n | (1 - exp (- 1 — )(1 + )) 

(also shown in Figure 5-17), except, of course, for some additional delay 
introduced by the LP. 

For trace B and a fortiori for trace C, the frequency swing is large 

I A w in ) = j 

W E, LP 

and non-linear effects appear clearly. Just as for the ideal case mentioned in 
the Appendix A and discussed in Reference 9, a marked overshoot appears 
(symmetric for this type of input swing which is itself symmetric about the center 
frequency). The experimental overshoot of trace A is approximately 10% and 
some rather lengthy calculations indicate that this overshoot is quite compatible 
with the theory, given the uncertainties in the quantities measured. 
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FIGURE 5R.8 DEPENDENCE OF B. P. ON MAGNITUDE OF ITS 
INPUT FREQUENCY SWING (i.e., MAGNITUDE 
OF INPUT VOITAGE SWING TO S. C. 0. ) 

ISOLATED CHANNEL /6; 7-5$ B. P. , 125 Hz 
CONSTANT DELAY L. P. 


TRACE 

A 

INPUT SWING 

2V 

to 3V 

TRACE 

B 

INPUT SWING 

IV 

to 4V 

TRACE 

C 

INPUT SWING 

OV 

to 5V 

TIME SCALE 

10 ms/DIVISION 
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5. 8 


Discriminators 


EMR 229 employs a phase lock-loop discriminator. The dynamic effects are 
described mainly by the loop filter and the total loop gain. A linearized model 
of the discriminator is shown (compare with Reference 7, Chapter 8). 



- * 
, i 


-n" c- 


(K K ) has dimension ohm 
' p v' 


Ka 


is dimensionless 


volt 

volt 


Fig. Linearized Model of the Discriminator 

l + R Ci s 


With these notations one has 

1 


LF (Cj+C^Js 

v out ( s ) 1 + RCj S 

L w in (s) 


1 + RC 


C 2 


1 Ci+C 2 


1+^2 ^*^ 1^2 

K v (1 + < EC i + irr) s+ “ 


C ' lC, 


K A K P K 


s2 ) 


(5-8) 


(5-9) 


P v 
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The loop gain is generally sufficiently large such that 


K A KpKy 1 

Ci +C 2 RCj 


Therefore, the resulting transfer function takes the standard form 


K a Kp 


K, 


R C, C; 


-i-L = R Cj 


w. 


n 


A w in (s) 




1 + 2 


w. 


n 


1 + 2 (“ - — ( 1 + 

1 W 


Cj + c 2 
C ! R KpK v K A 


) 


(5-10) 


A w i n ( t ) i s measured from the center frequency of the corresponding channel. 


The phase error in the loop (see diagram above) is given by 

s , C 1 + c 2 


^ e( s ) 


w (s) 

error 


4 w in (s) 


K A Kp K v 


4 ^ in <s) 


1 + 2 ^ — — 


w n 


w 2 n 


The static phase error at band edge A w ^ n = is given by the 

manufacturer as less than 2° = 0. 035 rad. 


Therefore, 


0. 035 C 1 + c 2 

W E,LP K A K p K v 


1 

Loop Gain (rad/ sec) 


For Channel #6, as an example 

K a K K 27i- 0. 075 F_ 

A P v x £ 

Ci + C 2 7 0. 035 


800 

= 2. 3 x 10 4 rad/sec 

0. 035 


During normal operation, 

90° = ZL rad 

2 


the dynamic phase error should be less than 
(no lock-loss). 
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For infinite loop gain, this amounts to the following condition, derived from 
the transient response of g(t) 


a I , . - Arc sin f ) 

^ Win 1 max, normal exp(- \ K 2 _) < ' rr /2 fif ("<1) 

w n v' i - r iA 


= w-n T -d “ (800 rad/ sec for channel #6) 

m I max, normal 


With ( A 

one finds (if ^ < 1) 


v n > 0. 64 W E> LP exp (- 


( — - Arc iinp) 




One has also 
2 \ w n = 


K A K p K v 

Cl + c. 


C l + C 2 

C 2 


C 1 +C 2 


Loop Gain ^ Loop Gain 


The preceding two inequalities indicate that the discriminator has proper 
frequencies of an order of magnitude higher than those of the BP filter (+ 7. 5%) 
and thus responds about ten times faster, i. e. , even faster than the S. C. O. 
Therefore, within the accuracy of the experimentation, the discriminator 
dynamics may be neglected entirely . (Transfer function is unity). 


5. 9 The Low Pass Filters (LP-F) 

The output from the discriminator still contains spurious components such as 
a ripple at twice the frequency of the subcarrier. The purpose of the LP-F is 
to remove such unwanted components, while distorting the true signal as little 
as possible. 


Basically, the 210 C-01 LPF’s in the EMR 229 have the schematics shown 
below. v 



v. 


Low Pass (Output) Filter 
Diagram 
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The transfer function is a third-order function 


V out( s > 


V in (s) r i + r 2 


1 + a^s + a 3 s 2 + a^s 3 


a, = 


R 1 R 2 
R l +R 2 


Cj + (r 3 + R 4 ) c 3 + 


c 3 

RltEj 3 


with 


= 


C 3 C 2 R 3 R 4 + C 3 C 1 p — < r 2 r 3 + R 3 R 4 + R 4 R 2^ 


C 1 Cz C 3 R 1 R 2 R 3 R 4 


R l +R 2 


Through an appropriate choice of component values, two basic types of filters 
are realized by the same diagram. 

1) Constant Amplitude or Butterworth Filters, with general 
transfer -function (normalized). 

By (s) = A (5- 11 ) 

1 + 2 T bU s + 2 T 2 bU s 2 + T 3 bU s 3 

1 

(1+T bu s) (1+T bu s + t 2 bU s 2 ) 

2) Constant Delay or Linear Phase or Bessel Filters , with general 

^normalized) transfer-function. 


B e (s) 


15 

15 + 15 T be , + 6T 2 be s 2 + T 3 be S 3 


(5-12) 
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The preceding filters are discussed in References 8 and 10 and useful transient 
(including step) responses are shown in Reference 11. 

Determination of Parameter Values 


Everything is reduced to an estimation of two time -constants T BU’ t be 
( provided the real filter behavior approaches the theoretical one). 

The filters were subjected to rectangular inputs of sufficiently low repetition 
rate to constitute nearly perfect step inputs. 

For the Butterworth filter, the manufacturer specifies the cutoff-frequency 

w^ = 2 rr by an attenuation of 0. 5 db (a factor of -J ). This gives 

co co 1. 059 


T BU w co 


0.704 


The following table gives a comparison of theoretical and experimental results, 
based on the preceding relation. 


TABLE 

Step Response of Butterworth Filter 
Channel #6 ^ co = 25 Hz (standard filter) 
w co = 157 rad/ sec 


T BU “ 4. 5 msec 



Theoretical 

Experimental 

Overshoot 

8.15% i 

(8. 1 + 1%) 

Time of Overshoot 
Time full-scale is 

4. 8 T-gjj= 21. 6 msec 

(22 + 1) msec 

reached first 
Time full-scale is 
reached 

3.8 Tgy=17. 1 msec 

(17 + 1 ) msec 

second time 
Time full-scale is 
reached third time 

7. 2 T bu = 32. 4msec 

( 30 + 2) msec 

(essentially settled) 

10. 5 Tgu= 47. 4 msec 

(46 +3) msec 

1 
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The agreement is excellent, and the general formula above for Tgy is 
valid for all channels. 


For the Bessel-filter, the manufacturer 1 s amplitude response curves give 

3. 35 db at w (cutoff), 12. 3 db at twice cutoff and 20. 9 db at three times cutoff, 
c o 

This leads to 


BE 


1 . 8 


The table below gives again a comparison of theoretical and experimental 
results. 


TABLE 

Step Response of Bessel-Filter 


Channel #6 


f = 25 Hz 
co 


w =157 rad / sec 
co 



T^ =11.5 msec 


Theoretical 

Experimental 

Time half-scale is 
reached 

0. 96 = 11. 0 msec 

(9. 6 + 1.5) msec 

Time of settling 
(within accuracy) 

2. 3 Tgg = 26. 4 msec 

(26.6 + 1.5) msec 

Overshoot 

0. 75% 

Not noticeable 

1 


Other experiments for other cutoff-frequencies gave similar satisfactory 
agreements such that the general formula above for Tgg is an acceptable 
formula. Experimental step responses for both Butterworth and Bessel filters 
are shown in Figure 5-19 A Sc B, and the analog computer verification is 
exhibited in Figure 5-20. 
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FIGURE 5-19A STEPRESPONSE OF 25 Hz CONSTANT AMPLITUDE 

LOW PASS FILTER . TIME SCALE 20 ms/DIVISION 



FIGURE 5-19B STEPRESPONSE OF 25 Hz CONSTANT DELAY LOW PASS FILTER 
TIME SCALE 20 ms/DIVISION 
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5. 10 Complete Telemetry Chain 


After the individual components were tested and their characteristics determined, 
the transmission of signals through the entire cascade of T/M components was 
studied. For analog computer verification, the separate models of each 
individual successive component (P. S. D. , S. C. O, B. P. , L. P. ) were simply- 
cascaded. 

Because of interaction of successive stages, the entire chain might a 
priori have a transfer function different from the product of the individual 
transfer-functions. Such interaction was observed for the pressure-transducers. 
However, for the T/M chain, the interaction turned out to be negligible and 
the cascading of individual component models is a valid procedure. Moreover, 
for standard filters, the LP filter is by far slower than the other components, 
except the P. S. D. (where used) and therefore the entire chain will behave very 
much as the LP-F with some relatively small additional lags superimposed, or, 
in the case of channels #6, 7 and 8 the entire chain will behave very nearly as 
the LP-F combined with P. S. D. (if used). 

Experimental results and the analog simulation corroborating the statements made 
above are shown: 

a) For trapezoidal inputs in Figures 5-21 A 8c B, 5-22. 

The agreement is very good. Figure 5-23 shows the response of the 
LP alone to the same input, proving indeed that the LP is by far the 
major contributor of dynamic distortion in the T/M chain. Results for 
other channels can be simply obtained by correct time scaling. 

b) For triangular inputs in Figures 5-24 A & B, 5-25 and 5-26. The 
same comments apply as in the preceding case of trapezoids. 

c) For rectangular (step) inputs, in Figure 5-20 seen before and also 

in Figures 5-27 A & B, 5-28 and 5-29, with the same comments as above. 

d) For rectangular (step) inputs, but now with P. S. D. included, see 
Figures 5-30 A & B, 5-31 and 5-32. It is obvious that indeed, as 
mentioned before, the P. S. D. constituted a noticeable part of the dynamic 
lag, comparable to that of the LP itself. 

NOTE: Experimental data for other channels were obtained, but 

they offer nothing new, since these results can be obtained simply by 
rescaling times, inversely proportional to the center -frequency. 

Peculiar effects in the total T/M chain are now discussed in the next sub- 
sections 5. 1 1 and 5. 12. 
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FIGURE 5 -21 A RESPONSE OF CHANNEL #6 TO TRAPEZOIDAL 
INPUTS (P.S.D.. EXCLUDED) 

7.5$ B.P. , 25 Hz CONSTANT AMPLITUDE L. P. 
TIME SCALE 20 ms/DIVISION 



FIGURE 5-21B SAME AS FIGURE 5 -21 A EXCEPT THAT CONSTANT 
DELAY L. P. WAS USED 
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Figure 5-22, Responses of Channel 6 to Trapezoidal Inputs 
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FIGURE 5 -24A 


RESPONSE OF CHANNEL $6 (P.S.D. EXCLUDED) 
TO TRIANGULAR INPUTS 

7.5$ B.P., 25 Hz CONSTANT AMPLITUDE L. P. 
TIME SCALE 20 ms/E>IVISION 



FIGURE 5 -24b SAME AS FIGURE 5-24A, EXCEPT TEAT CONSTANT 
DELAY L. P. WAS USED 
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FIGURE 5 - 2 7-A STEPRESPONSE OF CHANNEL #6 (EXCLUDING P.S.D. ) 
7.5 $ B.P. , 25 Hz CONSTANT AMPLITUDE L. P. 

TIME SCALE 20 ms/DIVISION 



FIGURE 5-27B STEPRESPONSE OF CHANNEL #6 (EXCLUDING P.S.D. ) 
7,5$ B. P., 25 Hz CONSTANT DELAY L. P. 

TIME SCALE 10 ms/DIVISION 
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Figure 5-28, Step Responses of Channel 6 
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Ordinates: Normalized Responses 


esponses of Channel 6 
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Time (msec 










FIGURE 5 -3 OA STEP INPUT OF CHANNEL #6, INCLUDING P.S.D. #1 
7.5# B. P., 25 Hz CONSTANT AMPLITUDE L. P. 
TIME SCALE 20 ms/DTVISION 

TRACE 1 INPUT TO P. S. D. 

TRACE 2 OUTPUT OF P. S. D. 

TRACE 3 OUTPUT FROM L. P. 



Figure 5 -3 OB SAME AS FIGURE 5 - 3 0 A EXCEPT THAT CONSTANT 
DELAY FILTER WAS USED. 
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5.11 Choice^ of Filters, and Channel Interference 


The discussion of Section 5. 1 mentioned that the standard frequency response 
of each channel (as listed in column 5 of table 5. 1) was based upon the assumption 
of a deviation ratio of 5. It was also pointed out that this standard frequency 
response (such as 25 Hz for channel #6) is not necessarily compatible with the 
transducers used (in this case a yaw rate gyro of natural frequency 26 Hz, which 
can certainly transmit signals higher than 26 Hz without excessive attenuation). 
Therefore, the user might think of using non-standard filters (both BP and LP) 
to fit the transducer bandwidth more closely. To illustrate this point. Figures 
5-33 A through F, 5-34 A through F show various step responses of channel #6 
(with all other channels under power, and zero inputs) for 7. 5% and 15% BP and 
for LP going from 15 Hz to 125 Hz cut-off frequency (NOTE: 25 Hz is the standard 

JLP). It is clear that the use of wider band LP permits a faster response and a 
Smaller (dynamic lag and distortion, but the figures show that there is a limitation 
to that procedure. Indeed, as the LP becomes much wider than the standard LP, 
a ripple appears (see Figures 5-34D and 5-34C), which becomes very objectionable 
for very wide band LP (see the 125 Hz LP filters). It will be shown further that 
this ripple is due to interference from adjacent channels. Pictures such as 
Figure 5-33 and 5-34 permit the user to make a trade-off between desired speed 
of response (lack of dynamic distortion) and interference -ripple (or noise). 

The user’s choice would, of course, depend on his particular preferences and 
on the expected types of signals. Comparison of the pictures also show that, for 
the same LP bandwidth, the ripple increases with bandwidth of the BP. 

The Figures 5-35 A & B, giving step responses of the same channel #6, but now 
isolated (i. e. , all other channels are completely removed) show no ripple for 
a 125 Hz LP (to be compared with Figure 5-34 F) and only very little ripple 
for the extremely wide band filter with 500 Hz cut-off frequency. This proves 
that the ripple discussed here is indeed due to interference from other channels. 
These pictures of Figure 5-35 further indicate, that if the user has to transmit 
only a few signals, the response can be speeded up by the use of wider LP 
filters, provided the channels used (and powered) are spaced farther apart. 

In order to show the influence and importance of the ground equipment used, the 
pictures in Figure 5-36 A & B, give again the step responses of channel #6, but 
now with a different discriminator (the EMR Model 210 Fixed Discriminator). 
Comparison with Figures 5-33 shows less ripple for the fixed discriminator, 
due to the fact that the BP filters differ from the ones of the tunable discriminator 
used in all other experiments. (The fixed discriminator BP has a much sharper 
cut-off slope than the tunable BP). 

Finally, Figure 5-37 A through D proves that among all channels, it is the 
nearest neighbors which contribute most or all of the interference ripple. 
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FIGURE 5-33 

STEP RESPONSE OF CHANNEL 
6, 7. 5% B. P. , CONSTANT 
AMPLITUDE L. P. , (NO P. S. D. ) 
TIME SCALES 10 ms /DIVISION 

A) 15 Hz L. P. 



B) 25 Hz L. P. 
(STANDARD) 


C) 50 Hz L. P. 
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FIGURE 5-33 CONTINUED 


D) 75 Hz L. P. 




E) 100 Hz L. P. 


F) 125 Hz L. P. 


i 

■ _ 





FIGURE 5-34 



STEP RESPONSE OF CHANNEL 
#6, 15% B„ P. CONSTANT 
AMPLITUDE L. P. (NO P. S. D. ) 
TIME SCALES 10 ms /DIVISION 

A) 15 Hz L. P. 



B) 25 Hz L. P. 
(STANDARD) 
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FIGURE 5-34C0BTINUED 


D) 75 Hz L. P. 


E) 100 Hz L. P. 


F) 125 Hz I>„ P. 
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FIGURE 5-35 



I 


STEPRESPONSE OF THE ISOLATED 
CHANNEL jfS 

15 io B. P. , CONSTANT 
AMPLITUDE L. P. 

(NO P. S. D. ,) 

TIME SCALES: 10 ms/DIVISION 
A) 125 Hz L. P. 
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FIGURE 5-36 

STEPRESPONSE OF CHANNEL 
jj6, 7.5 B. P., CONSTANT 
AMPLITUDE L. P., (NO 
P.S.D. ) 

TIME SCALES: 10 ms/DIVISION 

FIXED DISCRIMINATOR 


A) 25 Hz L. P. 


B) 110 Hz L. P. 





FIGURE 5-37 

STEPRESPONSE OF CHANNEL 

#6 

TIME SCALES: 20 ms/DIVISION 

7.5$ B. P, CONSTANT AMPLITUDE 

L. P. 

A) 25 Hz L.P. (STANDARD) 

CHANNEL #5: OV D. C. 
CHANNEL #7: 5V D.C. 

(BEST CASE) 



B) 25 Hz L. P. (STANDARD) 

CHANNEL #5: 5V D. C. 
CHANNEL #7: OV D. C. 
(WORST CASE) 
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FIGURE 5-37 CONTINUED 


C) 125 Hz L. P. (WIDE) 

CHANNEL 5 : OV D. C. 

CHANNEL 7 : 5 V D . C. 

(BEST CASE) 



D) 125 Hz L. P. (WILE) 

CHANNEL 5: 5V D. C. 

CHANNEL 7: OV D. C. 

(WORST CASE) 
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This interference is easily explained theoretically. Consider the channel #6 
and its two neighbors #5 and #7. 


Channel 

Lower Edge 
(Hz) 

Center 

(Hz) 

Upper Edge 
(Hz) 

5 

1202 

1300 

1399 

6 

1572 

1700 

1828 

7 

2127 

2300 

2473 


In Figure 5-37 D (worst case with wide LP Filter) one observes on one side of 
the square wave a ripple with frequency of approximately 190 Hz 
(- 1572 - 1399 = 173) and on the other side a ripple at nearly 290 Hz 
(= 2127 - 1828 = 299). With the figure below one has (neglecting phase shifts 
in the BP, which are immaterial here) 



Constant | 
Signal 


Constant 


SCO 

Neighboring 

N 

m 2 (gain 
at Wj+ A Wj ) 


BP 

Channel 

?\ 

1 

? 


V2sin(w^+A W})t 




an output signal from the BP given by 

s (t) = (g(wj ) m j VJ ) sin w^ t + (g(w^ + A w^Jir^v^) sin (wj Wj)t 
Because of the selectivity of the BP, 

| g (wj + i Wj)|« | g (Wj) I 

must be satisfied. 
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Therefore, one has, since 
g (w x + A Wj) m 2 V 2 

g (wj) mj 


O (mj), Vj = 0 (V 2 ) 
£ j le I << i 


We can restrict our attention to the input signal 


si(t) = sin (wjt) + £ sin (wj + A w-^t 


to the discriminator. 


For neighboring channels 


1 < 1 , but not necessarily small. 

W 1 


Many discriminators, and particularly, the ones used in this study, basically 
detect phase differences or frequency differences by changes in zero-crossings. 
For the EMR 229, e. g. , phase detection is based on comparison of pulsed 
signals from a limiter and from a V, C„ O. (multi- vibrator ) in the feedback 
loop (See Section 5. 8). 


If £ = 0, zero crossings t n of si(t) are given by t n = nff /v/\ (n = 0, 1,2,...) 
and spacing t n + 1 - t n = *nr /w^. 

For £ non-zero, but small relative to unity, zero crossings are shifted by an 
amount A t n given by s^ (t n + A t n ) = 0 

Using the fact that |£| « 1 and hence t n j<< 2 rr , this leads to 

wi A t n + £ Sin ( - Wl nTT ) = 0 
W 1 


Hence the new spacing between zero crossings becomes 


Wl + ^ 4 n+l 

- (t 

n 

+ A 

4 n> : 

rr 

Z £ 

sin 

A w 

" ~i ' 

W 1 

w 


rr 

~T~ 


cos 


A wj 

W 1 


(n + j- ) rr 


T47 


Since the signal (frequency shift) given as output by the discriminator is 
inversely proportional to the length between zero crossings, one gets 


f (output discriminator) 
f (same for £. = 0) 


or, for running time t ( ft ^ ^t) 

w ! 


TC /wj 

TEL - 2 ^ — sin 4LLL — _ cos L (n+T) 7T 

wj wj Wj 2 W| 2 


f (t) 

f ( £ =0) 


= 1 + (TA sin ^ !■ TL ) cos (A w 1 t + L L 


Aw, 





One obtains indeed a ripple with frequency 


A w x 

— and with amplitude 

2 TT 


i_i_ sin A Wi zn 
7T w 1 2 

Due to the presence of £ , the ripple amplitude depends on the bandwidth of 

the BP filters (for narrow BP, £. decreases and there is less ripple amplitude). 
Moreover, the LP filter following the discriminator also helps in suppressing 
the ripple. Those statements are confirmed by inspection of the Figures 5-33 
and 5-34, showing the influence of BP and LP respectively. 


5. 12 Power - Inter ruption 

In order to observe the effect of power interruptions on the signal transmission, 
several tests were run under a variety of conditions. Some results are shown 
in Figure 5-38 A through E. Inspection of these pictures demonstrates that 
power interruptions can give rise to signals at the discriminator output resembling 
normally transmitted signals. One possibility of distinguishing between genuine 
transmitted signals and power -interruption signals exists in the fact that the 
latter category would show up as a simultaneous disturbance in all channels. 
However, in the case of catastrophic events sensed by the transducers, such 
simultaneity may also arise in regular signals not due to power interruption. 

Figures 5-39 and 5-40 give the output frequency of the SCO’s in channels 6 and 15 
as a function of power supply voltage. Clearly there is a wide range of supply 
variations (from 22V to 37V) where the output frequency is indeed independent 
of supply voltage, indicating the good quality of the SCO's. 
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FIGURE 5-38.A 

POWER FUIEE FROM 32V TO GROUND 
TIME SCAIE: 1 ms/DIVISION 


FIGURE 5-38B 

TRACE #1 CHANNEL #15 

TRACE #2 CHANNEL #6 

TRACE #3 INPUT SIGNAL 0-4V D. C. 

TRACE #4 POWER SUPPLY 24 V - 32V 

TIME SCALE: 20 ms/DIVISION 


FIGURE 5-38C 

TRACE 1 CHANNEL #15 

TRACE 2 CHANNEL #6 

TRACE 3 INPUT SIGNAL AT 86 Hz 

REPETITION RATE 

TRACE k POWER INTERRUPTION 32 V 

TO GROUND 


TIME SCALE: 5 ms/DIVISION 






FIGURE 5-38D 

TRACE 1 CHANNEL #15 
TRACE 2 CHARNEL #6 
TRACE 3 INPUT SIGNAL AT 
8 6 Hz REPETITION 
RATE 

TRACE 4 POWER INTERRUPTION 
32V to 2kV 

TIME SCALE 5ms/DIVISION 



FIGURE 5-38E 

TRACE 1 CHANNEL #15 
TRACE 2 CHANNEL #6 
TRACE 3 INPUT SIGNAL AT 

7.8 Hz REPETITION 
RATE 

TRACE k POWER INTERRUPTION 
FROM 32V to 2^V 

TIME SCALE gO ms/DIVISION 
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Figure 5-39 

V. C. O. Channel #6 Power Supply Variations 


1.6 


In 


\.s 





Power 






Figure 5-40, S. C. O. Channel 15 Power Supply Variations 







6. THE INVERSION PROBLEM 


6. 1 Introduction 

In previous sections 2, 3, and 4 analytical models were derived for the three 
categories of instruments; rate gyros, accelerometers, and pressure trans- 
ducers, including appendages such as mounting brackets and connecting tubings. 
These models permit direct calculation of the instrument outputs for given input 
signals, say, by means of analog simulation. However, in order to find out what 
really happens in the vehicle carrying the instruments, the inverse problem of 
calculat i ng input signals for known Jr ecorded) output signals is what matters. 

This report describes an inversion process for the above-mentioned instruments, 
under the following assumptions and rules. 

A. The recorded output data are available as equidistant , sampled quantities 
in engineering units, presented in tables or as graphs. (These data 
come from continuous traces having finite bandwidth). 

B. All instrument models are, where necessary, simplified to single-input, 
linear time -invariant systems (zero order model). This means that 
parasitic inputs such as accelerations for the gyro and pressure 
transducer, and transverse acceleration for the accelerometer, are 
neglected. The neglected secondary effects are small. 

C. The inversion is performed on a digital computer. It is an optimum 
inverted filtering, in a minimax sense. 

D. The reconstructed input signals are obtained as equidistant/ sampled 
quantities in engineering units, displayed in tabular form or as computer 
plots. Error bounds are given. The inversion process can be used to 
study the uncertainty with which the reconstruction is found, and the 
sensitivity of the results to changes in the system parameter values. It 
can also yield information on statistical effects by means of so-called 
Monte Carlo simulations. 

The testing of instruments, coupled with theoretical analysis, has resulted in 
the following typical (zero order) transfer -functions: (x is input signal , 

y = output signal, s = transform variable). 

Rate Gyros 


y(s) 

X(s) 


K 


1 + 2 £ 

n 


s 




(6-1) 
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The bracket effects are negligible for the main input. The electro-mechanical 
output device is considerably faster than the mechanical response of the 
gyro; moreover, the device forms an entity with the phase- sensitive demodulator 
(not included in the inversion task statement). Therefore (6-1) reflects the 
dominant part played by the mechanical response only. 

Accelerometers 


y(s) 

x(s) 


K 


1+2 f 


n w. 


n 


The bracket-effects are negligible for the main input. The damping 
to decrease as inputs vary faster (a non-linear effect). 


Pressure Transducers 


Giannini 


500 


lb 


in fc 


transducer 


(6-2) 


f tends 
' n 


y(s) 

x(s) 


K 


(1+2 f 

n 


s 

w r 


w 


) (l+2f -5- +. 


or 


(6-3A) 


The parameters f , w n represent the basic Bourdon tube, while ^ Qr , w or 
represent the connecting orifice. However, this interpretation is somewhat 
artificial, since both elements interact and (6-3A) is to be regarded as an 
entity which should not be split, in physical interpretations. The damping £ 
tends to increase as the inputs vary faster (a non-linear effect). 


Bourns 800 


lb 

. 2 
in 


y(s) 

x(s) 


transducer 


K 


1 + 2 


s 

a w a 


w^ 


d+2l — +^-)(l + 2l^-+JL_ ) 

n w n w 2 n w 2 b 


(6-3B) 
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The parameters f n w n represent the basic Bourdon tube; while f a , w a 
and f w-^ model a low frequency approximation of the T -tubing and orifice. 
Again, there is interaction and (6-3B) should also be regarded as a single 
unit. Of course, for purely mathematical purposes splitting is allowed. 

The preceding results show that all instruments (as was obviously expected) 
smooth out the input wave form x(t) into more gradually varying output waveforms 
y(t). Inversion therefore requires reconstruction of a rougher signal shape 
from smoothed data, and is in fact equivalent with differentiation of recorded 
signals. All the problems attending diff erentiation do show up in the inversion 
process . 


Observe that, with one exception, all transfer functions to be inverted are, 
second order systems. The one exception is the factor 1 + 2 f , 


w a 


in the T-tubing of the Bourns transducer, which, when inverted, becomes 
a classical, smoothing second order dynamic system. 


The next sections discuss the basic inversion relations, the error-analysis, 
interpolation and shifting, the cascading, the inversion of the numerator, etc. 


6. 2 Theory 

6.2,1 Basic Inversion Relations of the Second Order System 
Consider the normalized system (K = 1 ) 


y (s) 

x(s) 


1 + 2 


corresponding with the relationship 

y" (t) 


y (t) + \ y'(t) + 

W J 


x(t) 


n n 


n 


(6-4) 


The output data y(t) are available as N equidistant samples (directly measured 
or obtained by interpolation; see section 6.2.2). 


Yi = y (t Q ); y 2 = y (t o + T);. . . ; y^ = y (t Q + (k-i)T);. . . 

... y N = y + (n- 1 ) 1 ') 


1 
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We approximate the continuous second order operator (6-4) by a discrete 
second order operator: 


x k = *yk+l + {* Y k + r y k- 1 


(6-5) 


where x^. = x (t Q + (k-l)T) is the input sample to be reconstructed, at the 
same sampling instant as y^.. If y(t) is decomposed in its harmonic components 

e iwt , one finds the conditions 

^ TT 


* e iwt +/3 + Y e " iwt 


1 + i 


w 


w 


W r 


and, for zero-frequency accuracy, 
+ /3 + =1 

The foregoing conditions lead to: 

(^) 2 

* *r - 


1 - cos (w T) 


(w n T)* 


(- 


W T/2 


sin (w T/2) 


-) 


(6-6) 


? — 
* __ W 


c* 


ar 


sin (w T) 


2^ 


w T 

sin (w T) 


or 

(A 


1 

(w n T) 2 


/ w T/2 \ Z + ) n 

sin (w T/2) / w T 
n 


w T 

sin (w T) 


r 


i 

(w n T) 2 


w T/2~ \ 2 _ f n 

sin (wT/2) w n T 


w T 

sin (w T) 
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Clearly ct and y* depends on frequency w. 
(Taylor - development around w = 6) gives 


o 


1 + t » <». T > 

(w„ T)^ 


A low-frequency approximation 


Po 


1 - \ (w T) 
J n _ n ' 

(w n T) 2 


1 - 


o’ <T 0 


<- n T) 2 - 

( w n T ) 2 


2 


(6-7) 


These relations correspond to the choice of symmetric differences to approximate 
derivatives. Formula (6-5) is not valid for the two extreme data points. For 
the leftmost point we use 

X 1 = c< I^ y 3 + ^ t Y 2 + #L y i (6 ’ 8) 


and for the rightmost point 

X N = y N + f 3 R y N- 1 + (f R y N-2 


(6-9) 


The corresponding coefficients are given by 

2 

, 1 / in T / 9 \ 


ft L 

ar* 


(w T) 
n 


2 cos w T 

<»» T l 2 




T/2 


sin (w T/2) 


^ n 


w T 
n 


ft 


T/2 V * In 
sin (w T/2 )) w n T 


w T 


sin (w T) 


w T/2 
tan (w T/2) 


r 
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Hence 


/Z p = ^ replaced by - ^ ^ 

^ B = * L W “ h ^ n replaced by - ^ 


* K = l - p p -r - , 


The low-frequency approximations are: 

* Lo = L 

-2 + 4 f n ( Wn T) 


z 3 


Lo 


< w n T ) 

^Lo = 1 ‘ ^Lo " Z 3 Lo 


and 


/* 


Eo 


-2 -4 f (w T) 

J n n 

"2 


(w n T) 


y- . 1 + L<'"» T > 

0 Ro 


(W n T) 

^ Ro = 1 Z 3 Ro _ (Tr 0 


(6-10) 


(6-11) 
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6.2.2 Interpolation and Shifting of Sampling Instants 


References 12, 13 and 14 are very useful foj this section, 
output data to be N s equidistant samples 


>r this f 

(h • 


<yk = ^(k-1) T s ) , l s <k«N s 


Assume the original 


(6-12) 


While these original output samples (yk) are found at instants t = o, t = T g , . . . 
t = (N s -1) T s , the error-analysis (see next section 6.2,3) or the preference of 
the user may dictate a different origin and/or different spacing (see section 6. 2. 1). 


t = t Q , t = t Q + T , ...» t = t Q + (N-l ) T 

We may assume that the signals y(t) are bandlimited for practical purposes, 
with finite energy (if necessary, after subtracting the average value). This 
means that if the Fourier transform of y(t) is 

PO 

Y(w) = ! y(t)e _iwt dt 

— ^ 

then Y(w) = 0 for |w | > a minimum value which we can set equal to |w_ w n . 


Then y~ is a number, often larger than one. The energy is 


E y = 


r , i /><w n ? 

I y 2 (t) dt = / | y (w)f dw < 

- J oo Jy*'"n 


CO 


For the samples 


jykf , 


one has approximately (exactly, if T s is sufficiently 


small; T s < -II — and if y(t) is truly band -limited): 


2 u w. 


¥ w n 


E = T e 

y ^ s 


N 


k = l 


y k^ 


The bandwidth w n can, if necessary, be estimated by a discrete Fourier- 
a nalysis (preferably the Fast Fourier transform techniques) on { yk } . The Nyquist- 
Shannon sampling theorem says that, in order not to lose information concerning 
y(t), the sampling interval T s must satisfy 


T s 




rr 

K w n 


If this inequality is not satisfied, aliasing occurs and even an ideal (low pass) 
filter cannot reconstitute y(t) exactly from f yk ( . Since the reconstruction is 
not performed in real time, a non-causal interpolating operator can be used. 



Therefore, one uses here the interpoiating function (corresponding to an ideal 
low pass filter for w^wjjyj-p)* 


sin WjnT t rr 

- with f* w n s < w INT s < 2 — h w n 


rr 


One has the interpolating operator with the narrowest bandwidth if = p.w n , 

For any time t the interpolated value y of y(F) is taken as 

— — -_ s -*/ sin WTMT (F - (k-1) T s ) 

y (t ) = 5. y<fe-i) T _) — - s 


k=l 


_ir. ( t - (k-i) t s ) 

S 


The interpolated and shifted values y^ are now given by 

N * Sin W INT (t o + ( j’ 1)T ' (k_1) T s> 


V 


>\ 

y k 


k=l 


IT 


(t Q + O' 1 ) T - (k-1) T ) 


T can be larger, equal or smaller than T s . In some cases, to avoid the Gibbs 
phenomenon, the interpolation formula might be changed slightly to produce a 
Fejer sum. 

This interpolation formula can be used for a variety of values of t Q in an interval 
of length T to obtain shifted sample instants of the output (yk) and hence input 
(xk), for the same set of original (y^.). 

6.2.3 Error Analysis for Single Second Order Dynamic Systems 

The equation (6-4) models the rate gyros and accelerometer s (but not the pressure 
transducers). Looking at the discrete relationship (6-5) used for inversion, 
one observes three immediate sources of e rror. 

A. Quantization Error 


Disregarding the inaccuracies in (X , ^ there will be an error 

<£ q X in Xk due to quantization, roundoff, reading error and noise in { y k j • 

The following reasonable model is used to describe those errors in j yk \ . 
One assumes that the errors have a mean value £ and further that 
the deviations around this mean have an upper bound h^. 
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In other words, for each k : 

average of y k - (y k) true =. t y 

\ y k - (y k )true - £ y 1 ^ k y 

As a result, because of the interpolation formula established in section 
6. 2. 2, the | yk( values have errors with mean £ y = £ y and with an 
upper bound on deviations around that mean which is strictly speaking 
infinite. However, use of the statistical formulas derived in the 
Reference 12, combined with the fact that £ y = £ y already takes care 
of systematic bias, shows that hy may also be taken as a realistic "bound 11 
on deviations around the mean for ^ y^. ^ , except for infrequently occurring 
highly unfavorable error distributions. 

Then, the x^- samples have errors with mean 


-A 



because of condition (6-6), and deviation around the mean bounded by 


A 



where 

1+ l n < w n T > + I l ~ ^ w n T >l + 1 2 -< w n T > 
(w n T) 2 

Clearly c^ increases as (w n T) decreases to zero and c^ becomes 
infinite for w n T = o. 

(see also Figure 6-1). This consideration imposes a lower bound on 
(w n T). Since the purpose of this error analysis is only to establish a 
compromise value for T and no exact knowledge of the error is required, 
we will use here for the quantization error the following relation: 

^qx^l^yl +c q ky 

a conservative bound, reached only in rare cases with a very unfavorable 
distribution of noise errors. 


<=,- Woi + m»i + iir 0 
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B. 


Operator Error 


Besides the quantization error, there is still an operator error 
£ Dx * n c * ue to re pl acement of the operators (X , ^ 

by the low frequency constant multipliers Q , Q , Q * Again 
Reference 12 gives background information for the following developments. 

For the system modeled by equation (6-4), the true inverse operator 
is: 


G (w) = 1 + i 2 f n (w/w n ) - (w/w n ) 

The approximate inverse operator is: 

D (w) = ^ o ex P ( iwT ) + ft O + tf 0 ex P (-iwT) 


= 1-26. 


T) + i 2f si 


sin w T 


(w n T) z 


w n T 


One has 


€ < y 

DX * 


f w n 


G(w) - D(w) I d w 


*^ w n 


The operator error £ ^ increases when T increases (for fixed pt w n 
of course), unlike £ which decreases under these circumstances. 

Indeed, to find £ £)X one nee ds an expression for 


f f* w n 

|&(w) - D{w) 

-r w n 



G (u.w n ) - D(u-w n )j 


2 

d u 


= w n 1 ^ (P» ^ ) 


where p = w T 
L n 
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Then 


■F 


I^tp. \ )=2 f |[-| (l-cos pu. ) -u. 2 ] 2 + 4^{£H12_“-u.) 2 \ d 


8 sin (p-p) + 4 (p-p) cos (p.p) + sin (}^.P) cos (p.p) 

<fP) 5 2 / .5 , 3 _ 

+ 2 (}*pr sin (pi p) + -UiE2 - _|_ (pip) ] 

4 f 2 w.3 

+ r— |_ (|*P) ‘ 4 sin (p-P) + 4 (p-P) cos (p'-P) ~ sin (p-P) cos (p P) 


( pP) : 


+ | (p-p) 3 ] 


A series development shows that 

y p ’ t > = -rIV « r p)4 + - • • > + ?58«r-p) 4 + ---> 


63 
, E w 

now e 2 < y n 

DX 21 r 




1/2 


Figure 6-2 shows the increase of (p, p ))with p = w n T, for a fixed 

value of ( h- = 2). The bound £ is also very conservative, being 
reached only by the rarely occurring most unfavorable combination 
of signal and transfer function. 


Because of the linearity of the system model the worst total error on 
x(t) is given by the sum 


£ qx + £ Dx 


Due to the different behavior of and £ with varying T, £ x , will 

reach a minimum for a finite Optimum value of T, say T Q . 

The preceding discussion on how to find T Q is based on the idea that 
for this particular application a minimax formulation: 

find: min \ max error on x(t) j 

T any time t in the 

record 

is more appropriate than, say, a least integral squared error criterion. 
Indeed, one wants to discover what exactly happened in the vehicle in 
case of an accident, where large peak signals can be expected and one 
wants to find these peaks as accurately as possible, rather than obtain 
a good average fit (which would give too much importance to more trivial 
portions of the signal x). 
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Since both 8 and E are conservative, £ x will be a very con- 
servative es tkmate, reached only in rarely occurring most unfavorable 
combinations. Actual errors can be expected to be quite smaller in 
most examples (see examples in this report). However, regardless of 
the conservatism of £ x , the formulas for this bound are useful in 
establishing a best compromise T Q . 

It is possible to find less conservative error estimates by different 
means. These means include "Monte Carlo" runs to establish expected 
values for the quantization error, and approximate calculation of the 
spectrum y(w) of y(t) by means of fast Fourier transforms, followed by 
numerical calculation of the exact formula for 8 If the probabilistic 

structure of the errors is known, more refined statistical modeling 
becomes possible. 

Note that £ x bounds the error at input sampling instants on j x^. ^ . 

The inversion results are recorded as samples ^ x^ [ and no interpolation 
is performed on x. Note that if y(t) were truly band limited, then x(t) 
would also be truly band limited, with the same bandwidth. However, y(t) 
is only approximately bandlimited, and therefore the approximate 
assumption of finite bandwidth is better satisfied by y(t) than by x(t). 
Therefore, it is preferable to avoid making any assumptions about the 
spectrum of x(t). Values of x(t) at new instants of time are better obtained 
by interpolation and shifting on | yk | . 

C. Parameter Error 


The third kind of error is due to errors in the values taken for the system 
parameters w n> ^ n , and, more generally errors due to replacement of 
the true system equations by (6-4). One can get an idea of the importance 
of this error by making parametric runs with neighboring values for 
wn, f n. Note that one obvious reason for parameter error is the fact 
that we use average parameters for all instruments of the same class, 
say all (-3g, +3g) accelerometers, while tests have shown that there may 
be a noticeable variation (more than 10%) among individual members of 
that class. 

There is also a fourth type of error, less obvious than the preceding 
ones. 
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D. Inaccuracy Due to Truncation of Data (Edge Effects) 

If the true output samples -J yk [ contain non-zero elements for t<o 
or for t >(N s _i)T s> then the interpolation formula produces an error 
in the values of y. due to truncation. Because of the shape of the 
interpolating function, this truncation effect is largest at the edges of 
the record and very much suppressed toward the middle of the record. 
A simplified estimation of this truncation effect is obtained below. The 
influence of the first omitted sample 


y 0 at t = -T s (left) or 

y Ns+ i at t = N g T (right) is the most important. Assume that \ y Q | = 
and | yjsfs-i-i | = 1 • U se an avera g e for the absolute value of the 

interpolating function. 


y 1 


sin wjnT t 


Tf 


1. 5 re t 


av 


Then, at a distance ] A t | from the edges, one finds an error 

T 


( I y t l or 1 y N 1 ) ' 


1.5TC (T + |4 tD 


If Wjj^rp « Tr /T g several truncated values will contribute; however, 
there will also be more correlation of errors. For this reason the error 
can be estimated as 


(I yi I 


y Ns I ) T 


r T s 


51 r (T s + I a t l ) 


with 


<T - 1 if W INT T s > 


X ~ 3 if W INT^' S 


In order to have compatibility with the rest of the inversion, the truncation 
error should be of the order of h^. But, due to the worst case 
assumptions of the error analysis, one may allow truncation error 
estimates an order of magnitude larger than hy, say 10 h^. Then one 
finds 


* A t Heft or 
right 



1 
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as an estimate of the distance from both edges where results for ( x k | 
must be treated with caution. If no interpolation is used in the inversion, 
the truncation error does not arise. In many cases of truncated records, 
extrapolation by eye of | y k j before inversion can be used to decrease 
the truncation error. 


6. 2. 4 Inversion of a Cascade of Second Order Syst ems (Giannini Pressure 
Transducers) 

Now we take up the case of a cascade of two systems like (6-4), as exemplified 
by (6-3A). , or part of (6-3B). Such a system can be decomposed in two sub- 
systems, at least for the purely mathematical part of the reconstruction as 
shown in the normalized diagram below: 


x(t> 


Hence 


Z(s) 

X(s) 



1 

z(t) 

1 

1+2 > 
2 

s + s 2 

=> 

1+?. > s + s 2 

w 2 w 2 2 


1 W 1 W 2 1 


1 




y(t) 


1+2 f 


s . S L 


2 w 2 w 2 


Y(s) 

Z(s) 



W 1 


For each subsystem the formulas of Section 6. 2. 1 can be applied immediately. 
However, the error analysis is somewhat different. First, in order to have 
the intermediate signal z(t) as smooth as possible, one chooses the following 
order for inversion: 


First: 


wi = max (w n , w or ) or max (w n> w b ) 


h 

followed by: 


is the corresponding damping coefficient 
W 2 = min (w^ w Qr ) or min (w n , w^) 
is the corresponding damping coefficient. 
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Hence, W 1 > w 2 . If = w 2 , then ^ ^ is taken as the lower damping 

coefficient. 

We can assume a bandwidth of y(t) in the form 


pL w 2 

requiring a sampling interval T s for (y k ) given by 


T 


s 




Tr 

p.w 2 


Then the intermediate signal z(t) has approximately the same bandwidth w 2 
and the same interpolating function can be used for both z(t) and y(t). 

Assume that a sampling interval T 1 is used from y k to z k , and an interval 
T 2 from z k to x k . We now calculate the error bounds. 

A. Quantization Error 


One gets in the same manner as in Section 6. 2. 3 that 


A 


£ qx ^ I ^ y * + ^ C q^2 h y 


where 


(C q } l 


1 + f j (w ! T 1 ) + \l - f ! (wjT^l + \Z - (wjTQ 2 


(wi T 1)^ 


and (c ) 2 is given by a similar relation, where index 1 is replaced by 
index 2.. 


This is again a conservative bound. 

Operator Error 

Proceeding as in Section 6.2.3, one finds (see Reference 12) 
± «DX ■ ( K ”’ Z Cc 2 °1 - D 2 D ll Y e d " 

i - Y -'”2 


= L (G 2 -D 2 ) G 1 + (G^Dj) D 2 ] Y e iwt dw 

w 2 


then 


£ 4 

DX 


2TT 


I (G 2 -D 2 ) G x + (G 1 -D 1 ) D 2 \ dw 


T' w 2 


i 
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or 


£ .$ 

DX 


Ir-L.) 

\2 TV J 


in . 


2 1/2 


MG-DaJGjI dw) 


-|A W2 


K w 


( I \ (G 1 -D 1 ) D 2 I dw) 1 | 


- |t^W 0 


While the preceding expressions can be calculated exactly, the results 
are fairly complicated and more precise than is needed here. Therefor 
we replace them by a simpler approximate relation (with some 
compensation of errors) 

1 

\ G r D il 



e. 



- P ' W 2 


T w i 


In terms of the pre- calculated expression for I (p, (** ), one can rewrite 

M" \ 


the last inequality as 
E 


DX ^ 
approx 


2 rr 


ill 

\S 1/2 ^ 1/2 \ 

)(J r Vf 2 I 1 J w 2 T 2« 't 2>1 + C W 1 X r (w i T l, t l>i ) 


The minimum of £ x = £ qx + £ can be found by calculating both 
sources of error for a net of values of Ti and T 0 , The bound £. for 
the cascade will be even more conservative than for the single second 
order system. 

The spacing Tj for the first inversion can be obtained if necessary by 
interpolating (y^.) into (y^) with T g satisfying 

T s s <TT /( f ^ w 2) 

This gives (z^) at a spacing T]_. The distance T 2 for the second inversion 
can be obtained if necessary, by interpolating these (z^-) into a new set 
(z^) with distance T^ provided T^ ^ T g . 


\ 

\ 
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If T i > T g , another set of (z^.) is to be calculated, say by shifting 

the origin (and sampling times) of the new (z^) to half way points of 

the sampling times of the old j z^ j . This process can be repeated 

until one has a combined set of I \ values with spacing smaller than 

T . 1 

s 

For many applications one can expect larger inversion errors for the 
cascade than for a single second order system, due to the need for 
higher order derivatives in the inversion process. This observation 
implies that, unless the transducer models are exactly known, and there 
is very little noise, it is preferable by far to restrict oneself to low order 
approximating models and that the supposed gain in knowledge of the real 
physical device by means of a high order model may be illusory in 
practical applications. Therefore, if a physical system consists of a 
large number of cascaded devices, it may be preferable, at least for the 
purpose of inversion, to take only the dominant blocks, i, w. , the slowest 
ones, and, where possible, to lump several fast blocks into a single lower 
order approximating block. Such long chains do indeed occur in the 
Scout System, where the entire telemetry transmission follows the 
transducer. 


In fact, a general theory exists for such inversion problems, 
6, 2. 5 Inversion of a Numerator in the Direct Transfer Function 


(T -Tubing of Bourns pressure transducers) 
The numerator factor 


Y(s) 

V(s) 



in the transfer-function (6-3B) gives rise to a conventional second order dynamical 
system when inverted. Here V(t) is defined as an intermediate (mathematical) 
signal, such that inversion of V(s)/X(s) becomes the problem treated in Section 
6.2,4. The order of inversion, with Y/V inverted first, then followed by inversion 
of V/X is again chosen to find the intermediate signal v(t) as smooth as possible. 


Assume that y(t) has a bandwidth 
min < w n, w b> 
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Then v(t) will have approximately similar bandwidth. For the inverstion of 
V/X that bandwidth implies a spacing in v(t) - samples T g v satisfying 


T s,v^ 7 r/( jumin( w n > w b ) ) 


A similar spacing T s> y = T Sj v is needed for j y k j to avoid aliasing. 

If h a (t) is the impulse response connected with the transfer function of V(s)/Y(s) 
then one has for the forced part of v(t) 


v(t) 


,t'=t 


forced 


h a (t') y (t-f) dt' 


t 1 =o 


(A free solution with 2 arbitrary constants could be added, see further). 

If now one has a set of ^ yk | with spacing T 
yk = y (t Q + (k-l)T), t Q = o 

one can replace the integral for v(t) forced by the discrete approximation: 

V K = 1 Z. k y k+1 -n < H n- H n-2 ) ' \ Vx (* k - H^) 

c n=l 

where f*' 

H = h (V) dt' = H (jT) j> o 

J d J 

t 1 =o 


Hj = o if j 4 . o 

i. e. , H.(t) is the step response, and T is sufficiently small such that y(t+T)-y(t) 
is nearly zero over the record-length, (i. e., y(t) is nearly constant over a 
length T). The latter condition is satisfied if T. w. m i n (w n , W] 3 )«l (See 
Reference 12). As an example, one might take 


T ^ 

' K . min < w n> w b > 

while, for simplicity, T is also taken as a submultiple of T g . For H. one has 
further, explicitly 

e ^ a ^ cos ( V 1 - f w j T - Arc sin V ) 

tt = 1 _ \aa J ‘a 7 
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Once / | is found, the rest of the inversion proceeds as outlined in Section 

6.2.4. 1 1 


However, there is one particularity already mentioned before. Since only the 
forced part of v(t) was used, assuming v(t = o) = o and v(t = o) = o, one can add 
a free solution corresponding to 

JL + 2f V = o 

2 'a w 

w a 

a 

This means that as far as x(t) is concerned, one can add a free solution of the 
same differential equation with x replacing v. 

If somehow x(o) and x (o) are known, then 


x free 
x free 


(o) = x(o) - x forced (o) 
(o) = x(o) - x forced (o) 


The free solution to be added, at the sampling times of xk is: 

- f = „ I V 

x free (t) = e C x free (o > cos V 1 " * a V 

+ < fr ;; (C) + ^ a x fre e (o) > sin ( \A "Ha W a t}/ V 1 ‘ ^ a i 

In some cases x(o) and x (o) can be estimated. Say, that before t = o a steady 
regime is established, x(t) is a constant or a ramp for t <o 


x (t) = x(o) + x(o)t J t ^ o 


then y(t) is also a constant or a ramp with 
y(t) = y(o) + y(o) t j t ^ o 
where, for K = 1, 

k(o) = y(o) \ \ \ 

x(o) = y(o) + y(o) 2 ( _2_ + — J 

w n w b w a > 

Xforced^ 0 ) * s obtained directly from j x k J , if no shift is made in origin, while 
^forced(°) 2. X ^~ X ^ (T x |t w n « 1) or some other relationship, but this is subject 
T x 

to large errors ! 
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Other known boundary conditions can also occur, e. g. , at the tail end of the 
record to be inverted. In all cases two boundary conditions are needed to 
specify x(t) completely. Superposition can be used to eliminate the need for 
specific formulas for other boundary conditions. 

The error analysis for this case is not complicated. One needs an estimate of 
errors in v(t) based on the errors in y(t), in order to be able to use the results 
obtained for the simple case (Section 6.2.4) of the Giannini models. 

Using the formula for v^, one finds immediately for the bias in the error 

^ vk “ ^ y H k-1 

Hence 

1 f Vk I < l£y I I H k-1 I = i (1 + overs h°°t) 

y max y 


A good value, for k not very small (neglecting boundary effects) is then 


£ 


vk 


£ 

y 


As for the unbiased part of the error, one obtains the bound 


h 

v 




h 

y 


using the fact that the step response j h a dt was normalized. Finally, 
of the normalization of h a (t) again, one may take very approximately 


because 



for the signal energy. 

A 

Hence, one simply uses the error description of y also for v, without any change 
at all. 


6. 2. 6 Some Remarks on Inversion 

This section discusses a few problems arising in the inversion process. 
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A. 


Order of the Discrete Inversion Operator 


In Section 6.2. 1, a discrete second order operator (i. e. , second order in the 
basic delay operator e“^s) replaced a continuous second order operator (second 
order in the basic differential operator s). A higher order discrete operator 
could, in fact, approximate the ideal operator better. However, better inversion 
results generally would only be obtained if y(t) were known exactly, without any 
error or noise, or with very little noise. With noise present, a higher order 
discrete operator leads to a larger coefficient Cq (see the results for the cascade) 
and hence larger errors. 

To offset this, one would have to take a larger value of T and this would then 
negate the apparent gain produced by the higher order inversion operator. In 
fact, if signal and noise description are known, one can find the optimum order 
of the discrete inversion operator, but it would often be rather small in practical 
cases. 

There exists a more general theory covering this aspect. We refer to the end of 
Section 6. 2. 4 for a brief discussion of what to do with high order models, and the 
desirability to colLapse long chains of devices into a low order approximation, 
in cases where noise and parameter uncertainties are significant, as is usually 
the case. 


B. Interesting recent general results of a theoretical nature on inversion 
of linear systems can be found in the references listed under 15, 16 and 17. 

C. The present formulation was made with the models of the Scout 
instruments as specific applications. Therefore, the inversion of the fourth 
order Giannini model was reduced to successive application of the second order 
accelerometer model. However, one can also directly obtain inversions for 
higher order models, without the use of factoring. 

D. For very fast inputs (i. e. , fast compared with the response of the 
transducer), the reconstruction of the actual waveform is very dubious and is not 
even to be considered a well-posed physical problem. What counts is the total 
area under the fast pulse which acts more like a mathematical impulse. 

This effect is also demonstrated in Figure 6-3 for an accelerometer with 
w n “ 850 rad/ sec, = 0. 55 where all the inputs have equal area. Since output 
area is nearly proportional to input area, the strength of the impulse can be 
calculated, without any need for inversion by direct simulation and area- 
measurement with a planimeter (after the fact has been established that this 
is indeed a very fast input). 
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Figure 6-3, Uncertainty for Nearly Impulsive Inputs 
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E. In the present version of the digital program, the bandwidth of the 
transducer output signal is estimated by the user through the input y. . A 
poor choice of can lead to erroneous results (especially in interpolation). 

It would be useful, in an enlarged version of the inversion program, to have the 
computer estimate |y- by means of a spectral analysis. Such an analysis can be 
performed fast by the use of the fast Fourier transform. 

F. Edge-Effects 

The fourth type of error discussed in Section 6. 2. 3, due to truncation of y^. data 
can be (partially) removed in some ways. One method is to use such sampling of 
the original y(t)-record that interpolation becomes unnecessary. Another method 
consists of using a different interpolation scheme near the edges, such as spline 
functions. However, this approach requires further study to estimate best choice 
of order and type of interpolating function. A simple, but efficient way of 
alleviating edge-effects was mentioned at the end of Section 6. 2. 3; it can be used 
in the present program. 


6. 3 Org anization of_the Digital Inversion Program 

The theory of the inversion process used was given in Section 6.2, While the 
purely computational (numerical) aspects are very simple and straight-forward, 
the logic is a little more complex. Therefore, the present section describes 
the logical organization of the digital program for the reader who wishes to know 
its details. 


Inversion Program 


Inputs: 


L so 

N s 
T _ 


(y kf , 14 M N g 


gain K 

-A >\ 

Errors £ , h_ r 

y y 


Bandwidth 




optional 
option 

Choice W INT (IWINT) 


to 1 

t max l optional 


Option to stop after error calculation (ISTOP) 

Type of instruments: Class CA: Gyro, accelerometer 

Class CB: (Giannini) pressure transducer 

without T-tube 

Class CC: (Bourns) pressure transducer 

with T-tube 
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Special Inputs 

If CA: w n , ^ n ; jT \ ; jT ch \ 

If CB: w n> f n ; w or , f orj {T 1 } . (t 2 \;\? 1 ch> T £ ch 

If CC: wn, ( n ; w b> \ b; w a> \ a ; |Tj_ | , ( T 2 } = 

i T l ch' T 2 ch 1 
oL (optional) 

y Q (optional) 

y (optional) 

o 

Notes: 

1. { T ch} : ( T i ch^cif” °P tional (NOPT1) 

If they are used, they should have values covered by 

J, (T 2 ^ if one wishes to have an error calculation. 

2. One should have 

t £ t / t 
so * o * max 


or 
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Main Program 


Calculate E y 

If t Q not listed: t Q = t go 

If t not listed: t = t + (N c i )T 

max max so ' s-i' s 

Check if Class CA, CB, CC, then go to appropriate subprogram (CA), (CB), 
(CC). 

* Calculate call results J x^ j 

Print A t-^, A tp (edge effects) if applicable 
Print t Q , T, N 
Print / Plot j t^, x^ ' 

If Class CC, moreover, (for formulas see section 6.2.5) 

Head y Q , y Q (set zero, if not given) 

Calculate x(o), x(o) 

Calculate x free> x = x free (o) , x forced (o) = ^ 

*free, 1 = *free (°) * *forced(°) = Xz ~ Xl -. 

T 

x 


Calculate ( x f 1 

j (t D + (k-1) T)j 

call results j x free> fc | 
Print |t k , x free> k J 
Print/Plot I t k , x free> k + x k j 


1 4 k 4 N 


Stop 
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Subprogram (CA) 

For j j ||U. : Calculate I (p f ^ n) 


p = w n T 
£ DX 


Calculate c , £ , £ 


qx 

Print J_jc_ T „ ( call results £ ) 

K r x 

Check ISTOP (stop 

or 

(continue 
If continue, then 


Choose for each jw. , unless overridden 

by | — ► | T ^ is ultimate choice 


- If T 1 = T and t = t 

S o so 

Then [t 0 , N, T, t k , x k ] = BASINV [ t Q , N g , T g , £k, w n> } n ] 

(see subroutine BASINV further on) 
and return to MAIN at point* 

Otherwise [ t Q , T, N, t k> yk] = INTER P [ t gQ , N g , T s , yk, y*. , w Q , T', 
V t max J W INT 3 

(see subroutine INTER P further on) 

tV N » T » V *k^ = BASINV t V N, T, yk, w n , t n 1 


and return to MAIN at point* 
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Subprogram (CB) 


Calcuiate wj, W2 and ^ \ % 

according to section 6. 2. 4. 
For ^ , ( 1*2 | 5 j K* | calculate 


Pi , 


1 I (P 2 . ^ 2 > 


DX 


Calculate (c K, (c ) 2> £ 


qx 


Calculate £ 


Print — f T 7 , (call results £ ) 

K. ^ I x 


Check ISTOP (stop 

or 

(continue 


If continue, then 


Choose t , T ^ Q p t for each p. , unless overridden by 

^T lch , ch | |t* 1 » T '2 | is ultimate choice 

Now, for each p , enter (CA) at point 
with f T' = T'j, w n = W lf 




W2 

Wi 


pc > other symbols are correct 

input 


and leave at 


with 


>\ 


N 1 T 1 = T » z k ~ x k 

1 other symbols unchanged 
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Then, check whether Tj^7r/(^W2) 
If yes, go to P 2 
If not, determine whether 
If yes, go to P 2 


>‘P 


If not, find smallest integer 1, such that 2 ^ 7T /(ptW 2 ) 

and theja for j = 1, 2, . . . , 2^ - 1 go to E^ with t Q replaced 
by t 0 +j Tj2“P and otherwise same conditions as first entry at Ej 
above. 


Leave again each time at L^ with z k, j = x k* j = ^k an< ^ ot herwise 
conditions as above. 

Now, reorder all j z^. j j in order of increasing time 
^ t o > 1 2 <• • • • 1 0 - • 


A 

Zi 


z- 


1 "1, 1 Z l,2 Z 2 "2, 1 “ 


Call new series 


i Z k’ *k \ 


(Still for each ^ ) 


Enter (CA) at point E^ 


with 


' o, W = W9 

^ n c 


’ In " I 2 
P 


so 


t G , N s = N X or Nj 2 


Values for w TT v T rp, t , t , u. are correct 
INT* max o r 


1 


< y k 




T s - T : or T l 2-P 


Leave at point Lj with N, T 


Return to MAIN at point * 
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Subprogram (CC) 

T 

Take T = — — , with m smallest integer (m>, 1) for which 
v m ' 

^s ^ <*L 

m v max ( p. ) • min (w n , w^) 

(if <p{ is not input, use cK. = 0 . 2 ) 

Print T 

v 

Enter (CB) at point E^ , leave at point 

Check if T = T . 

v s 

If yes, proceed to Pj 
If not, then for each pc 

[t so , T g , N g , t k , yk"] = INTER P N g( T g) yk, IL 

w l’ T v ’ ‘so’ max’ W INT 1 
P^ Calculate ^ , see section 6.2. 5 

Set new | yk j = j vk j , and enter (CB) at E^ 

Leave at L 3 

Print "forced solution" 


Return to MAIN at point * 



Subroutine for Basic Inversion - ,r BASINV M 


Inputs t , N, T 

o 

, {y k l ’ w n’ ^ n 
Calculate tk 

* o 

* Lo I^Lo X Lo 

* Ro /’Ro^ Ro 

(k^l, k^N) 

x l» X N 

Outputs f t Q> N, T 

I f V x k ] 

Symbolically: 

[ t Q . N, T, t k , = BASINY [t Q , N, T, yk, w n , 1 
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Subroutine for Interpolation - "INTER P M 


Inputs 

‘so' N s> T s 

|v k ( 


H-* w n 


‘o’ T ’ ‘max 


W INT ; °P t i° na ^ 

Verify: 

T s ^ ft /(yw n ) 


If no: warning "sampling interval too large" 

If yes: proceed 

W INT = ^ w n » unless input (IWINT) commands 
W INT = ^ / T s (wideband) 

N largest integer such that 


t n + (N-l)T v < (t max ; t so + (N.-l) T_) 


t k = t s o + ( k -D T s belongs to yk 
tj = t Q + (j _1 )T belongs to 


Interpolate 

N 


v s * 
Yj = yk 

J k=l 


sin w INT (tj - t k ) 


'TT - ' , ^ 

^ <V*k> 


1 C j<N 


Outputs 


t 0 . T > N 

i t k* y k ( 


Symbolically: 

][ ^o* r ^ , » 1 = 

INTEBP [ t ao , N s , T g , ^k, y , T, t Q> t max , 


INT 
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DIGITAL COMPUTER INPUT I 
REQUEST FORM 


PROBLEM NO. MEMO NO. 


CONTINUATION SHEET 


DESCRIPTION OF INPUT 


Input Common to Classes CA , CB , CC 


Number of points of Y (500 max. ) 
Time spacing of points of Y 


*NMU 


Number of values of (10 max. ) 


*ICLASS 


Values of 1, 2, 3 choose CA, CB, CC respectively 


*N0PT1 


Value 1 allows program to choose its own time spacing; 
value 2 imposes choice of spacing by user 


*NPL0T 


Value 1 allows for plotting; value 2, no plotting 


*IST0P 


Value 0 bypasses determination of (KX); value 1 
calculates (KX) . Values 0 and 1 both allow for error 
calculations 


Initial value of time for Y 


Initial value of time for (KX) (preset to TSO) 


KG A IN 


Value of gain K of instrument 


TMAX 


Value of maximum time for (KX) (preset by program) 


ALPHA 


Value of (preset to 0. 2) 


*IWINT 


Value 1 allows for narrow band , value 2 allows for 
wide band- 


Value of ^ 


A 

Value of h 


*NTCH 


Number of chosen values of spacing imposed by user 
(10 max. ) 


YBK 


Values of Y (500 max. ) 



DIGITAL COMPUTER INPUT 
REQUEST FORM 


PROBLEM NO. MEMO NO. 

2894 


CONTINUATION SHEET 
PAGE 4 OF 5 PAGES 


WN 

ZETAN 

MU 


Value of O 


Value of 


Values of ^ (10 max. ) 


Input for Class CA Only 

TCH Values of chosen time spacing imposed by user (10 max. ) 

*NT Number of values of T (10 max.) 

T Values of time spacing , T (10 max. ) 

Input Common to Classes CB , CC s 


T1 

Values of T (10 max. ) 


T2 

Values of (10 max.) 


*NT1 

Number of values of (10 

max. ) 

*NT2 

Number of values of (10 

max. ) 

T1CH 

Values of T lch imposed by 

user (10 max. ) 

T2CH 

Values of imposed by 

user (10 max. ) 

Input for 

Class CB Only 


W0R 

Value of CJ 

or 


ZETA0R 

Value of P 

S or 


Input for 

Class CC Only 


WA 

Value of 

a 


ZETAA 

Value of 



9 




2-0064 

>-63 


PROBLEM NO. 


I MEMO NO. 


[section no. 


DIGITAL COMPUTER INPUT 
REQUEST FORM 


2894 


CONTINUATION SHEET 
PAGE 5 OF 5 PAGES 


WB 

Value of CO, 
b 


ZETAB 

Value of 


YO 

Value of Y 

o 

(preset to 0 # 0) 

YOD 

Value of Y 

o 

(preset to 0. 0) 

* Must be 

integer 



Warnings 

(1) TSO< TO TMAX 

(2) Program stops if number of interpolated points exceeds 500 and prints 

"Number of points are too great" 


(3) MU, WN, W0R , W A, and WB must be positive and have exponents 
between -3 and +10 in order for the plotting routines to work. 



6.4 


I llustrative Examples 


First Example: 

Take the accelerometer with w n = 850 rad/ sec = 0. 85 rad/msec, = 0. 55; 
triangular input. The computer plotting routine plots the transducer output 
(y) and the reconstructed input (x). The plotter also connects x- samples by 
straight lines, to facilitate reading. The true input has been added (by hand) 
to provide an immediate criterion of the quality of the reconstruction. All 
necessary parameter values are given in the computer printout: among them, 

bandwidth fX = 2, error h^ = 0. 025, sampling time T g = 0. 5 ms. 

Figures 6-4 and 6-5 give the reconstruction for smallest total error <f x , 
corresponding with T = 0. 6 msec; for interpolation with w^-p = /a. w n = 1, 700 
rad/ sec (narrowband) and Wjjsjt = 'K = 6,280 rad/sec (wideband) respectively. 

T 

s 

Figure 6-6 shows the reconstruction for T = T s avoiding interpolation. All three 
results are quite acceptable, indicating that in this case the bandwidth was 
properly estimated. 

Second Example: 

This is a Giannini pressure transducer with = 1, 300 rad/ sec, f n = 0. 2, 

w or = 300 rad/sec, ^ Qr = 0. t ^ Le bandwidth factor is taken as either 1 

or 2. Results are shown in Figures 6-7 and 6-8 for T^ and T^ selected for minimum 

total error. Clearly, = 1 is a better choice here; the edge effect is very 

noticeable at the right end. Figure 6-9 pictures the reconstruction for 

Tj = T 2 = T g = 0.625 msec, avoiding interpolation: obviously no edge effect appears. 

Third Example: 

This deals with a Bourns pressure transducer, together with its T -tubing. 

One has w n = 400 rad/sec, w a = 1, 000 rad/sec, w, = 200 rad/sec, 

I n =0 - 3 > fa = 0 ' 8 ’ ^b = 0 - 8 

A 

T^ = 1 msec, h = 0. 015, [-<-* =1 and 2, y^ = 1. 0 (initial value) 

The results for the total solution (forced plus free) x are shown in Figures 6-10 
and 6-11 for K- = 1 and 2 respectively, with narrowband interpolation (wjj s y-p = pc w^). 
There is a very strong edge effect. Figure 6-12 gives a much better reconstruction 
with wj nt = Tr' /T g (wide band interpolation) indicating that ja should be 
chosen larger here than 2. 0. 


A 

m 
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NOTE: 


All examples show that the calculated error bound is indeed quite 
conservative, as predicted. However, the choice of spacing 
intervals found by minimizing this worst error is quite good and 
this is the main purpose of the error calculation. It would be 
possible to find smaller error estimates by use of statistical 
methods. However, to maintain generality, a non- statistical 
error calculation was preferred here. This made possible avoidance 
of additional hypotheses about signal and error distribution. 

The computer printout corresponding to the plots follows after 
Figure 6-12. 
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INPUT OUTPUT 
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INPUT OUTPUT 




INPUT OUTPUT 








INPUT OUTPUT 
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INPUT OUTPUT 
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INPUT OUTPUT 
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INPUT 

CARES READ 

14 .44 * 2 2 • 0 c 

DATA* 

DATA* 

NS 15 TS 0.50C NMU 1 ICLASS 1 N0PT1 
NTCh 1 -IWI-NT 1 KGAIN l.CC<| FB Y ■ ft-.-flO'T 

1 I STOP 0 TSO 3.0 CO TCH 0.500"" 
may n.^^nn vr^ n.Ann /-j.rpRrr ... 

DATA* 
DAT A* 

0.15DG 0.4500 0.9250C 1.475DC 1.90DO 
l.onn 0.55DP 0.2GDO -0.05DC HU -2-»0D?i 

2.1500 2.1400 1.8700 1.4500 
N T 4-0 T o . 1 ns'- { 0 . 1 no ) 1 . r> CO 

DATA* 

nAT AtSc 1 

WN U. 85000 ZFTAN 0.5500 NPLOT 1 




T.D i IF? TN THP TNTFRVA I ( TSO ,TPAX) 


MU 

- 

d. ?Q MnC 

_01 _ 

J 

— 

r .i Muon 

nn 

_ El_ 

= 

0. 1 38 29(1 

1L? 

MU 

T 

0 . 20CCOD 

01 

T 

= 

C. 230000 

DC 

ex 


0. 34862D 

0 1 

..M.U 

= 

o. ?nocnn 

01 

j 

= 

t. . 30000 n 

on 

FX 

= 

O . 1 A ? 7 1 h 

in 

MU 

= 

0.200000 

01 

T 


0.400000 

t?C! 

EX 

= 

U. 1041 RD 

01 

MU 

= 

0. 20 CO OR 

_OJL 

j_ 

— 

o . 5 fl fi ' > c n 

0 3 

ex 

= 

3. fi41 r AD 

:*> ij 

MU 


0. 200 000 

01 

T 

- 

C. 600000 

00 

EX 

= 

0.804470 

■■VO 

MIL 

- 

0.200000 

01 

j 

— 

r.. mnr en 

ms 

_ FX 

— 

U. 85 544D 

<1/? 

MU 

= 

0.200 con 

0 1 

T 

* 

C. 8 000 CD 

no 

EX 

= 

0 • 9 6 0 fi 3 D 

j>0 

.. MLL. 

— 

■ ■Q.20QC0D 

.ill 

I_ 

- 

C . 90 OB 00 

cc 

EJL_ 

= 

0.11 0 37 D 

0 1 

MU 

= 

0. 20CG0C 

01 

T 

= 

c . mo one 

<n 

EX 

= 

C. 12741 D 

01 


Tfl * O.f: TF = r.60000D 00 N =12 PU = 0.200000 til 


DTLE.fi T =-ft « 5CCCC D 03 DTRIGHT = -u .47 8780 00 


TSUB ( 1) = 0.0 X = 0. 62 67 3D 00 

JL.-6L n C -J-D S P X = _S . 11.01 6C_t?3 

T SUB ( 3) = 0. 120000 €1 X = 0.18856D 01 

T.SUB< .4) = .0L.18MCQL.G1 X_*_2*23J6iD.i>JL... 

TSUR ( 5) = 0. 240000 Cl X = 3.233480 Cl 

TSUBl 6) = n. ci _x..= 0 . 1 7743 c 01 

TSUB ( 7) = 0 . 3 6tO(;D Cl X = 0.96412D ‘'0 

. .XS.U.BJ . JLL.= n.42tf"D Cl X = 0. 26 J32 D - .M 

TSUB ( 9) = 3 .4 80 COD Cl X =-'1.655160-01 

_T S U 3 { 1 0 ) = O.540CCU) £ 1 . X =-C. 246310-11 

TSUB(ll) = O.fcOCrCD Cl X = C.19974C Gfl 

TSUB (12) .= C. 6.6CC.CP_C1. X = 0.20 30 70 ^ 


FOR F i <£r 6 " 4 
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r 


INPUT CARCS READ - - - - . 14.44. 26.01 

DATA* IWTNT 2 

DATA*! .... . ... ... 


TO L IES-.LN— THE ^INTERVAL (XSO , TMAXT 


MU = 0.20CC0D 01 
MU = 0.20CC0D 01 


MU = Q.20QC0C 01 
MU = C . 200 COD 01 


MU. =. 0. 


MU = 0.200000 01 
MU = 0.200CC.D 01 _ 

mu = o. 20 ocnn 01 
MU = . C..22OC0C 01 
MU = C.20rC(!C Cl 


T = C.10000D 00 EX = . 0 .138 29_D_._0.2 

T = 0.2 0000 D 00 EX = 0.34862D 01 

T = C.30090D DC _.E_X = 0.16271.0, 01 

T = 0.400000 00 EX = Q.10418D 01 

T = C. 50000D 00 EX = 0.841040 00 

T = C.6O00CD 00 EX = 0.80447D on 

T.t 0.7G030D 00 - __£X_= 0^85544D DO 

T = C.8DO00D 02 EX = C. 960830 00 

I = C. 900000 00 .. _ EX .= C, 110370 01 

T = C. 100000 01 EX = '5.127410 01 


TO = 0.0 


TF = 0.600000 CO 


N =12 


MU = 0. 20000 C 01 


DTLEFT =-Q.5CC0CD 0Q DTRIGHJ.. =-0.4.78780. Q.Q. 


TSUBI 1) = 0.0 X = 0.5690 1C 00 

TSUBC 2) = g.fcOCOOD (!0_ . . X__r C.98741C. DO 

TSUBI 3) = C.120COD Cl X = 0.19343C 01 

TSUBI 4) = 0.180000 51 . X..= . 0. 26 8'j7D.. 'll 

TSUBI 5) = 0.24COOD Cl X = 0.21464D 01 

I SUB .! 61 . =...J3 ...3.C.0CDD. £1 X_«.X>.a7655 D 01 

TSUBI 7) = 0. 360000 01 X = S.90034D DO 

TSUBI 8) . = 0.42CCQQ Cl __X_=_ 0._18293D 0_D 
TSUBI 9) = 0.480000 Cl X = 0.25484D 00 

TSUBI 10) = 0 . 54 O' CD Cl X =r 0.135490 00 

TSUBI 11) = 6<-'CCCD Cl X = 0.4306 7C-0 2 

T SUB I 1 2J =..n, 6 6f .CCD . Ql_. _. X =-0.136650 OP 


Ir QR FI 6 . <3-Z-5. 
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INPUT- -CAR-OS— READ 

DAT A* NOPT1 2 

DATA*! 


TO — LJ-i-S— IN J XE INTERVAL (TSC . TNAX) 


MU 

— 

r.pnntn r cn 

J 

= G.iaaacD or 

FX 

- 

0 . 1 38290 

:i ? 


MU 


0.2GCC0C Cl 

T 

= <G. 2 00100 00 

EX 

= 

C. 348620 

oi 


MU 

= 

O.ZOC^DC _.Q1 ... 

J_ 

= o.3 none p or . 

FX 

— 

£L, 162710 

_31 


MU 

= 

0.230000 01 

T 

= 0.40C30D 00 

EX 


0.10418D 

01 


MU 

- 

o . por non 01 

T 

= 0j_500'?ilD 1C 

FX 

— 

0.841 £40 

no 


MU 

= 

' i j • 2 u C FO ? C 

T 

= 0 .600 CCD 00 

EX 

= 

(,.6’>447D 



-J4U 

- 

€*2511.001) 01 

_T_ 

= r.700/ ion on 

FX 

- 

_J1*8554AJD_ 



MU 

= 

0 • 200 CO D 01 

T 

= C.8C0C0D OC 

EX 

= 

7.96C83D 

00 


„ MU- 

- 

0.20CC0D ,ijl_ 

T 

= c.9woo no 

FX 


... 5,110.3 7B_ 

.'3 1 


MU 

= 

0.20CC0D 01 

T 

= c.inoooD oi 

EX 

= " 

"'• 127410 

:?i 



TC = G.C 


TF = C.50*)CPD OG N =15 PU = 0.2CC00C 01 


.ISUM._1) = JL.L 

TSUB ( 2) = 0.5PC0CD OP 
T S U B ( 3) = 0. ICC COD 01 
TSUB C 4) = C.15CC0D 01 
_TS UB i ... 5 L = C.2CCC0 D 01 
TSUBI 6) = ('.25PCOD Cl 

t$.ub.( _7_> .=..a,.=nfi.c.an_:Ci 

TSUB ( 8) = 0.3 5GCGD 01 
TsiiBt 9) = c.4pn<:en c i 
TSUB (1C ) = C .4 50 COD Cl 
TSUB Ml) = 0.50CC CD Cl 
TSUB( 12) = C.55COOD C-1 
..JSUail3j_= C.feC flOCD Cl 
T S U B ( 1 4 ) = G .6 5CC0D Cl 
TSUB ( 15 ) = C.70CrCD 01 


X 

X 

X 

X 

X 

X 

x 

x 

x 

X 

X 

X 

x 

X 

X 


= Q.77275C Cl 
= 0.1668SC 01 
= 0.242180 Cl 
= 0.266670 01 
= 0.2044 7D 01 
= 0.18047D Cl 
= 0.1021 1C 01 
= 0.338200 00 
■= 0.1466 ID 00 
= 0.158030 CO 
=-0. 16471C 01 
= 0 . 6833 SD-01 
=-C . 228370-1 1 
=-0.140140-01 


FOR FIG. 6-6 



input 

CARDS READ 

14.44,27. n_ 9 _ 

3/03/70 

DAT A* 
n aja* 

NS 17 TS 0 .62500 NMU 2 ICLASS 2 

7PTA0R n.?nr5__A!-M ? NT? ? NTTH 1 

NCPTl 1 I STOP 0 WOR 1.3D0 KGAIN 1.000 
44CH 0 .62^0'“ -72 rw C,^2^nn iwlNT 1 

♦ DATA 

*pata 

DATA* 

DATA* 

TSC O.ODf- EBY 0.000 HBY 0.62500 
AS.nr^ in7,Bnn ] a*, 171.7800 

YBK O.CCO 0.( 00 1 * 8800 11.2500 3C.0C0 
r , 7 son iQi.78nrs iftft.orr 

♦ DATA 

♦ HAT A 

DATA* 
OAT A* 

186.000 187. 5Df 193.75DC 2CC.CD0 
t i o.^nio fi.Apcno r.7nrt t? n.sno 

V|j 1.009 2.0 DO WN 1.300 Z ET AN C-.2C0 
fl-A2 5n(! 0.700 _ .. .... 

♦ DATA 

♦ DAT A 

DAT A* 1 



♦ DATA 


TO LIES IN THE INTERVAL (T$C,TMAX) 


MU 

- 

0. 100000 

01 

T 1 

= 

C. 500020 

00 

T2 

= 

0.50*Q0D 

00 

EX 

= 

0 • 5 5 5 5 8 D 

02 

MU 

— 

0. 1G0C0R 

01 

T 1 

- 

C. 5C0QCD 

CO 

T 2 

- 

. C. 625000 

ns 

FX 

- 

0. 4*3 80 ID 

02 

MU 

= 

0. 100C0D 

01 

T 1 


0. 5000QC 

00 

T 2 


0.700000 

on 

FX 


0 • 365840 

02 

MU 

= 

C« lQG'uOD 

Cl 

T 1 

= 

C. 62 500 D 

no 

T 2 

- 

n. 500000 

..00 _ 

EX 

= 

0.A0150D 

0 2 

MU 

= 

0.10000D 

ni 

T 1 

= 

C.6250CD 

CO 

T 2 

= 

*.625*00 

00 

EX 


0 • 32 654D 

02 

MU 

= 

o.ioooon 

01 

_ 11 

- 

wmmzzm 

CO 

T2_ 

= 

ravzsznnaa 

EuEi 

EX 


WVUt'Ud.1 

_£2_ 

MU 

= 

0.100C0D 

01 

■9 




B 

= 

0.500*00 

00 

EX 

= 

0.354810 

02 

MU 

= 

0.10000D 

01 

■n 


Rn SI 


mn 

- 

C. 625000 

E3v-1 

EX 

— 

BRUM* 

02 

MU 

= 

0. lOCCOD 

?1 


= 

C. 70CCC0 

00 

T2 

- 

0.70C00D 

00 

EX 

= 

0.299470 

02 


MU 

= 

C, 200000 

01 

B 


C. 50000 C 

00 

T 2 

= 

C.5000GD 

00 

EX 

= 

0.22285D 

03 


- 


01 

sn 

- 

ItHAfeUihl 

CO 

T 2 

= 



EX 

— 


03 

MU 

= 

0. 2QOOQD 

Cl 

T 1 

= 

C. 5000 CD 

on 

T 2 

= 

0, 7OC00D 

00 

EX 

= 

0 * 27855D 

C3 

MU 

- 

0. 200000 

Cl 

T 1 

- 

C. 6250*0 

on 

_ T 2 

= 

0.50 or 00 

no 

EX 

- 


03 

MU 

= 

0.20CC0D 

01 

T 1 

= 

'it* 6 2 5 u v D 

oc 

T 2 

= 

0.625000 

no 

EX 

= 

0.28641D 

03 

MU 

= 

0. 20 GOOD 

•01 

T 1 

= 

0 • 62500D 

00 

T2 

- 

0.70000D 

00 

EX 

= 

0. 314910 

03 

MU 

= 

0. 20000 D 

■01 

T 1 

= 

Q • 700 CCD 

CO 

T2 


o. socroc 

no 

EX 

= 

0.272150 

0 3 

MU 

= 

Q. 20 GOOD 

01 

T 1 


7 00 QC n 

on 

T2 

— 

3. 62 BO-DO 

on 

FX 

- 

0,312240 


MU 

= 

0* 20000 D 

01 

T 1 

= 

C. 700*00 

no 

T2 

= 

O.7C0«00 

00 

EX 

- 

0.341680 

03 


• 

II 



TFl 

= V • 

700000 ac 

TF2 = 

O'.7C0Jt:D 

n-j 

N 

= 15 

MU = C.1O0OOD 01 

T $U3 ( l) 


c.u 


X 

= 7.196280 

0 2 





T SUB < 2) 

= 

r.70 fOCD 

r r 

X 

= 0.319800 

C 2 





T S U R ( 3) 

= 

C. 14COOD 

Cl 

X 

= C.472C3C 

n 2 





T SUB ( 4) 

= 

0 *210000 

a 

X 

= C.72317D 

02 





TSUB ( 5) 

- 

<3.280000 

a 

X 

= 0.935990 

*2 





T S U B ( 6) 

= 

3. 350CCD 

Cl 

X 

= 0.124310 

*7 3 





T S U B ( 7) 

- 

0 . 4 2000 D 

Cl 

X 

= 0.148120 

0 3 





T SUB ( 8) 

- 

; ^socod 

c?. 

X 

= C.1685CC 

C‘ 3 , 






„TS_UBJL_5J_ 

= 

0 . 5 6CGCD 

Cl 

X 

- C. 183180 

03 ! 

[for 

F 1 CS 

Cd 

- 7 


TSUR(IO) 

= 

(..63CC0D 

n 

X 

= 0.189350 

0 3 1 






_T_suR.a.u_ 

= 

7* 77CCCD 

Cl 

X 

= 0.184560 

03 





TSUB( 12 ) 

= 

C .7 7CC no 

Cl 

X 

= 0 . 16 8070 

n 3 





TSUILUL 

= 

C.64CCCO 

Cl 

X 

= 0. 141670 

03 





TSUB < 14 ) 

= 

C.SKCCD 

Cl 

X 

- C.94873C 

02 





T SUB ( 1 5 ) 

- 

r •90TCOD 

ni 

X 

= 0.351700 

02 






■t c .. - t uo. 


TFl - 0. 530,SCD. ac TF2 = O.5t)U3CU.Q0 N =20 MU = q.2nOQf;Q m 


TSUB ( 1) = C.n X = 0 . 203440 02 

rsue( 2 ) -=-Q..scccnri cn x— --0.26413D. 0 2 

T SUB ( 3) = O.r'SOOD 21 X = '1.125260 02 

xs.uai—^3— s_c-.X 5 f-fir.ti. u x_=-iL. 5i425n 02 

TSUB ( 6) = f.20CO?D Cl X = 9.945150 02 

XS.UB.L. 6 ) = ■-■.? 51-1:00 Cl X = 0 . 1C 45 BO 0 3 

TSUB ( 7) = 0.200000 Cl X = C.8995CC (22 


TSUB ( 8_)_? -0. 2 50000 01 X = 0.973520 02 

TSUB ( 9) = C.4OC0CD Cl X = 0.146060 03 

TSUB HO) = 0.4 5n;'no ci X = 0.1 935713 C 3 


TSUB (11) = C .5 KCCD Cl 


X = 0.188220 

0, 


for r i <$ , 3 - S 


T S U B ( 1 3 ) 

= C .6-0 (.'COD 

ri 

X 

= 

C. 133710 

1 ---=3 

T.S_U-6Q4.) 

= <.j65vLTOJD 

f * 

x_ 

- 

0.1 9 C 0 5 n 

3 7 

TSUB( 15) 

= c.7r?rco 

Cl 

x 

= 

0*254130 

03 

J.S.UB-UL6J 

= 0.7 5 C T 0 D 

01 

X 

— 

<1. 23759n 

0 3 

TSUBI 17) 

= r.80 r ^00 

Cl 

X 

tr 

".136800 

H3 

.TSUB (18). 

= o.85rcon 

01 

X 

— 

3.481470 

0 2 

TSUB ( 19 ) 

= C* 90 C CCD 

a 

X 


0*215330 

0 3 

T.S.U812DL 

s r. A c g f ( p 

Cl 

x_ 

= 

0.234260 







INP.uT--CJUU^S£AO 

DAT A* NOPTI 2 
DATA*.l 


1 ^ . 44 . 28 . 69 


3 / 0 3/70 
*OATA 
tOATJL 


T r LIES I N TH E— IN TERVAL- C TS Cf-T M AX) 


MU 

= 

c.iaorso :n 

T 1 

- 

C. 500CiiLC 

CO 

T? 

- 

n • 500 HOD 

00 

FX 

= 

0 . 5 55 5 PD 

02 

MU 

= 

0 • iDCD^r o\ 

T 1 

= 

r. 5‘jTOro 

A o 

T2 

= 

rt .625K‘C 

no 

EX 

a 

0. 408010 

02 

MU 

- 

o. loorno "i 

T 1 

= 

■: . s^-'oco 

■oo 

T 2 

= 

_*.U 700000 

00 

EX 

= 

0.365840 

02 

MU 

= 

r.lOCM-D M 

T 1 

= 

r .625?0C 

ro 

T2 

= 

D. 50000 D 

oo 

EX 

= 

0.40150D 

02 

MJJ_ 

= 

o.iooooc :i 

Tl 

- 

i! • 6 ? 5 0 C D 


T2 

= 

0.6 2 500 D 

no 

EX. 

— 

0 . 3265 AD 

02 

MU 

= 

0* IOC r . DC “il 

T 1 


C.625CQD 

r* n 

T 2 

= 

'j. 7 00 COD 

0 c 

EX 


0.31055D 

02 

MU. 

= 

"J. 

Tl 


•■'.700000 

re 

T? 

- 

n.soonnn 

00 

FX 

- 

0.354810 

OP 

MU 

= 

nr coo ~ 7. 

Tl 


C. 7S00CC 

CD 

T2 

= 

- • 62 5 r jD 

00 

EX 

= 

0.30603D 

02 

. MU 

= 

locrr.o ci 

LJL 

= 

r. 700 OGO 

00 

12 

= 

C.7000PD 

no 

EX 


0.299470 

02 


MU 


C.200CUC 

^1 

Tl 


C. 5030CD 

nr 

T? 


0.500000 

00 

EX _ 


0.22285D 

03 

MU 

= 

::.2'-c :cc 

-;;i 

Tl 

a 

si^roo 

nr 

s# V 

T2 

= 

‘ .62 5CDD 

PV 

EX 


0.253C6D 

C 3 

MU 

— 

l? • ? D C V 0 D 

"1 

Tl 

= 

s*"ccn 

A **. 

_ T 2 

- 

70 0*V3n 

no 

EX 

— 

0.278950 

03 

MU 

- 

rs. 2 00 COD 

01 

Tl 

= 

T.6250CD 

no 

T? 

- 

C.50DDCD 

no 

EX 

= 

0 .248940 

03 

MU, 

- 

r.zcr-rcD 

: i 

Tl 

= 

V. 625000 

00 

T 2 

= 

n.A2 c 00C 

on 

EX 

= 

C.28641D 

C;3 

MU 

"T 

C.2300OD 

•;u 

Tl 

= 

r.6250PD 

oo 

T 2 

= 

3.79» , r-00 

or. 

FX 


0.314910 

03 

. __MU. 

= 

0. 20CQ0C 

_Cil 

Tl 

= 

c;. 70000c 

nn 

T 2 

= 

rt .50C0CD 

00 

EX 

= 

0.27215D 

C*3 

MU 

= 

". 200 

01 

Tl 

= 

0. 7PD0CD 

00 

T2 

- 

0.6? 500 D 

00 

EX 


0 .312240 

03 

MU 

- 

U.2CC'^C 

?l 

Tl 

- 

f .7CD1CD 

■m 

T 2 

- 

r.7'30''C0 

J22 

EX 

a 

0.34168D 

JL2 



TO = 0.0 


TF1 = O.6250CD 00 TF2 = 0.625«0D 30 N =17 MU = C.290QQD 01 


KJ 

O 


TSUR < 1) = 0 . 0 X — a.. - 0 . 16 . 53 90 02 

T SU8 ( 2) = 0.625CCD 00 X = C. 261480 9? 

TSUB 1 3) = n . l25COD 01 X = 0.32 3350 02 

TSUB ( 4) = 0.187500 Cl X = 0. 8^4280 02 

■ T SUB < 5) = 0 . 2 SCOOP 01 X-=- £-..792-020 02 

TSUR ( 6) = 0.212 500 f.l X = 0.113530 03 

T S U R ( 7 1 a 0_2.7-.5S0n Cl X— =--0.1 36 240 7 3 

TSUB ( 8) = 0.437500 Cl X = 0.155600 03 

T S U R ( 9) = 0.500000 01 X = 0.167600 0 3 

TSU8U0) = 0.562E0D Cl X = 0.17592D T3 

TSUB( 11) = 0.625000 Cil X = 0.179350 (?3 

TSUB ( 12 > = ^ .6 8750D 01 X = 0.192650 03 

TSURM3) = t .750000 m x = 0.1R7440 03 

TSUB ( 14 ) = 0.812500 Cl X = 0.200260 03 

TSURM5) = 0.87 5000 M X = 0.190810 03 

TSUB ( 16 ) = 0 .9375 CD 01 X = 0.208350 03 

TSUB < 17 ) = O.IOISOD f'2 X = ". 219050 0 3 


FOR 

Fig. 

6 - a 










INPUT-CAfi CS RE A D 14.44.31.29 3/03/70 

DATA* NS 41 TS l.ODS NMU 2 ICLASS 3 N0PT1 2 ISTOP C Ml 1 M2 1 *DATA 

-DATA* ZETAB, 0.800 WB C.2PO 1WIM 1 K G A 1 N l,.:\BQ.TSii .D.KC.O EBY * Q A T, A 

DATA* HBY 1.5D-2 YBK C.CD? 0. "DO C.C2DO P.C6CC 3.120'? C.195CG C.29C^ *DAT A 

.. DAIAi n .39Bnn r.Binr 0^41*02 tUMJK. C. flnr L^&6Sm <L*23m LUSt&na 0 1 33 00 *nflTA 

DATA* YBK{ 17-41) 1.C0C MU l.CDC 2. IDO NTCH 1 T1CH 2.30C T20 4. >DO *DATA 

—HAUL* lnN f UA-Ed. .ZE.T.AN il. 3 D C > H A 1.0 CO *DATA 

DATA* ZFTAA O.8D0 YO 1 .CDC T1 2. <?CO T2 4. ODD *DATA 

— D-AXA&l *nA T A 


TO LIBS IN THE INTERVAL ( TST.TMAX) 


TV = Q.53CDCD Of; 


MU = O.IOCCOD 01 T 1 = C. 20.0000 01 T2 = '.40i?<:OD 01 EX = 0.562540 CO 


MU « 0 . 200 ODD Cl T 1 = C.20C0CD 01 T2 = 0.400000 01 EX = 0.236190 01 



TC = 0.0 


TF1 = 0.200C0D 01 TF2 =0.400000 01 N =11 


OTLEFT =-0 . 1CCC0D 01 DTRIGHT = 0.224410 01 


MV = 1.6 


T SUB ( 1) = O.C X = G.41162D 00 

TSUB( 2) = 0. 400000 Cl X =__£_• 70 J. 4.6 C _00 

TSUB ( 3) = 0.80 C COO Cl X = 0.976520 0,0 

T SUB ( 4) = C.12CC0D 02 X = 3.116670 01 

TSUB ( 5) = 0 . 160COD 02 X = 0.12158C 01 

TSUB ( 6 ) = 0 . 2QCC00 C2 X = 0.110 570 SI _ 

TSUB ( 7 > = 0 , 240C0D 02 X = 0.864340 00 

TSUB ( 8_1_=_ 0.2 80 COD (?. X = 0. 556250 SO 

TSUB ( 9) = 0.32CCOD 02 X = 0.259450 00 

TSUB (10) = 3.360000 0? X = 0.384330-01 

TSUB ( 11 ) = 0.4OCC00 C2 X =-0.213030 00 


FCRCEC SOLUTION 


T S U B ( 1) = O.C XFREE = 0.589380 00 
TSUB ( 2) = S .AfiC CCiD 01 . _X_FRE_E « 58960 0-0 3 
TS UB ( 3) = C.ecerr.D Cl XFREE =-C. 1C131D-02 
TSURI 4) = C.12CC0D 02 XFREE = 0.59922D-O4 


._j-.SUB-( 5) =-C,4-60.T:G0-.02- XFREE - = -0 . 191900-7 5- 

TSUB( 6) = C.20nCCD 02 XFREE = 0. 157960-17 

— 7 SUB4 - 7-) — =— £ . 2 4 C C C D £2 XFR.FF ■= 0.223890-. 18 

TSUB ( 8) = 0 .2 8000-0 02 XFREF =— : . 16C84D- 19 

-T 5UB< ■ 9) . a..32- C - 000 -0;2 X F - RE - E - = -04-5-9 4880-11 

TSUB ( 10 ) = 0.360000 f2 XFREE =-0.903750-13 

- TSUB< UJ. -= rU-4CtC^O...*2 XF-REE- =-i?. 44 514D-14 


TSLBl 1) 
TSUB ( 2) 

TSUB I. .3) 

TS UB ( 4) 
T SUB ( 5) 
TSUB ( 6) 
TSUB ( 7) 
TSUB ( 8) 
TSUB ( 9J 
TSUB(IO) 
TSUBlll) 


— r* r 

= r.ATCfCD 01 

—rr — 0.-* Q&Z&&D-LL-. 

= c.izrcro t? 

= C.lfctCCD -12 

= c.2ir^o r*2 

= .G.24CGCD LZ 
= C:.2 8CC^O ?2 
= £+32££l}£}..J22.. 
= 0.3 6f f CCD 02 
= C.VjCVZD 02. 


X.SUM. 
X SUM 
X SU M 
XSUM 
.XSUM 
XSUM 
XSUM 
XSUM 

XSUM 

XSUM. 


JD^IOGSOD. 01 
O.7C205D CO 
,0.. S-75-5JLD 


0.1166 80 r< 1 
-0 * 1.2 1.5.8 Q Cl 


-F.-0.K ./TJ 6. . 


0.11C57D “1 
n.8.fc434C..Ca 


0.55625D nn 


G . 3 ffA 3 3D— r 1 
£ ■ 213C3.Q. 03 
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TO = 0.0 


TF1 


0 . 20000 D H 1 


TF2 = O.4GO0OD 01 


TSU B ( 1) 
T SUB ( 2) 
TSUB ( 3) 
TSU8 ( 4) 
TSUB ( 5) 
TSUBi 6) 
TSUB ( 7) 
T SUB f 8 J . 
T SUB ( 9) 
TSUB(IO) 

Tsusni) 


T S U 8 ( 1 ) 

TSUB ( 2) 
T S U B ( 3) 
T S U B ( 4) 
T S U B ( 5 ) 

T S UB ( .6 ) 
T SUB ( 7) 
...TSUB .! 81 
T S U B ( 9) 
TSUB (1C) 
T S UR < 11 ) 


-TSUBt X 
T S U B { 2) 

- TS U B ( 3) 
T SUB C 4) 
T SUB { 5 ) 
T S U B ( 6) 
..I SUE I 71. 
TSUB ( 8) 
T SUB ( 9) 
TSUB(IO) 
TSUB( 11 ) 


A 

1 


N =11 


DTLEFT =-0.100000 01 CTRIGHT = 0.324410 01 


M V - 1.0 


= 0.0 

= P»40Ci00D Cl 
= 0.80000D Cl 
= 0.12CCCD 02 
= 0.160000 C2 
= 0.200000 02 
= 0. 24G0CD C2 
= 0.28000 0 G ? 
= C.32000D C2 
= 0,360000 C? 
= 0.400 COD 02 


X = 0.44209C 00 
X = 0 .9 0 9980 PC 
X = C. 113740 01 
X = 0. 98710 C 00 
X = 0.85247D 00 
X = 0.10410D 01 
X = 9.12234D 01 
X = Q.89135D O P 
X = 0.17760D 00 
X =— 0 .265930 00 
X =-0.918820 00 


FORCED SOLUTICN 


= 

0 .0 


XFREE 

= 0.55791D r>Q 



= 

0 . COD 

Cl 

XFREE 

=-0.165580-02 



- 

C • 800000 

01 

XFREE 

=-0.827460-03 



= 

0.12CCC0 

c? 

XFREE 

= 0.524950-04 



= 

c . lefcoo 

02 

XFREE 

=-Q.178f 9D-:5 



= _ 

r.2vc COP 

C2 

XFREE 

= 0. 198350-07 




n. 24COOD 

C2 

XFREF 

= 0.17666D-O8 



- 

n,28C CCD 

.02 


=-0.139160-^9 



= 

0 . 3 2 C COO 

r? 

XFREE 

= 0.543020-11 



= 

i .secern 

C2 

XFREF 

=-C. 9 5 22 30-1 3 




r .40.0CCD 

0? 

XFREE 

=-C. 32983D-14 




XSUM = 




C.4CCCCD 01 
c.acccoa ci 
0 . 120 C C D 02 
_C* 1 60 CCD . .112 


XSUM = 0 • 9 C 8320 00 
XSUM = 0.1 1366D 01 
XSUM = C.98716D r*fi 
X SUM . = Q .852470 13 0 


- 

l-.2C-:CCD 

c? 

XSUM 

= 0.104100 01 

- 

c t 2 4 n r n p 

0? 

XSUM 

= 0.12? un n\ 


C .28CCOD 

C2 

XSUM 

= 0.891350 CO 

= 

O.32CC0D 

C2 

XSUM 

= G.177.6CD. on . 

= 

C.36CC0D 

'2 

XSUM 

= -0.265930 ,v .) 

= 

0.40CC00 

.C‘2 . 

XSUM 

=“0.918820 .DO 


FOR F I 6 . <3 *" J i 
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TC 


0 . 0 


TF1 = O.2O0n0D 01 


TF2 = C. 400000 01 


N =11 


... CTLEFT =-C.lCCQOD 01 DTRJGHT r= O^Al 471 D 00 


H U - 1.0 


T SUB ( 1) = 
TSUB( 2) . = 
T S U 0 ( 3) = 
TSUBC 4) = 
T S U B { 5) = 
T SUB ( 6 ) = 

T SUB ( 7) = 
. T.S..UB (. .8 ) = 
T SUB ( 9) = 
.TSU,B(.10)..5 
TSUB(ll) = 


0.0 

0.4CCCCD Cl 

C.ercCfiO Cl 

0.12Cft0D C2 
C.16Cr00 C2 
3.2OC0OD C2 
ft.24CftnD C 2 
0.2 800.00 02 
0.2 20000 C2 

0.4CCC0D 02 


X = 0.47119D 00 
X = (1.95238C 00 
X = 0.101860 01 
X = C.98256D Of; 
X = ft . 1 0 3 7 8 D 01 
X_= 0.96 29 7D 00 
X = 0.99853D 00 
X = 0.10CG1D 01 
X = c. 999990 CO 
X = 0.10000D 01 

x = n.iooocD oi 


FORCED SCLUTICN 


T SUB ( 

1) 

= 

0.0 


TSUB( 

-2 ) 

= 

C .400 COD 

01 

TSUB ( 

3 ) 


0.80^ 000 

Cl 

T_SUB.(. 

A). 

= 

0..1.2CC0.D. 

C.2 

TSUB ( 

5) 

= 

C .16CGCD 

C 2 

TSU3( 

6 ) 

= 

0 . 20 f 'C00 

\2 

TSUBC 

7 ) 

= 

0 . 2 A r CCD 

12 

T SUB ( 

8) 

= 

C.2 9CC'"D_ 

02 

T SUB ( 

9 ) 

= 

U • 32CCCD 

n 2 

TJS.UB(12J. 


C • 3 6 C C ft D . 

t!2__ 

TSUB( n ) 

= 

ft.AOCCC D 

T2 


XFREE =-C’« A7119P On 
XFREE =-0. E6553D-02 
XFREE = 0.1 3032D-02 
X FREE =— C . 63 96 3D— 54 
XFREE = C.16798D-05 
XFREE = ft . 5_296 7D-0 8 
XFREE =-0.310950-^8 
XFREE = 0.178130-09 
XFREE =-C. 55417D-11 
XFREE = ft « 371 71D-1 3 
XFREE = n . 69 7 3 2D-1 A 


TS 4 JB( ~ 1 ) 

= - 

- 

. 

-X-SUM 

= 

0.0 - - 



TSUB ( 2 ) 

- 

ft . ACCCOD 

01 

X SUM 

= 

0.943720 

00 

— IS U at — 34 

— 


Ol 


_ 

n. i m QQn 

f*k 1 

TSUB( 4 ) 

- 

C. 12 CCCD 

C 2 

XSUM 

- 

0 . 9825 CD 

00 

- TSUfi ( 5) 

-=- 

c . L 6 CGC 0 

C 2 — 

- XSUM 


_ 0 ^_lJt 37 aD — 0 1 

TSUB ( 6 ) 

= 

G . 20 CCCD 

02 

X SUM 

= 

0 • 962970 

00 

_ TSUB( - 7 ) 

_ = 

U. 24 CC 60 

iI 2 — 

X SUM 

- =. 

-G«S -9 853 D 

tin 

T SUB ( 8 ) 

= 

0 . 28 CC 0 D 

C 2 

X SUM 

= 

0 • 10 C 0 ID 

01 

— I-SUB.t_. 9 J. 

= 



JLSd IM 

- 

n.QQCjqqn 

n 

TSUB ( 10 ) 


0 . 36 C 00 D 

02 

X SUM 

- 

0 . 1 0000 D 

01 

TSUB< 11 ) 

= 

C.AGCCLD 

._£ 2 _. 

X 

_= 

o.iatsuoD 

Cl 


FOR F-l <S . <Z-l 2 
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7. CONCLUSIONS 


The results obtained in the present study indicate that it is indeed possible to 
construct successful models of the transducers and telemetry elements included 
in this report. Furthermore, it is possible to compute most of the necessary 
coefficients by performing some relatively simple experiments, some of which 
may destroy the instruments (as happened with the opened transducers here). 
However, if one models a large series of similar instruments, destructive 
testing of one or a few is not objectionable. The determination of some other 
coefficients required more elaborate tests. In all cases, it is desirable to test 
the model obtained by means of more complex experiments (as was done here) 
to verify assumptions and the range of validity. For the individual categories the 
following more detailed conclusions can be drawn: 

A. Rate Gyros 

The equations of the complete model, in all generality, are given by: 

Mechanical part of gyro: equations (2-6\ (2-7), (2-8), (2-5) 

Electrical output circuit: equations (2-9) through (2-13) 

Bracket: equation (2-14) 

The corresponding numerical values of all parameters are listed in Sections 
2.2.4 and 2. 1.2. 

However, the response of the entire gyro system to the principal input (angular 
rate about the input axis) is very well described by the simple second order 
linar system equation (2-6), where only the first term on the right hand side is 
retained, i. e. , the mechanical response of the gimbal along the output axis is 
the dominant element in the gyro model. Therefore, a considerable economy is 
achieved in describing the gyro behavior mainly in terms of only three parameters: 

The steady state gain (volts per deg/ sec) 

The natural frequency Wgy/ ( 2 tt) 

The damping factor > 

* & y 

The correlation between actual tests in the laboratory and simulations on the 
analog computer using the best-fit parameter values in the theoretical model 
was good, and this for a wide variety of signals being applied to the gyroscopes. 

The overall accuracy of the three main parameters is of the order of 5%, while 
it is somewhat lower (10% to 20%) for the secondary parameters. 


B. 


Accelerometers 


The equations of the complete model are given by 

Accelerometer only: Equations (3-2), (3-3) 

Bracket: Equation (3-4) 

The corresponding numerical values of all parameters are listed in Sections 
3. 1. 1, 3. 1. 2, 3.2. 1. 

The response of the entire system to the principal input (linear acceleration 
along the sensitive axis) is very well described by a simple second order system 
(the first term is equation (3-2)), with the same three parameters as for the 
gyros: 


Steady state gain 
Natural frequency 
Damping factor 

The correlation between tests on the transducers and simulations based on the 
theoretical model was good. The overall accuracy for the three main parameters 
was of the order of 15%, mainly due to spread among individual instruments of 
the same category. (See Section 3.2. 1), An important, and surprising, result 
was: 

a) The low value of the damping coefficient ^ n (see again section 3. 2. 1), 
as opposed to the higher values obtained by the manufacturer for slower, 
more "normal" signals. As a consequence, the response of the 
accelerometers exhibits a marked overshoot for very fast transient input 
signals. 

b) A non-linear effect was observed, where the damping factor ^ 
depends on the speed of variation of the input signals, with f decreasing 
for faster inputs (due to various non-linear effects in the internal 
structure, such as eddy-currents, etc). 

C. Pressure Transducers 

The complete theoretical itd del is described by: 

Equation (4-2) for the Bourdon tubes 

Equations (4-3) and (4-4) for the T-tubing (Bourns transducer) 
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The corresponding numerical values of the parameters are listed in Sections 
4. 1. 1 (Bourdon tubes), 4. 1. 3 (T -tubing). The agreement is very good 
between experimental data and theoretical predictions for the Bourdon tubes; 
this is quite satisfactory as no other resonance calculations were found in a 
literature search for the Bourdon tubes and the present treatment was a first 
attempt for this problem. 

For practical purposes, a sufficient model was found by cascading two second 
order systems, equation (4-6), depicting Bourdon tube and orifice effect, or 
else by an added second order numerator, equation (4-7), to include the effect 
of the T-tubing. Therefore, modeling is mainly reduced to finding 

A steady- state gain 

Two or three natural frequencies 

Two or three damping factors 

Because of interaction, these models are to be taken as a whole, and should not 
be split (except for purely formal operations). 

As for the preceding transducers, the correlation between tests in the laboratory 
on instruments and analog computer simulations, using best fit parameter 
values was good. The overall accuracy of the parameter values obtained was 
of the order of 10% to 20%, again mainly due to spread among individual 
instruments. As for the accelerometers, a surprising resuLt was the low damping 
for fast input signals, leading to a considerable overshoot in the responses. 

This contrasts with the much higher damping obtained by the manufacturer for 
slower, more regular signals. There was also some non-linear effect in 
damping (now, with damping somewhat increasing for extremely fast inputs) due 
to non-linear dissipation mechanics in the turbulent fluid flow. 

D. Telemetry 

A first observation, substantiated in section 5. 1, is that some of the trans- 
ducers have a bandwidth larger than that of the channel assigned to them (under 
standard IRIG conditions); this holds typically for the Giannini transducer in 
channel 11 and the roll rate gyro in channel 7. 

The following conclusions can be drawn for the present airborne equipment. 

The phase-sensitive demodulators, used in conjunction with the rate gyros, are 
modeled by a second order system, equation (5-1), whose damping factor depends 
on the polarity of the output voltage swing, due to a biasing by means of resistors. 



The demodulators introduce a dynamic lag of roughly the same magnitude as the 
gyros themselves; and, therefore, their presence is important. There was 
also a noticeable spread in damping among individual devices. 

The subcarrier oscillators are modeled by a linear second order system, 
equations (5-2) and (5-3). For practical purposes, the lag due to SCO's is 
small compared with the lag of other components and may be neglected in many 
cases. 

The mixer can be considered as an ideal transducer without any lag. The same 
remark applies to the transmitter. 

For the ground equipment used in this study, the following conclusions hold: 

The receiver can be considered as an ideal element without any lag. The 
300 kHz and 1MHz tape recorders introduce no discernible distortion for the 
channels studied here under controlled laboratory conditions. 

The bandpass filter is modeled as a critically damped linear second order system, 
which is the same as a cascade of two identical first order systems (equations 
(5-5), (5-6), and (5-7)); for practical purposes their dynamic lag is small compared 
to that of other components when standard filters are used. Further analysis 
showed a pronounced non-linear dependence of the output on the magnitude of 
the input, which c an be approximately described as a lowering of the damping 
factor (from unity) when the magnitude of the input swing increases. However, 
with standard filters, this theoretically interesting non-linear effect is 
rather unimportant in practice since it^nearly entirely suppressed by the low 
pass filters because of the relative speed of the bandpass response. 

The discriminator, in the narrow sense, is modeled by the linear second order 
model, equation (5-10); however, its response is so fast that, for practical 
purposes, the discriminator lag is entirely negligible. 

The low pass filters can be modeled as linear third order systems, equations 
(5-11) and (5-12), of two types (constant amplitude or constant d elay), each 
completely described by a sin gle. time constant. These low pass filters are the 
slowest and hence most important components in the ground equipment; their 
single parameter is determined numerically with excellent accuracy (better 
than 5%). 

The entire telemetry chain can be modeled simply by cascading individual com- 
ponents, if standard filters are used, because there is no interaction then. The 
agreement between tests and analog computer simulations is very good. 
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If very wideband, nonstandard bandpass and iow pass filters are used, faster 
response is obtained; but interference from adjacent channels appears, as 
discussed in section 5. 11 with a theoretical explanation. 

Power interruptions can create traces resembling useful signals; therefore, 
one has to monitor different channels as a help in locating power failures. 

In what preceeds the tolerance on parameter values varied from a few % to 20%. 
This accuracy reflects the combined effects of: 

a) Spread among individual instruments of the same class. 

b) Limit of accuracy of experiments (the generation of precisely 
known mechanical transients -acceleration and angular rate - is not easy). 

c) Inherent limitations in system identification (theoretical limitations). 

Note that in several cases, a shift of say 10% in parameter values produces a 
very much smaller change in output signal. This is a desirable feature for 
instruments since it obviates the need for fine tuning of components. 

E. Inversion Program 

The inversion method developed here gives good results when applied to the 
three kinds of transducers. It appears desirable sometime in the future to add 
to the present computer program a simple digital spectrum analyzer in order to 
remove the necessity to have the user estimate signal bandwidth. The present 
scheme gives the best obtainable results in a minimax sense. It includes an 
error analysis and can be used for studies of parameter sensitivity, Monte Carlo 
simulations of random noise effects, etc. In terms of the inversion capabilities 
in general, one can consider four categories of signals: 

a) Signals that are extremely fast, but without large magnitudes (having 
vanishingly small products of magnitude times duration) will not appear 
at the output of the transducers and hence will pass unnoticed. 

b) Extremely fast signals with large amplitudes (similar to mathematical 
impulses) are observable; exact input wave form reconstruction is 
impossible (as it is with any method) but also unnecessary, since such 
inputs act as mathematical impulses whose strength can be easily 
determined by the measurement of the area under the output curve. 
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c) For other fast signals (slower than the prece-ding categories), 
the inversion program is fully applicable. The results obtained were 
very satisfactory after determination of the proper signal bandwidth. 

Of course, perfect reconstruction is impossible, but under the conditions 
given here, best possible reconstruction is achieved. 

d) Finally, really slow ("normal") signals pose no problem since 
the transducers transmit them with little or no distortion. 

It is interesting to note that for a real system, where noise is unavoidable, a 
relatively low order, simplified model may give better results in the inversion 
process than a completely detailed, high order model, which would only be 
useful for an idealized zero-noise inversion. In fact, the optimum order of the 
inversion model can be determined in terms of the signal-to-noise ratios for a 
given system. 


Avco Corporation 
Wilmington, Massachusetts 

June 30, 1970 
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APPENDIX A 


FREQUENCY TRANSFER -FUNCTION OF A BANDPASS FILTER 


Consider a BP filter, with input signal e^ n (t) ( as defined in section 5. 7) where 
Ak is considered constant. (This is an approximation in itself, as seen in the 
discussion and illustrations of the SCO). Then, the output is 
t V 


e (t) = A 
out' 7 


l f BP ( 4 ' 4 ') cos ( w o 4 ’ + j A w in (t") dt")dt' 


= ZA [ f DLP ( 4 “ 4 ') cos < w o cos K 4 ' + A w in( 4 ' ')dt' 1 )dt' 


- 

This can be rewritten as 
t 

e out^ = A ( f DLp( t-t ') cos < w o 4 - 


cO 


o t -'B + j A 


w. (t")dt")dt' 
in 


t t* 

+ A / fnT.P ( 4 - 4 ') cos ( w 0 ( 4 * 24 ') - < E" j A w in (t")dt")dt' 

o 


/ 


DLP 

CO 


Since \A w^ n \ w Q for the telemetry channels, the second term in the preceding 
formula is nearly zero for a narrow band BP (see explanation related to the 
"first property" in section 5. 7) since the high frequency oscillations 2w 0 t' 
approximately cancel contributions from (t-t 1 ). 

Therefore, t 

e out (t) ~ A Be ex P (i( w o t_ ? )) j 4 DLP ( 4-4 ') ex P ^ j ^ w in (t")dt")dt' 


C^O 


t- IX 


A Re exp (i(w Q t -<g)) J f DL p( “*) exp (i [ A w in (t M )dt")d il 

o 


~ A Re exp (i(w t -<£)) M(t) exp (i 0(t) ) 


which also defines the real functions M(t) (amplitude) and j?" (t) (phase). 
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For A w in = o, i. e. , no disturbance in frequency, one has 

oo 

M(t) exp (ijT(t)) = y^f DLp ( u. ) d u. = F DLP (°) 
o 

(for a physically realizable filter) 
giving ( M (t) = Fqlp (o) 

1 y(t) = o 

This is the expected steady state response to a sine wave. A cos (w Q t 
Now, for A w in (t) £ o t - u- 

OO 

— [ M(t) exp (i (t))] = i f f DLP (^) exp (i ( A w. (t")dt") 

dt / 

o o 

A w. (t - n) d UL 
in ' ' 

Now, provided j A w in ) a ^- so satisfies the more stringent condition 

< I A w in | average ) « edge of D. L. P. 

or ( I 4 w. I J// bandwidth of BP filter, 

v I in I average 

one has 

(M(t) exp (i (t)) = M exp (i £ ) i + exp (i <§ ) 

OO 

= -^ + F DLP (o > 1 if = 1 f f DLP ( " > A w in^ t_ ^ ) d ^ 

o 

or dM = 0 

dt 

CTO 

^ =_i /f nT p(** ) 4 w. (t-u-)du. 

dt f dlp<°> / 

o 

which proves the M Second Property" of Section 5. 7. 
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If l A wi n | approaches the band edge of BP for an extended time, this is not 
longer valid. This fact is clearly illustrated, for the case of an ideal bandpass 
filter and a symmetric step in w. n (t), in the above mentioned reference (in a 
figure shown there, taken from Malinger 1 s article in Proc. I. R.E. 30, 378-S&3, 
August 1942). Therefore, for frequency swingsjA w^ n ] avera g e which are 
comparable with the bandwidth of the BP, pronounced non-linear effects appear. 
They can be calculated numerically for each particular case, by the methods 
used in this appendix if the first formula of this appendix is employed and all 
subsequent simplifications based on the smallness of \ A w * n \are omitted. An 
example of such a calculation is given in section 5. 7. 
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APPENDIX B 


GENERAL DESCRIPTION OF TESTS ON TRANSDUCERS 


Open Instrument Testing; See Figures B1 through B4 on Pages 247 and 248 for 

pictures. 

The three types of instruments were opened to further evaluate their physical 
dimensions, electrical circuits and undamped natural frequencies. These 
dimensions were compared with available vendor information and where necessary 
vendor information was not available the model was updated from actual 
dimensions. The electrical circuit for the accelerometers was also traced out. 

Accelerometer: 


The cover of the accelerometer was machined off. The circuit was traced from 
the two control circuit boards. The power supply board was not traced. 

The output of the signal generator was measured as a function of the angle of 
the pendulum. The high frequencies used in the signal generator precluded the 
use of a scope probe and so the torquer current was measured as a function of 
pendulum displacement using a dummy load. 

The torque of the torque generator was also determined as a function of input 
current. 

Pressure Transducer: 


A small clean hole was machined in the pressure transducer housing. The moving 
element was moved away from its null point and then the restraint was removed. 
The moving member then moved to its normal position. This mechanical step 
function showed the response of the pressure transducer very well. 

The dampening fluid was then removed and the process was repeated. The 
undamped material frequency was clearly indicated. 

The pressure transducer was then disassembled and the dimensions and weights 
of the moving elements were determined. These figures were compared with 
the vendor’s data where available and the inertias were compared to refine the 
mathematical model. 

Gyro: 

A small clean hole was machined in the gyro housing. The entire gimbal 
assembly was moved away from its null point and then the restraint was removed. 
The gimbal then returned to its null. This mechanical step function clearly 
showed the response characteristics of the gyro. 
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The dampening fluid was then removed and the process was repeated. The 
undamped natural frequency was clearly indicated. 

The gyro was then disassembled and dimensions of the various elements were 
recorded to further refine the mathematical model with the actual dimensions. 

Carco Table Testing: 

The Carco table is a three -axis flight simulation manufactured by Carco 
Electronics. It is a precision electrohydraulic angular positioner. (See 
Figure B-5). 

Because of the limits of the response of the table to a given position input, a 
special velocity transducer (tachometer) was incorporated into the table. This 
device had a frequency response well above the response being investigated and 
therefore was suitable to compare to the response of the instrument under test. 

A ramp, blast and square wave function were fed into the position input of 
the simulator. A sine input was also fed into the system with the function 
superimposed. The output of the velocity transducer was the input to the instrument 
being analyzed. 

For the accelerometer the tangential motion of the table was used to simulate 
the acceleration. The output of the velocity transducer was differentiated to 
compare the input of the accelerometer with the output to the accelerometer on 
the analog computer. 

Over Range Testing: 

The three types of Scout instruments were subjected to inputs that were higher 
in magnitude than the limits of the instrument but of durations that were a fraction 
of a time constant for the response. To accomplish these inputs special fixtures 
were made. The boundaries of the systems were varied to include mounting 
brackets and plumbing configurations similar to the Scout vehicle. 

Pressure Transducer: (See Figures B-6 and B-?) 

Two types of inputs were used for the pressure transducer. A shock tube was 
used to generate very fast rise inputs. The input would decay, however, as 
the shock wave reflected up and down the tube. The applied pressure was 
monitored with a fast response transducer (piezo-electric) and this output was 
compared with the output of the Scout instrument. 

To vary the input rise time a second fixture was fabricated. This fixture 
consisted of a large volume, a burst disc, a variable orifice, and a small 
chamber to which the Scout pressure transducer and the reference transducer 
were connected. 
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figure b-6 


PRESSURE SHOCK TUBE 
WITH HEADCAP 
SIMULATOR 


FIGURE B-7 

SMALL PRESSURE 
FIXTURE WITH 
SIMULATED SCOUT 
TUBING 




The rise time was varied by initiating the burst time of the diaphragm and 
then later by restricting the flow into the small chamber. 

When it was applicab le the plumbing was simulated and then tests were con- 
ducted on close coupled instruments. The pressure path of the motor endup was 
also simulated. 

Gyro: (See Figure B-8) 

The Scout gyro was tested on a rotating fixture that was specially constructed 
for the purpose. The fixture also contained the gyro mounting block and the 
mounting block to vehicle bracket. 

Impulses were applied to the moving member with various size and material 
hammers. The table was also stopped with variable deceleration rates. The 
deceleration was accomplished by driving a shaped wedge into various materials 
such as bees wax, paraffin, and polyvinyl slugs. 

The inputs to the system were monitored by moving through a known angle in 
a period of time that was monitored by a piezo-electric accelerometer. 

Accelerometer: (See Figure B-9) 

The Scout accelerometers were tested on the fixture that was fabricated for 
testing the gyros modified to accept the accelerometers and their mounting 
bracket. They were mounted to sense the tangential velocity of the rotating arm. 

Various sizes and material hammers were used to provide an impulse of various 
rise times and amplitudes. Different types of stops were used to decelerate 
the moving element. 

The input to the system was monitored by a fast response accelerometer and 
the output of the Scout accelerometer was compared to this input. 


231 


FIGURE B-8 


FIXTURE FOR TESTING 
GYRO AND 
ACCELEROMETER 
(WITH GYRO) 


FIGURE B-9 

FIXTURE FOR TESTING 
GYRO AND 
ACCELEROMETER 
(WITH ACCELEROMETER) 



APPENDIX C 


EXPERIMENTAL TEST SETUPS FOR THE SCOUT T/M INVESTIGATION 


This appendix describes the various test setups used to gather the experimental 
information needed for the analytical descriptions of the Scout Telemetry Shelf, 
(Tele Dynamics Model 1799A S/N 101). 

The test setups were based on the T/M shelf documentation received from 
NASA/Langley with the exception of the mixer input circuit. The actual shelf 
has a 6. 98 IK resistor from the signal input to signal ground and this was 
incorporated in the test. Each test is described separately in the following pages. 

Subcarrier Oscillator Steady-State Test 

Each subcarrier oscillator (SCO) tested was removed from the shelf and tested 
individually as shown in Figure C-l. The D 0 C. standard supplied the input 
voltage which was checked with a digital voltmeter. The output frequency of 
the SCO was determined by a counter and monitored on an oscilloscope. The 
input voltage was varied between 0. 0V and 5. 0V except for some SCO where the 
voltage was varied between -15V and 30V to determine the out-of-band 
characteristics of the SCO. The data recorded was the d. c. input voltage with 
the corresponding output frequency and voltage. 

Subcarrier Oscillator Transient Response Test 

The transient response of an individual SCO was determined by modulating the 
SCO (See Figure C-2) with a square wave (OV to 5V) of a sufficiently long period 
to allow for the SCO settling time. The square wave was approximately changed 
for each SCO tested. The counter was used to verify the bandedge frequencies. 
The actual data was taken on a Polaroid photograph displaying the modulating 
frequency and the SCO output waveform. The delay trigger function of the 
oscilloscope was used to look at the rise and fall times of the input and output 
waveforms. The Plug-In was operated in its chopped mode. 

Mixer Tests 


The same test setup was used for all three mixer tests (See Figure C-3). The 
first test determined the gain of the mixer as a function of frequency. This 
information was obtained using a sinusoidal input and measuring the R. M. S. input 
and output voltages. The second test determined the gain linearity. The input 
level of a sinusoidal voltage was varied and the corresponding output voltage was 
noted. The gain of the amplifier at the various input voltages was then cal- 
culated. The last test used triangular waveforms at various repetition rates. 

The response of the amplifier to these waveforms was recorded on Polaroid 
film from an oscilloscope display. 
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Item 

DoC. Standard 

Digital Voltmeter 

Power Supply 

Counter 

Oscilloscope 

SCO 0 


Equipment List 


Manufacturer and Model jj= 

Electronic Development 
Corp o , MV 100NK 

Hewlett-Packard 3^0A 
with 3bk^A Plug-In 

Harrison Labs 6204A 

Hewlett-Packard 522B 

Tektronix 535 with Type 
CA Plug-In 

Tele-Dynamics 5001AA-3 


Figure C-i- Sub-Carrier Oscillator (SCO) 
Steady-State Test 
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Item 

Waveform Generator 

Power Supply 

Counter 

Oscilloscope 

SCO 


Equipment List 


Manufacturer and Model # 

Hewlett-Packard 33COA with 
model 3 30 1+A Plug-In 

Harrison Labs 6204A 

Hewlett-Packard ^22B 

Tektronix 535 with 
Type CA Plug-In 

Tele Dynamics 5001AA-3 


Figure C-2, Sub-Carrier OsciiLator 
Transient Response 
Test 
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Equipment List 


Item 

Waveform Generator 

Attenuator 
Voltmeters 
Power Supply 
Mixer 

Oscilloscope 


Manufacturer and Model # 

Hewlett-Packard 3300A 
with model 3304A Plug-In 

Helipot T- 104A 

Hewlett-Packard 400FL 

Harrison labs 620kk 

Tele -Dynamics 119TA 

Tektronix 535 with Type 
CA Plug-In 


Figure G3- Mixer Tests 
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Phase -Sensitive Discriminator, Steady-State Tests 


The steady- state tests on the Phase Sensitive Discriminator (P. S. D. ) were 
performed using the test setup of Figure £-4. A test jig was made up (See 
Figure C-5) to facilitate the P e S. D. testing* The test jig allowed the reference 
voltage phase to be reversed while low-level signals were applied to the P. S. D 0 
The power amplifier was necessary to supply the specified reference voltage 
level which could not be obtained from the waveform generator. 

The voltmeters were used to obtain the input-output voltage relationships; the 
oscilloscope was used to monitor the output voltage. 

Phase -Sensitive Discriminator, Transient Response Test 

The transient response of the P. S. D. used the test setup of Figure C-6 which 
is essentially the same as Figure C-5 except for the Avco Low Frequency 
Generator, The L. F. Generator was used to modulate the signal voltage. A 
schematic of the L. F. Generator is shown in Figure C-7. The generator was 
originally designed to have a high power square wave capability of amplitude 
and V 2 . The response of the P. S. D. was displayed on an oscilloscope and 
recorded on Polaroid film. 


System Shelf Test - Configuration 1 

For this test the whole shelf was powered but only one channel was modulated. 
The modulation waveforms were square wave, triangular and trapezoidal. The 
latter was supplied by using a saturating transistor amplifier, see Figure C-9, 
between the waveform generator and the SCO. The shelf output, see Figure C-8, 
was sent to discriminator for demodulation. The waveforms were observed on 
an oscilloscope and recorded on Polaroid film. 


Saturating Transistor Amplifier 

The transistor amplifier was as shown in Figure C-9. Since the waveform 
generator used in the tests has signal offset and signal amplitude controls, the 
desired waveform output was a matter of selection the proper levels of the two 
generator contacts. The power supply to the amplifier was set initially at +6V 
and then find adjusted for the proper bias operation. 


Variable 

Input 



Figure C-9 - Saturating Transistor Amplifier 
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Equipment List 


Item 

Waveform Generator 

Avco P.S*Do Test Jig 
Power Amplifier 

Phase Sensitive Discriminator (PoSoD*) 
Power Supply 
Voltmeter, Output 

Voltmeter, Input 
Oscilloscope 


Manufacturer and Model # 

Hewlett-Packard 3300A with 
330bA Plug-In 

Special Test Circuit 

Eico HF-22 

Tele -Dynamics 1353& 

Harrision Lahs 620kA 

Hewlett-Packard 52^5L with 5265A 
Plug-In 

Hewlett-Packard 400FL 

Tektronix 535 with type 
CA Plug-In. 


Figure C4- Phase Sensitive Discriminator 
Steady-State Tests 
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Item 


Equipment List 


Manufacturer and Model # 


Waveform Generator 

Avco Low-Frequency Generator 
Avco P.S.D. Test Jig - 
Power Supply 
Power Amplifier 
Power Supply P.S.D. 

Phase Sensitive Discriminator (P.S.D.) 
Voltmeter, Output 

Voltmeter, Input 


Hewlett-Packard 3300A with 
model 3304 A Plug-In 

Special Test Circuit 

Special Test Circuit 

Hewlett-Packard 721A 

Eico HF-22 

Harrison Labs 6204A 

Tele -Dynamics 135 3A 

Hewlett-Packard 5245L 
with 5265A Plug-In 

Hewlett-Packard 400FL 


Oscilloscope Tektronix 535 with Type 

CA Plug-In 

Figure C6- Phase Sensitive Discriminator 
Transient Response Test 
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+28 V 


Ground 


Sync 

Output 


B2 

B1 



VI 



Item 

Waveform Generator 

Power Supply 
Telemtery Shelf 
Discriminator 

Oscilloscope 


Figure C8- System Shelf Test 
Configuration 1 


Equipment List 

Manufacturer and Model # 

Hewlett-Packard 3300A with model 
3304A Plug-In 

Harrison Labs 6204 a 

Tele -Dynamics 1799A 

Electro-Mechanical Research Inc. 
Model 229 Tunable Discriminator 

Tektronix : 535 with Type CA 
Plug-In 
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System Shelf Test - Configuration 2 


This test was connected as Configuration 1 but with the P. S. D. driving the 
SCO. The P. S. D. test setups of Figure C-4 and C-6 were used to drive the 
P. S, D. The waveforms were observed on an oscilloscope and recorded on film. 

Power Interruption Test 

The test setup of Figure C-8 was used with a triangular input waveform. The 
T/M shelf was connected to the power supply/ supplies through the Avco L. F. 
Generator (See Figure C-10a). The rate and level of power interruption was 
changed with the results displayed on the oscilloscope with a Tektronix 1A4 plug-in. 
This displayed a high and low SCO channel, the type of power interruption and 
the modulating waveform. 

Figure C-lOb shows the method used to obtain a transient (short duration) power 
interruption with the same test setup. 

Discriminator Filter Tests 

Figure C-ll shows the test setup used to determine the response of the 
Discriminator Bandpass and Low pass filters. Figure C-12 shows another 
method that is useful in determining the transient response of the filters. 
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Figure 10a 

PS*/ 



PS#L 

Figure 10b 


Equipment List 

Item Manufacturer and Model # 

Po-wer Supply #1 and #2 Harrison Labs 6204A 

Power Supply #3 Hewlett-Packard 721A 


Figure CIO - Power Interruption Test 
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Equipment List 


Item 


Manufacturer and Model # 


Waveform Generator 
Voltmeter 


Hewlett -Packard 3300k 
with 330bk Plug-In 

Hewlett-Packard 400FL 


Osci lloscope 


Tektronix 535 with Type 
CA Plug-In 


Bandpass Filter 


Part of E'ectro-Mechanical 
Research Inc., Model 
229 Tunable Discriminator 


Figure Cll - Discriminator Filter Test 
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Item 

Waveform Generator 


Equipment List 


Avco Low-Frequency Generator 

Power Supply 

Voltmeter 

Band Pass Filter 


Oscilloscope 


Manufacturer and Model # 

Hewlett -Packard 3300A 
with Model 3304A Plug- IS. 

Special Test Circuit 

Hewlett-Packard 712A 

Hewlett Packard 400FL 

Part of Electro -Mechanical 
Research Inc., Model 229 
Tunable Discriminator 

Tektronix 535 with Type 
CA Plug-In 


Figure C12 - Discriminator Filter Test Transient Test 
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FIGURE B-l RATE GYRO 



FIGURE B-2 ACCELEROMETER 
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